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Abstract 

We present an up-to-date review of Big Bang Nucleosynthesis (BBN). We discuss 
the main improvements which have been achieved in the past two decades on the 
overall theoretical framework, summarize the impact of new experimental results on 
nuclear reaction rates, and critically re-examine the astrophysical determinations of 
light nuclei abundances. We report then on how BBN can be used as a powerful 
test of new physics, constraining a wide range of ideas and theoretical models of 
fundamental interactions beyond the standard model of strong and electroweak 
forces and Einstein's general relativity. 
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1 Introduction 



A remarkable scientific achievement in the second half of the 20**^ century 
has been the establishment of "Standard Models" of Particle Physics (SMPP) 
and Cosmology (SMC). In particular, the latter has been possible thanks to 
an incredibly fast growth of the amount and quality of observations over the 
last couple of decades. The picture revealed is at the same time beautifully 
simple and intriguingly mysterious: on one hand, known gauge interactions and 
Einstein's general relativity seem able to explain a huge wealth of information 
in terms of a few free parameters specifying the composition/initial conditions 
of the Universe; on the other hand, these numbers are not explained in terms 
of dynamical processes involving the known fields and interactions. This is the 
case of the "dark energy" density (consistent with a cosmological constant), 
of the non-baryonic dark matter, of the baryon-antibaryon asymmetry, the 
flatness, homogeneity and isotropy of the universe on large scales, etc. 

The very success of the cosmological laboratory is thus providing much indi- 
rect evidence for physics beyond the SMPP. On the other hand, advances in 
particle physics (a very recent example being the phenomenology of massive 
neutrinos) have an impact at cosmological level. This interplay has proven ex- 
tremely fertile ground for the development of 'astroparticle physics', especially 
since many theories beyond the SMPP predict new phenomena far beyond the 
reach of terrestrial laboratories, but potentially testable in astrophysical and 
cosmological environments. In this respect, the nucleosynthesis taking place in 
the primordial plasma plays a twofold role: it is undoubtedly one of the obser- 
vational pillars of the hot Big Bang model, being indeed known simply as "Big 
Bang Nucleosynthesis" (BBN); at the same time, it provides one of the earliest 
direct cosmological probes nowadays available, constraining the properties of 
the universe when it was a few seconds old, or equivalently at the MeV tem- 
perature scale. Additionally, it is special in that all known interactions play an 
important role: gravity sets the dynamics of the "expanding cauldron" , weak 
interactions determine the neutrino decoupling and the neutron-proton equi- 
librium freeze-out, electromagnetic and nuclear processes regulate the nuclear 
reaction network. 
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The basic framework of the BBN emerged in the dec ade bet ween the seminal 
Alpher-Bethe-Gamow (known as af3'y) paper in 1948 (jAlp48l ) and the essential 
settlement of the paradigm of t he stell ar nucleosynthesis of elements heavier 
than ^Li with the B^FH paper teurSTI ). This pioneering period — an account 
of which can be found in (IKraQGl ) — established the basic picture that sees the 
four light-elements ^H, '^He, ^He and ^Li as products of the early fireball, and 
virtually all the rest produced in stars or as a consequence of stellar explosions. 



In the following decades, the emphasis on the role played by the BBN has 
evolved significantly. In the simplest scenario, the only free parameters in 
primordial nucleosynthesis are the baryon to photon ratio rj (equivalently, the 
baryon density of the universe) and the neutrino asymmetry parameters, 7]^^ 
(see Section 6.3). However, only neutrino asymmetries larger than rj by many 
orders of magnitude have appreciable effects. This is why the simple case 
where all fjujs are assumed to be negligibly small (e.g., of the same order 
of 77) is typically denoted as Standard BBN (SBBN). Since several species of 
'nuclear ashes' form during BBN, SBBN is an over-constrained theory whose 
self-consistency can be checked comparing predictions with two or more light 
nuclide determinations. The agreement of predicted abundances of the light 
elements with their measured abundances (spanning more than nine orders of 
magnitude!) confirmed the credibility of BBN as cosmological probe. At the 
same time, the relatively narrow range of 77 where a consistent picture emerged 
was the first compelling argument in favor of the non-baryonic nature of the 
"dark matter" invoked for astrophysical dynamics. 



The past decade, when for the first time a redundancy of determinations of 
7] has been possible, has stressed BBN as a consistency tool for the SMC. 
Beside BBN, one can infer the density of baryons from the Lyman-a opacity 



in qua s ar spec tra due to intervening high redshift hydrogen clouds (jMei93 



Rau96l : IWei97l ) ; from the baryon fraction in clusters of galaxies, deduced from 



the hot x-ray emission (jEvr97l ): most importantly, from the height of the 
Doppler peak in the an gular po wer spectrum of the cosmic microwave back- 
ground anisotropy (see (jPunOSl ) for the latest WMAP results). These deter- 
minations are not only mutually consistent with each other, but the two most 
accurate ones (from the CMB and BBN) agree within 5-i-10%. While losing to 
CMB the role of "barometer of excellence" , BBN made possible a remarkable 
test of consistency of the whole SMC. 



It comes without surprise that this peculiar 'natural laboratory' has inspired 
many investi gations, as testified by the numerous reviews existing on the sub- 
ject, see e.g. f lMal93l : ICop95l : lOliOOl : ISar96l : ISch98l : ISteOTl l Why then a new re- 
view? In the opinion of the authors, a new BBN review seems worthy because 
at present, given the robustness of the cosmological scenario, the attention of 
the community is moving towards a new approach to the BBN. On one hand, 
one uses it as a precision tool in combination with other cosmological informa- 
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tion to reduce the number of free parameters to extract from multi-parameter 
fits. On the other hand, BBN is an excellent probe to explore the very early 
universe, constraining scenarios beyond the SMPP. The latter motivation is 
particularly intriguing given the perspectives of the forthcoming LHC age to 
shed light on the TeV scale. A new synergy with the Lab is expected to emerge 
in the coming years, continuing a long tradition in this sense. Finally, a wealth 
of data from nuclear astrophysics and neutrino physics have had a significant 
impact on BBN, and it is meaningful to review and assess it. In particular, 
the recent advances in the neutrino sector have made obsolete many exotic 
scenarios popular in the literature still a decade ago and improved numerous 
constraints, providing a clear example of the synergy we look forward to in 
the near future. 

This review is structured as follows: in Section 2 we summarize the main cos- 
mological notions as well as most of the symbols used in the rest of the article. 
Section 3 is devoted to the description of the Standard BBN scenario. Section 

4 treats the status of observations of light nuclei abundances, which in Section 

5 are compared with theoretical predictions. The following Sections deal with 
exotic scenarios: Section 6 with neutrino properties. Section 7 with inhomo- 
geneous models. Section 8 with constraints to fundamental interactions and 
Section 9 with massive particles. In Section 10 we report our conclusions. Al- 
though this article is a review, many analyses have been implemented ex-novo 
and some original results are presented here for the first time. Due to the 
large existing literature and to space limitation, we adopt the criterion to be 
as complete as possible in the post-2000 literature, while referring to previous 
literature only when pertinent to the discussion or when still providing the 
most updated result. Also, we adopt a more pedagogical attitude in introduc- 
ing arguments that have rarely or never entered previous BBN reviews, as for 
example extra dimensions or variation of fundamental constants in Section 
8, while focusing mainly on new results (as opposed to a 'theory review') in 
subjects that have been extensively treated in the past BBN literature (as 
in SUSY models leading to cascade nucleosynthesis, the gravitino 'problem', 
etc.). Other topics, which for observational or theoretical reasons have at- 
tracted far less interest in the past decade in relation to BBN bounds, are 
only briefly mentioned or omitted completely (this is the case of technicolor 
or cosmic strings). Older literature containing a more extensive treatment of 
these topics can be typically retraced from the quoted reviews. In the follow- 
ing, unless otherwise stated, we use natural units h = c = ks = I, although 
conventional units in the astronomical literature (as parsec and multiples of 
it) are occasionally used where convenient for the context. 
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2 Standard Cosmology 



To keep this review self-contained, and fix the notation which we will be using 
in this paper, we summarize here the main aspects of the cosmological model 
which are relevant for our analysis. The standard hot Big Bang model is based 
on three fundamental astronomical observations: the Hubble law, the almost 
perfect black body spectrum of the background photon radiati on, and the ho - 
mogeneity and isotropy of the universe on large scales, see e.g. (jPeeSOl : |Pee93l ) . 
The latter, also known under the spell of Cosmological Principle implies that 
the metric itself should be homogeneous and isotropic, and singles out the 
Friedmann-Lemaitre-Robertson- Walker (FLRW) models. In comoving spheri- 
cal coordinates one has: 



ds^ = g^^dz'^dx" = df - a^{t) 



dr"^ 



+ r\de'^ + sm^ ed(f)'^) 



where a{t) is the cosmic scale-factor and k = 1,0, —1 the rescaled spatial 
curvature signature for an elliptic, euclidean or hyperbolic space, respectively. 



The Einstein field equations relate the energy-momentum tensor of the perfect 
fluid representing the matter-energy content of the universe. 



T, 



(2) 



with the space-time curvature R 



fiupa ) 



R 



Rfiu — QiMU — 8 TT Gjy Ti^ii, + Ag^i, , 



(3) 



where R^^, is the Ricci tensor, i?^,^ = g^"Rpp,av^ R the scalar curvature, 
R = g^'^R^j,, Gn the Newton gravitational constant, and A the cosmologi- 
cal constant. Substituting the FLRW metric (1) in the Einstein's equations 
(3) gives the Friedmann-Lemaitre (FL) equation for the Hubble parameter H, 



SttG 



N 



P- — 

a^ 



(4) 



The equation of state of the fluid filling the universe, P = P{p), specifying 
the pressure as a function of the energy density, along with the covariant 
conservation of the energy momentum tensor (which accounts for the entropy 
conservation if the fluid corresponds to a thermal bath of particle excitations). 



d{pa^ 
da 



-3 Pa' 



(5) 
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allows one to get the evolution of p as function of a, 



PM(xa ^ , (6) 
PR(xa~^ , (7) 
Pa oc const , (8) 

for matter (both baryonic and dark matter, Pb,Pdm ~ 0), radiation (Pr = 
Pr/3), or cosmological constant (Pa = ~Pa), respectively. 

As usual, the present values of radiation, baryon matter, dark matter and 
cosmological constant energy densities will be expressed in terms of the pa- 
rameters fli = p'^/pcr, i = R, B, DM, A, with pcr = SHg/ (SuGm) the cri tical 



density today and Hq = 100 h km s"^ Mpc-\ with h = 0.73l[]:|j3 dYaoOGf ). To 
quantify the baryon density parameter we will also use Ub = Vtsh? and the 
baryon to photon density ratio, r] = ub/ti^. The latter is also proportional to 
the initial baryon-antibaryon asymmetry per comoving volume produced at 
some early stage of the universe evolution. This ratio keeps constant after the 
— e~ annihilation phase taking place at a value of the photon temperature 
T ~ 0.3 MeV (see later). Moreover, at low energy scales there are no baryon 
violating interactions at work, thus the value of t] can be simply related to 
^Ib, see e.g. flSer04bl ) 



r/io = r/-10i° = — — — , (9) 



0.007Fp V To / V 



where Yp stands for ^He mass fraction (see Section 3.1) and Tq the photon 
temperature today. Note that the numerical factor multiplying Yp takes into 
account the effect of the ^He binding energy on the whole energy budget in 
baryonic matter. 

Matter and radiation fluids can be usually described in terms of a bath of 
particle excitations of the corresponding quantum fields. In particular, at high 
temperatures rapid interactions among them ensures thermodynamical equi- 
librium and each particle specie is described by an equilibrium (homogeneous 
and isotropic) phase space distribution function, 

-1 

(10) 



exp ± 1 



where £'j(|p|) = y\p\'^ + rnf is the energy, +/— corresponds to the Fermi- 
Dirac/Bose-Einstein statistics, and pi the chemical potential, which is zero 
for particles which can be emitted or absorbed in any number (like photons). 
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In the comoving frame, the number density, energy density and pressure can 
be expressed as follows 



n^{T) = g.J^^mPlT) , (11) 
p^(T) = gJ^EM) MIpIT) , (12) 



(2 7r)3 

^'.m^S./I^^J^/.dPl.r) . (13) 

where gi is the number of internal degrees of freedom. The BBN takes place 
in the radiation dominated phase, hence non-relativistic particles contribute 
negligibly to the total energy density, which therefore can be conveniently 
written in terms of the photon energy density = vr^T^/lS, 

Pr = g*Y ' f^^^) 
which defines (7*, the total number of relativistic degrees of freedom. 



where the first and second terms are due to all boson and fermion species, 
respectively. The possibility is left in the previous formula of different Tj for 
different species, accounting for pseudo-thermal distributions of decoupled flu- 
ids (like relativistic neutrinos at BBN times). 

Finally, we will also exploit in the following the definition of the entropy 
density, s(T), in terms of the phase space distribution function. For a given 
specie i one has: 

^.(T) = ^ = ^.y^^^^^/.(|p|,T) . (16) 

The total entropy density is conventionally written as 

siT) = gUT) n, = ^ gUT) ' ^^^^ 



where = (2(^(3)/7r^) is the number density of photons and 
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Use of Eq. (5) implies that entropy per comoving volume is a conserved quan- 
tity, s{t)a^ = const. 



3 Big Bang Nucleosynthesis 



3. 1 Overview 



Extrapolating the present universe back in the past, we infer that during 
its early evolution, before the epoch of nucleosynthesis, it was hot and dense 
enough for electrons, positrons, photons, neutrinos and nucleons, as well heav- 
ier nuclei, to be in kinetic and chemical equilibrium due to the high (weak, 
strong and electromagnetic) interaction rates. In particular, the initial val- 
ues of all nuclear densities are set by Nuclear Statistical Equilibrium (NSE). 
NSE implies that they constitute a completely negligible fraction of the to- 
tal baryon density, which is all in the form of free neutrons and protons. As 
expansion proceeds, weak process rates become eventually smaller than the 
expansion rate H at that epoch, hence some particle species can depart from 
thermodynamical equilibrium with the remaining plasma. This is the case of 
neutrinos which only interact via weak processes and freeze out at a tempera- 
ture of about 2-3 MeV. Soon after, at a temperature ~ 0.7 MeV, neutron- 
proton charged- current weak interactions also become too slow to guarantee 
neutron-proton chemical equilibrium. The n/p density ratio departs from its 
equilibrium value and freezes out at the value n/p = exp(— Am/To) ~ 1/7, 
with Am = 1.29 MeV the neutron-proton mass difference, and is then only 
reduced by subsequent neutron decays. At this stage, the photon tempera- 
ture is already below the deuterium binding energy ~ 2.2 MeV, thus one 
would expect sizable amounts of to be formed via n + p + 7 process. 
However, the large photon- nucleon density ratio r/"^, which is of the order of 
10^, delays deuterium synthesis until the photo-dissociation process become 
ineffective (deuterium bottleneck). This takes place at a temperature Tjv such 
that eyi]){B£i/Tjq)rj ~ 1, i.e. T/v ~ 100 keV, which states the condition that 
the high energy tail in the photon distribution with energy larger than Bo 
has been sufficiently diluted by the expansion. 

Once starts forming, a whole nuclear process network sets in, leading to 
heavier nuclei production, until BBN eventually stops, see Section 4. An es- 
timate of the main BBN outcome, i.e. ^He, can be obtained with very simple 
arguments, yet it provides quite an accurate result. Indeed, the final density 
n4Hc of ^He is very weakly sensitive to the whole nuclear network, and a very 
good approximation is to assume that all neutrons whi ch have not d ecayed 
at T/v are eventually bound into helium nuclei, see e.g. flKol90l : ISar96l ). This 
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leads to the famous result for the helium mass fraction Yp = ArnYic/nB 
1 + exp(Am/TD) exp(t(TAr)/r„) 



with t{Ti\f) the value of time at T/v and r„ the neutron lifetime. 

On the other hand, the determination of all light nuclei produced during BBN, 
and a more accurate determination of ^He as well, can be only pursued by a 
simultaneous solution of a set of coupled kinetic equations which rule the evolu- 
tion of the various nuclei, supplemented by Einstein's equations, covariant con- 
servation of total energy momentum tensor, as well as conservation of baryon 
number and electric charge. This is typically obtained numeric ally, alt hough 



nice semi-analytical studies have been also recently performed (jMuk04l ) . 



To summarize the general BBN setting, we start with some definitions. We 
consider Nnuc species of nuclides, whose number densities, rij, are normalized 
with respect to the total number density of baryons ub, 

Xi = — 2 =p, ^H, ^He, ... . (20) 

ns 



The list of nuclides which are typically considered in BBN analysis is reported 
in Table 1. To quantify the most interesting abundances, those of ^H, ^He, 
'^He and ^Li, we also use in the following the short convenient definitions 

2H/H = X2h/X„ ^He/H = X3He/Xp, Yp = AX,ne, 'Li/H = Xtli/X^, (21) 



i.e. the ^H, ^He and ^Li number density normalized to hydrogen, and the '^He 
mass fraction, Yp. Notice that, although the above definition of Yp is widely 
used, it is only approximately related to the real helium mass fraction, since the 
'^He mass is not given by 4 times the atomic mass unit. The difference is quite 
small, of the order of 0.5% due to the effect of "^He binding energy. However, in 
view of the present precision of theoretical analysis on '^He yield, this difference 
cannot be neglected, and one should clearly state if one refers to the conven- 
tional quantity Yp (as we do here, too) or the "true" helium mass fraction. In 
the (photon) temperature range of interest for BBN, lOMeV >T> lOkeV, 
electrons and positrons are kept in thermodynamical equilibrium with pho- 
tons by fast electromagnetic interactions and are distributed according to a 
Fermi-Dirac function /e±, with chemical potential ±/ie, parameterized in the 
following by the ratio (pe = fJ'e/T. The chemi cal potential o f electrons is very 



small, due to the universe charge neutrality (jLyt59K ISen96l ). 



^ ^ = ^ 10-10 . (22) 
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Table 1 

Nuclides which are typically considered in BBN numerical studies. 
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To follow the neutrino-antineutrino fluid in detail, it is necessary to write 
down evolution equations for their distribution in phase space, rather than 
simply using their energy density. This is due to the fact, as we will illustrate 
in details in the following, that they are slightly reheated during the e"*" — e~ 
annihilation phase and develop non-thermal momentum-dependent features. 
We denote these distributions by /^edpM); /i^edpM) ^^'^ 



fuu = fur = fu^{\p\,t) , fu. = fur = /p.(IpM) 



(23) 



In the standard scenario of no extra relativistic degrees of freedom at the BBN 
epoch apart from photons and neutrinos, the neutrino chemical potential is 
bound to be a small fraction of the neutrino temperature. This bound ap- 
plies to all neutrin o flavors , whose distribution functions are homogenized via 
flavor oscillations (jDol02al : lAba02l : IWon02l ). In this Section we focus on non- 
degenerate neutrinos, namely f^^ = fp^ and /^^ = fp^, while we will consider 
in details the effect of neutrino-antineutrino asymmetry in Section 6.3. 



The set of differential e quations ruling primord i al nucleosynthesi s is the fol- 
lowing, see for example (IWagGTl : IWagGOl : IWag73l : lEspOObl : lEspOOal ) : 



a 



ISttG 



N 



riB 

p = -3Hip + P) 



(24) 



(25) 
(26) 
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^ij^kl 



Ub X! ^3 ^3 = '"'e- -ne+ = L [ 0e 



dt d\p\J 

where p and P denote the total energy density and pressure, respectively, 



(27) 
(28) 
(29) 



P = P^+ Pe + Pu + PB 
P = P^ + P, + P,. + P, 



B 



(30) 
(31) 



Eq. (24) is the definition of the Hubble parameter if, whereas Eq.s (25) and 
(26) state the total baryon number and entropy conservation per comoving 
volume, respectively. The set of N^uc Boltzmann equations (27) describes the 
density evolution of each nuclide specie, Eq. (28) states the universe charge 
neutrality in terms of the electron chemical potential, with L {■me/T,(f)e) the 
charge density in the lepton sector in unit of the electron charge, and finally 
Eq.s (29) are the Boltzmann equations for neutrino species, with I^^ [f^^, fyj\ 
standing for the collisional integral which contains all microscopic processes 
creating or destroying the specie z/^. 

Since electromagnetic and nuclear scatterings keep the non-relativistic baryons 
in kinetic equilibrium their energy density pB and pressure Pb are given by 



PB 

PB = TnB 



Ub 



(32) 
(33) 



with AMj and M„ the i-th nuclide mass excess and the atomic mass unit, 
respectively. 



The pressure and energy density of the electromagnetic plasma (e and 7) 
receives a contribution at first order in the fine structure coupling a when 
considering QED finite temperature corrections which change the electromag- 
netic plasma equation of state via the appearance of a thermal mass for both 
photons and e^, which in turn change the particle dis persion relation Ef.^.y{p). 
This has been studied e.g. in ( lHec94l : lLop99l : lMan02l ). These corrections also 
enter Eq. (26), the expression of the expansion rate H, as well as the ther- 
mal averaged n-p weak conversion rates, which depend upon electron-positron 
distribution function (see later). It has been shown that at the time of BEN, 
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all thes e effects slightly inffuence the ^He abundance, at the level of per mille 
flLop99h . 



In Eq. (27) k, I denote nuclear species, A''^ the number of nuclides of type 
i entering a given reaction (and analogously Nj, N^, Ni), while the F's denote 
symbolically the reaction rates. For example, in the case of decay of the species 
i, Ni = 1, Nj = and J2 ^i^ki is the inverse lifetime of the nucleus i. For two- 
body colhsions Ni = Nj = Nk = Ni = 1 and Tij^ki = {o'ij^kiv), the thermal 
averaged cross-section for the reaction i + j ^ k + I times the i — j relative 
velocity. In Eq. (28), Zi is the charge number of the i—th nuclide, and the 
function L{^,uj) is defined as 

oo 

iK,-)-^/^CCV^(^j=^-^j^) . (34) 



Eq.s (24)-(29) constitute a set of coupled differential equations which have 
been implemented in nume rical codes sinc e the pioneering works of Wagoner, 
Fowler and Hoyle (jWagGTl ) and Kawano (iKawSSi : lKaw92l ) . Before discussing 
some details of this implementation and its present status, we describe two 
crucial phenomena which take place before the BBN, the freezing out of neu- 
trino distribution functions and of the neutron-proton density ratio, when the 
weak charged current n-p processes become too slow to ensure chemical equi- 
librium between the two nucleon species. At the time these two phenomena 
occur the synthesis of deuterium, and thus of the whole nuclear chain, is still 
strongly forbidden by photo-dissociation processes, and all baryon density is 
in the form of free neutrons and protons. This means that in principle they 
could be treated independently of the whole set of nuclear processes which 
leads to the proper BBN phase. In particular, once obtained the shape of neu- 
trino distribution functions, they can be used as given inputs when solving 
the set of equations (24)- (25), thus significantly simplifying the problem. 



3.2 The role of neutrinos in BBN and neutrino decoupling 



At early epochs neutrinos were kept in thermal contact with the electro- 
magnetic primordial plasma by rapid weak interactions with electrons and 
positrons, controlled by a rate F^ ~ (cr^w)ne± ~ GjpT^ x T^. When the tem- 
perature dropped below a few MeV, these weak processes became too slow 
compared to the Hubble expansion rate (F.^, < ~ ^/GnT'^) and the process 
of neutrino decoupling took place. An accurate estimate for the neutrino de- 
coupling tempera ture is 2. 3 MeV for the electron neutrino and slightly larger 
for z/^,^, 3.5 MeV flEnq92bl ). since the latter only interact with the electromag- 
netic plasma via neutral current processes. 
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Shortly after, the e"^ pairs began to annihilate almost entirely into photons, 
thus producing a difference between the temperatures of the relic photons 
and neutrinos. The MeV to 0.1 MeV range is crucial for BBN physics, and in 
particular for the neutron-proton fraction, so it comes without surprise that 
BBN is sensitive to the properties of neutrinos: it derives from the basic facts 
that the neutrino decoupling, the deuterium binding energy, the n — p mass 
difference and the electron mass all fall in the MeV range. In more detail, 
neutrinos enter BBN equations in three ways: 

I the momentum distributions of the z/g and entering the n-p inter- 
conversion weak rates; 
II their overall energy density content p^, which determines the Hubble 
expansion rate; 

III their overall pressure P^, which enters the energy- momentum conserva- 
tion law. Assuming the equation of state of relativistic species = Pu/<i, 
this effect is not independent from the previous one, and we shall not dis- 
cuss it further. 



The effect (I) is clearly model-dependent, in the sense that it can be quantified 
and parameterized only on the basis of a physically-motivated hypothesis. One 
popular example where the effect of neutrinos is dominated by this distortion 
is the case of a i^e ~ J^e asymmetry, see Section 6.3. Concerning the effect (II), 
i t is customar y to parameterize it via an effective number of neutrinos N^s 
(IShv69l : ISteTTI ) , defined by the relation 



off 



-pl^fi 



(35) 



where p,yfl is the energy density in neutrinos in the limit of instantaneous 
decoupling from the e"'",e~,7 plasma and no radiative or plasma corrections 
(see Eq. (40) below). Nowadays, precision electroweak measurements at the 
Z°-resonance pin down the numbe r of lig ht active neutrino species with high 
accuracy, = 2.9840 ± 0.0082 flAleOSh . consistent within ~ 2cr with the 
known three families of the SMPP. However, well before these measurements 
were available, BBN was already invoked to favor 3 light, thermalized (and 
thus probably active with respect to weak interactions) neutrino fami l ies, wi th 
a range not wider than 2-4 (see Ref.s in the reviews flSarQGl : ISchQsl : lOliOOh l 
While this connection between collider physics and cosmology has been histor- 
ically important, we would like here to clarify some common misunderstanding 
and ambiguities on the BBN bound on A^^efr. If the bound is derived by chang- 
ing only Ncs, but without including any effect of the type (I), it is incorrect 
to refer to it as a bound on neutrino properties. It is rather a statement on 
the rate of expansion of the universe at the BBN time, which may be indica- 
tive e.g. of the presence of additional (semi)relativistic species in the plasma, 
or of exotic thermal histories. Any reasonable BBN bound on the number 
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of neutrino species Ny requires a modification of the type (I) as well. Since 
no universal parametrization exists, this is often neglected. Yet, this should 
not diminish the importance of the point. For example, if we accompany the 
rescaling of Eq. (35) by an identical rescaling of the and distribution 
function (i.e. adopting a gray-body parametrization, keeping the temperature 
of the spectrum unchanged but altering the normalization), the range of Yp 
that for standard parameters corresponds to 2.8 < A'cfr < 3.6 (see Section 5) 
translates into a more severe constraint, 2.85 < < 3.12. It is worth noting 
that no role is played in BBN by possible heavy fourth generation neutrinos, 
whose mass must be larger than ~45 GeV to avoid the Z-width bound. Even 
if such a neutrino exists, it would be natural to expect it to be unstable and 
rather short-lived, for example due to mixing with the light neutrinos (neu- 
trino oscillations prove that family lepton-numbers are violated). Even if some 
protective symmetry prevents the decay (so that it accounts for a subleading 
fraction of the cold dark matter), no effect on the BBN is present, and we do 
not discuss it further. 



A major consequence of the settlement of the neutrino oscillation issue is that 
a very refined calculation of the neutrino decoupling is possible. The standa rd 
picture in the instantaneous decoupling limit is very simple (see e.g. f lKolQOl )): 
coupled neutrinos had a momentum spectrum with an equilibrium Fermi-Dirac 
(FD) fornlE with temperature T, 



/eq(|p|,T) 



exp 




-1 



(36) 



which is preserved after decoupling. Shortly after neutrino decoupling the 
photon temperature drops below the electron mass, mg, and annihilations 
heat the photons. If one assumes that this entropy transfer did not affect 
the neutrinos because they were already completely decoupled (instantaneous 
decoupling limit), using conservation of entropy per comoving volume it is 
easy to calculate the difference between the temperatures of relic photons and 
neutrinos and thus the eventual neutrino energy density p^^ introduced before. 

At high temperatures, T > 2 — 3 MeV, one can write the conservation of 
entropy per comoving volume in the form (we are considering a temperature 



It was noted in (|Cuc96l ^ (and more recently the topic reanalyzed in (|Do105I )^ that 



the ^He abundance in the early Universe is sensitive to the difference between FD 
or an exotic Bose-Einstein (BE) distribution of the (quasi)thermal neutrino bath, a 
conclusion which is expected according to the spin-statistics theorem but otherwise 
difficult to prove in laboratory experiments. They found indeed a good sensitivity 
to the quantum nature of the statistics, Fp(BE)-yp(FD)~ -(3 A%)Yp{FD) (or 
equivalent to A'cff — 2.4), with the present range oiYp thus disfavoring a BE spectral 
shape. 
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range well below muon annihilation epoch) 

{se±,y + Su) = const . (37) 

After i/-decouphng, one has instead two separate conservation conditions for 
neutrinos and for the electromagnetic plasma. Nevertheless, until photons are 
reheated by annihilation, both photon and neutrino temperatures redshift 
by the same amount and keep equal. If one specifies the entropy conservation 
laws at the two different epochs, ain well before annihilation, and aend well 
after this phase, one obtains 



"^e^,7('-^m)'2jn ^'^{.^end)Ojend i ('^^) 

where in the second equation one takes into account that both photon and 
degrees of freedom contribute at aj„, while only photons are present (and 
reheated) at aend- The ratio of these two equations using the expression of en- 
tropy density (16), gives the well-known asymptotic ratio of neutrino/photon 
temperatures after annihilation phase, 

' ^ ^ 1.401 , (39) 



T ^2 + 4x7/8/ Vll 



and the instantaneous decoupling expression of the neutrino energy density in 
terms of p^, 

7/4 

P.,o = 3-(^-j ■ ^40) 



However, the processes of neutrino decoupling and annihilations are suffi- 
ciently close in time so that some relic interactions between and neutrinos 
exist. These relic processes are more efficient for larger neutrino energies, lead- 
ing to non-thermal distortions in the neutrino spectra (larger for Vf. than for 
z/^^r, since also feel charged-current interactions) and a slightly smaller in- 
crease of the comoving photon temperature. Even in absence of mixing, a 
proper calculation of the process of non-instantaneous neutrino decoupling 
requires the solution of the momentum-dependent Boltzmann equations for 
the neutrino spectra (29), a set of integro-differential kinetic equations that 
are quite challenging to attack numerically. In the last two decades a series 
of works has been devoted to solving this system in an increasingly general 
and precise way, ultimately also including finite temperature QED corrections 
to the elect romagnetic plasma (a full list of the pre-2002 works is reported 
m flDoin2bh ). To give a feeling of the overall effect, it suffices to say that 
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the combination of these corrections leads to an effective neutrino number 
iVeff ~ 3.046 , while asymptotically one f inds T„/ T = 1.398, slightly smaller 
than the instantaneous decoupling value (IManOSl ) . 



The existence of neutrino oscillations imposes modifications to these correc- 
tions. The effect of fi, distortions on the Yp yield in the simplified case of 
two-neutrino mixin g, aver a ged mo mentum, and Maxwell-Boltzmann statistics 
was esti mated in f KhlSlI : HanOl ). All three approximations were relaxed in 
f lManOSl ) and a full calculation was performed with a density matrix formal- 
ism. The neu trino ensembl e is described by the momentum-dependent density 
matrices Qp f lDolSlI : lRaf92l : ISig 9i: lMcK94l l The form of the neutrino density 
matrix for a mode with momentum p is 



Qp{t) 



^ Qee Qefi 
Qfie Ql^ifJ. 
\ Qre Qt/j, 



(41) 



The diagonal elements correspond to the usual number density of the different 
flavors, while the off-diagonal terms are non-zero in the presence of neutrino 
mixing. There exists a corresponding set of equations for the antineutrino den- 
sity matrix Qp. In absence of a neutrino asymmetry (or of additional couplings 
flipping z/'s into u and vice versa) antineutrinos follow the same evolution as 
neutrinos, and a single matrix suffices to describe the system. The equations 
of motion for the density matrices are 

iQp = [Ho + Hi, Qp] + C[gp\ , (42) 



where the first commutator term includes the free Hamiltonian Hq and the ef- 
fective potential of neutrinos in medium Hi, while the last term is a coUisional 
term of order G'jp describing the breaking of coherence induced by neutrino 
scattering and annihilation as well as neutrino production by collisions in the 
primeval plasma. In a FLRW universe, Qp = {dt — Hp dp) Qp. The Hamiltonian 
can be written explicitly in the flavor basis as 



= diag(m^, m 



I). 



(43) 



where U is the mixing matrix which, assuming vanishing CP-violating phases, 
writes 



U 



C12C13 ■5l2Cl3 Sl3 

~Si2C23 — C12S23S13 C12C23 — S12S23S13 S23C13 
S12S23 - C12C23S13 -C12S23 - S12C23S13 C23C13 



(44) 



17 



Here Cij = cos% and Sij = sin% for ij = 12, 23, or 13. Apart for CP- violating 
phases, there are five oscillation parameters: Amgi = ml — ml, Am^g = mf — 
"^2; ^12) ^23 and ^13. Best-fit values and uncerta inties for these parameters (and 



upper bound for ^13) can be found in (lYaoOGl ) . 



Neglecting non-diagonal components of the effective potential, the matter term 
writes 

H, = dmg{V,,V„Vr) , (45) 



where, assuming T ^ m^. 



V„ 



Ml 



8V2GfP 

72 Pi^iJ.+i'i^ ' 



3MJ 
8V2G 



3m: 



fP 

2 Put+Ut 



(46) 
(47) 
(48) 



Finally, the collisions of neutrinos with or among themselves are described 
by the term C[-], which is proportional to G^. For the off-di agonal t erms it is 
sufficient to adopt simple damping coefficient as reported in (jDol02al ). Instead, 
for the diagonal ones, in order to properly calculate the neutrino heating pro- 
cess one must consider the exact collision integral 1^^ , that includes all relevant 
two-body weak reactio ns of the type z/q,(1) + 2 ^ 3 + 4 involving neutrinos 
and e^, see e.g. flDolQTl ). The kinetic equations for the neutrino density matrix 
are supplemented by the covariant conservation equation for the total energy- 
momentum tensor. Given the order of the effects considered, one should also 
include the finite temperature QED corrections to the electromagnetic plasma 
(lMan02h . 



We show in Figure 1 the asymptotic values of the flavor neutrino distribution, 
both without oscillations and with non-zero mixing. The dependence of the 
non-thermal distortions in momentum is well visible, which reflects the fact 
that more energetic neutrinos were interacting with for a longer period. 
Moreover, the effect of neutrino oscillations is evident, reducing the difference 
between the flav or neutr ino distortions. Fitting formulae for these distributions 
are available in (IManOSl ). In Table 2 we report the effect of non-instantaneous 
neutrino decoupling on the radiation energy density, Nf,ff, and on the ^He mass 
fraction. By taking also into account neutrino oscillations, on e finds a global 
change of AYp ~ 2.1x10"^ which agrees with the results in (IHanOll ) due to 
the inclusion of QED effects. Nevertheless the net effect due to oscillations 
is about a factor of 3 smaller than what was previously estimated, due to 
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Fig. 1. Frozen distortions of the flavor neutrino spectra as a function of the comoving 
momentum, for the best fit solar and atmospheric mixing parameters. R is the scale 
factor. In the case where we allow for 7^ consistently with present bounds (blue 
dotted lines), one ca n distin guish the distortions for z/^ and (middle and lower, 
respectively). From ( ManOSi ) . 



Table 2 



(ManO.'^l^. 


Case 






No mixing (no QED) 


3.035 


1.47x10" 


-4 


No mixing 


3.046 


1.71x10" 


-4 


Mixing, (9i3 = 


3.046 


2.07x10" 


-4 


Mixing, sin2(6'i3) = 0.047 


3.046 


2.12x10" 


-4 


Mixing, Bimaximal, ^13 = 


3.045 


2.13x10" 


-^4 



the failure of the momentum-averaged approximation to reproduce the true 
distortions. 



It is worth remarking that the precise computation of the effect of particular 
neutrino decoupling scenario on primordial yields can only be performed nu- 
merically. The neutrino distribution functions, once obtained by the solution 
of Eq. (42), have to be substituted in Eq.s (24)- (28) which will predict the 
primordial abundance. The process is particularly involved, it is in fact im- 
portant to follow the single neutrino distribution as a function of time because 
of the particular role played by the different neutrino flavors. In particular, the 
A'^eff reported in Table 2 is the contribution of neutrinos to the whole radiation 
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Table 3 

Compar ison of the exact BBN results with a fixed-AA'^cfj approximation. From 
iManoi). 



Nuclide 


Exact (No z^-oscillations) 


Fixed AiVcff = 0.013 




1.71x10"^ 


1.76x10-"^ 


A(2h/H) 


-0.0068x10-5 


+0.0044x10-5 


A(3He/H) 


-0.0011x10-5 


+0.0007x10-5 


AfLi/H) 


+0.0214x10-^° 


-0.0058x10-^° 



energy budget, but only at the very end of neutrino decoupling. Hence, not 
all the AiVeff there reported will be really contributing to BBN processes. In 
order to clarify this subtle point, we report in Table 3 the effect on all light 
nuclides, of the non-instantaneous neutrino decoupling in the simple scenario 
of no neutrino oscillation, and compare this column with the simple prescrip- 
tion of adding a fix AA^^efr = 0.013 contribution to radiation. Even though Yp is 
reproduced (by construction), this is not the case for the other nuclear yields. 



3. 3 The neutron-proton chemical equilibrium and the role of weak rates 



Neutrons and protons are kept in chemical equilibrium by charged current 
weak interactions, 



(a) Ue -\- n —>■ e~ -\- p , (d) -\- p —>■ -\- n , 

(b) e~ + p ^ + n , (e) n e~ +77^ + p , 

(c) e'^ + n ^ Ve + p , (/) e" + 77e + p ^ n , (49) 

which enforce their number density ratio to follow the equilibrium value, n/p 
=exp(— Am/T). Shortly before the onset of BBN, processes (a) — (/) become 
too slow, chemical equilibrium is lost and the ratio n/p freezes out for tem- 
peratures lower than the decoupling temperature To ~ 1 MeV. Residual free 
neutrons are partially depleted by decay until deuterium starts forming at T/v 
and neutrons get bound in nuclei, first in deuterium and eventually in ^He. 

The leading role of (a) — (/) in fixing the neutron fraction at the BBN, and 
thus Yp, simply means that to get an accurate theoretical prediction for ^He 
abundance requires a careful treatment of the weak rates. Large improvements 
on this issue have been obtained in the last decade, which we summarize in 



the following. Extensive analysis can be found in e.g. (jLop99l : lEsp99l ) 



At the lowest order, the calculation is rather straightforward, and is obtained 
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10 



1 



T/m, 



Fig. 2. The total B orn rat es, 
(dashed hne). Prom ( Esp99l ) 



lob^ for n ^ p (soHd hne) and p ^ n transitions 



by using V-A theory and in thejimit of infinite nucleon mass (we will refer to 
this as the Born limit), see (jWei72l ). The latter approximation is justified in 
view of the typical energy scale of interest, of order T ~ MeV, much smaller 
than the nucleon mass Mpf. For example, the neutron decay rate takes the 
form (neglecting the very small neutrino masses) 



UJB{n^e + l^e +P) = ^ (Cy"^ + J d\Pe\\Pe\^\Pu\^ 

XeGp.l) [l-/p,(|p.|)][l-/e(|Pe|)] , (50) 

where Cy and Ca are the nucleon vector and axial coupling, and \p,y\ = 



Am — y I Pep + rrig. The rates for all other processes {a) — {d), (/) can be simply 
obtained from (50) by changing the statistical factors and the expression for 
neutrino energy in terms of the electron energy. An average can be performed 
at this level of approximation over equilibrium Fermi-Dirac distribution for 
leptons, i.e. neglecting the effects of distortion in neutrino-antineutrino distri- 
bution functions, but taking into account the time evolution of the neutrino to 
photon temperature ratio T^/T. In Figure 2 we report the Born rates, ub, for 
n-p processes. The accuracy of Born approximation results in being, at best, 
of the order of 7%. This can be estimated by comparing the prediction of Eq. 
(50) for the neutron lifetir ne at ve ry low temperatures, with the experimental 



value r^^ = (885.7±0.8) s (jYaoOGl) 



21 



The Born calculation can be improved by considering four classes of effects: 



Electromagnetic radiative correcti ons. Th ese are typically split into outer and 
inner terms (for a review see e.g. (jWil82l ) ) . The first ones involve the nucleon 
as a whole and consist of a multiplicative factor to the squared modulus of 
transition amplitude of the form 



a 



l + —g{E,,E,) 



(51) 



where analytic expression for g{Ee,E,^), can be found in Ref. (jSirGTl ). On the 
other hand, the inner corrections a re deepl y related to the nucleon structure. 
Th ey have been e stimated in Ref. (IMarSGl ). and applied in the BBN context 
m flLop99l : [Ei^ l Furthermore, when electron and proton are both either in 
the i nitial or final state, o ne should also add the effect of Coulomb correc- 
tion flPicSi ICam82l : iBaiQOh . due to rescattering of the electron in the proton 
electromagnetic field and leading to the Fermi factor 



1 + an 



E, 

|Pe 



(52) 



Finite nucleon mass corrections. For finite nucleon mass M^r and at order 
l/M^r, the weak hadronic current receives a contribution from the weak mag- 
netic moment coupling. 



■ Gf /2 

'y2M 



N 



■Up{p)a^u (p- qYun{q) 



(53) 



where, from conservation of vector current (CVC), /2 = VudifJ'p — /in)/2 = 
l.SlVud- Both scalar and pseudoscalar contributions can be shown to be much 
smaller and negligible for the accuracy we are interested in. At the same order 
I/Mtv the allowed phase space for the relevant scattering and decay processes 
gets changed, due to nucleon recoil. Finally, one has to consider the effect of 
the initial nucleon thermal distribution in the comoving frame. All these effects 
are proportional to me/Mj^ or T /Mj^., and in the temperature range relevant 
for BB N, can be as l arge as radiative corrections. This has be en first pointe d 
out in flSec93l : lLop97h and then also numerically evaluated in (jLop99l : lEsp99l ) . 



Thermal- Radiative corrections. The n-p rates get slight corrections from the 
presence of the surrounding electromagnetic plasma. To compute these correc- 
tions on e may use Real Time formalism for Finite Temperature Field Theory 
(jLeB96l ) to evaluate the finite temperature contribution of the graphs of Figure 
3, for the n — >■ p processes. Inverse processes p — >■ n are obtained by inverting 
the momentum flow in the hadronic line. The first order in a is given by in- 
terference of one-loop amplitudes of Figure 3 b) and c) with the Born result 
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n 



P n 



P 



a) 



b) 





c) 



Fig. 3. The tree level Born (a), the one-loop (b),(c) 
sion/absorption diagrams (d) for n — > p processes. 



and the photon emis- 



(Figure 3a)). Photon emission and absorption processes (Figure 3d)), which 
also give an order a correction, should be included to cancel infrared diver- 
gences. Notice that photon emission (absorption) amplitudes by the proton 
line are suppressed as and can be neglected. 

All field propagators get additional on-shell contributions proportional to the 
number density of that particular specie in the surrounding medium. For 7 
and (neglecting the small electron chemical potential) we have 



iAf (fc) 

iSe{Pe) 



I 

i 



rrie 



271 5ie) f,{k) 
- 2iT S{p. 



/7V 
2 



(54) 
(55) 



with /-y the photon distribution function. The whole set of thermal/radiative 



Don84: 


Don85: 


,Toh8fil: Kei89a: 


Bai9nl: 


Kei89b: 


LeB9nl: Alt89: 


Kob85 


Saw96: 


ChaQTl: 


Esp98l'). Though thev aeree on the order of magnitude 



Don83 



Kob8fi 



is quite small — there is nevertheless no consensus on the detailed value, due 
for example to different ways of treating the wave function r enorma lization 



at finite temperature. Finally, it was correctly pointed out in ( iBroOll ) that a 
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(small) contribution to neutron-proton chemical equilibrium is also provided 
by photon-proton interactions (and inverse processes), 



7+P 



+ Ue + n 



Ue + n 'y + p 



(56) 



The corresponding thermal averaged rates can be found in e.g. flSer04bl l 



Non-instantaneous neutrino decoupling effects. Distortion of neutrino distri- 
bution functions changes the weak rates (a) — (/) which are enhanced by the 
larger mean energy of electron neutrinos. On the other hand, there is an op- 
posite effect due to the change in electron-positron temperature. Finally, since 
the photon temperature is reduced with respect to the instantaneous decou- 
pling value (11/4)^/^, the onset of BBN, via synthesis, takes place earlier 
in time. This means that fewer neutrons decay from the time of freezing out 
of weak interactions, and this in turn corresponds to a larger ^He yield. 

The effect of all correcti ons to the we a k rates discussed so far has been 
considered in details in fjLop99l : lEsp99l : ISer04bl ). The leading contribution 



is given by electromagnetic radiative corrections, which decrease monotoni- 
cally with increasing temperature for both p — > n and n ^ p processes, and 
by finite nucleon mass effects. Their sum changes the Born estimate for a 
few percent correction at the freeze out temperature T ~ MeV. Comparing 
the theoretical prediction for the neutron lifetime at zero temperature us- 
ing Gf = (1.16637±0.00001)_10liGe1/-2, Cy = 0.9725±0.0013 and the ratio 
Ca/Cv = -1.2720±0.0018, flYaoOGh . one finds rf = 886.5 s, quite an accurate 



result when compared with the experimental value (agreement is at the 0.1% 
level). It is worth commenti ng here on the fact that a recent measurement of 
the neutron lifetime exists (ISer05al ). which gives r„ = 878.5 ± 0.7stat ± O.Sgyst 
and results in a 5. 60" discrepancy from the previous most precise re sult, wh ich 
in turn is consistent with the other six determinations used in ( lYao06l ) to 
determin e the best fit. According to the Particle Data Group, [The result of 
lSer05a) l is so far from other results that it makes no sense to include it in the 
average. It is up to workers in this field to resolve this issue. Until this major 
disagreement is understood our present average of 885.7 ± 0.8 s must be sus- 
pect. Wh ile implications for BBN of this different lifetime have been explored 
( lMat05al ). in the following we shall assume that the best fit provided by the 
PDG is correct. Note that the PD G valu e also appears to be in better agree- 
ment with global electroweak fits (jSevOa ). Modulo this caveat, one might be 
confident that n-p weak rates are presently quite accurately computed, at per 
mille precision. Plasma corrections and finite temperature radiative effects are 
sub-leading, changing the rates at the level of (0.3 -f- 0.6)% only. Their effect 
is to slightly increase Yphj a. very small amount, AYp ~ 1-10~^, (IFieQSl ) . 



To conclude, we report a fit of the n-p rates w hich inc lude all effects described 
in this Section as a function of z = rrie/T, (ISer04bl ). accurate at the 0.01% 
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level, which the reader might find useful: 




1 

exp {-qnp/z) Y,ai z'^ 0.01 < T/MeV < 10 , (57) 



^ exp E/=i hi 0.1 < T/MeV < 10 

0.01 < T/MeV < 0.1 



with 



(5J 



ao = 1 ai = 0.15735 = 4.6172 

03 = -0.40520-102 a4 = 0.13875-103 as = -0.59898-102 
as = 0.66752-10^ 07 = -0.16705-10^ ag = 3.8071 (59) 
og = -0.39140 aio = 0.023590 an = -0.83696-10"^ 

= -0.42095-10"^ 013 = 0.17675-10"^ g„p = 0.33979 , 



60 = -0.62173 61 = 0.22211-10^ 62 = -0.72798-10^ 
63 = 0.11571-10^ 64 = -0.11763-10^ 65 = 0.45521-10^ 

(60) 

6g = -3.7973 67 = 0.41266 = -0.026210 

h = 0.87934-10-3 bio = -0.12016-10"^ Qpn = 2.8602 . 



3.4 Nuclear Reaction Network 



Nuclear processes during the BBN proceed in an environment very different 
with respect to the perhaps more familiar stellar plasmas, where stellar nucle- 
osynthesis takes place. The latter is a dense plasma where species are mostly 
in chemical equilibrium, the former is a hot and low density plasma with a 
significant population of free neutrons, which expands and cools down very 
rapidly, resulting in peculiar "out of equilibrium" nucleosynthetic yields. The 
low density of the plasma at the time of BBN is responsible for the suppres- 
sion of three-body reactions and an enhanced effect of the Coulomb-barrier, 
which as a matter of fact inhibits any reaction with interacting nuclei charges 
Z1Z2 ^6. The most efficient categories of reactions in BBN are therefore pro- 
ton, neutron and deuterium captures -(p, 7) {n,'j), (c?, 7)-, charge exchanges 
-(p, n)-, and proton and neutron stripping -{d,n), {d,p). 



25 



1 

-2 
10 

-4 
10 

-6 
10 

-8 
10 

-10 

10 

-12 

10 

-14 

10 

-16 
10 

-18 

10 



10 20 30 40 50 60 70 80 90100 

1/T (MeV'O 

Fig. 4. The evolution of some element abundances produced during BBN. The line 
M/H refers to the abundance of "metals" (see text). 

From a more technical point of view, accurate BBN predictions require a de- 
tailed knowledge of the nuclear rates entering the set of Equations (27): a 
physical understanding of the basic reaction chains helps however to imple- 
ment a code with a good compromise between reliability and computational 
time. In this Section we start with the phenomenology of nuclear processes 
related with the choice of such nuclear network, whereas a chronology of the 
codes devoted to solve the whole set of equations, and related numerical issues, 
is reported in Section 3.4.1. 

After the freeze-out of weak interactions, p(n, 7)^H is the only reaction able 
to synthesize sensible amounts of nuclei, since it involves the only two nuclear 
species with non-vanishing abundances, protons and neutrons. The very high 
entropy of the primordial plasma keeps this reaction in equilibrium down to 
energies much lower than its Q-value, which corresponds to the deuterium 
binding energy of Bp ~ 2.2 MeV. The deuterium bottleneck ends at a tem- 
perature of Td ~100 keV, given roughly by the condition exp(i?/5/T£))?7 ~ 1, 
which ensures that the high energy tail of the photon distribution has been 
sufficiently diluted by the expansion. At this point, almost instantaneously 
all the available neutrons are locked into deuterium nuclei. The effect of the 
bottleneck on the following nuclear processes is two-fold: the "delayed" pop- 
ulation of the deuterium specie results in BBN taking place at relatively low 
temperatures, with consequences on the efficiency of all the following reac- 
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Table 4 

The most relevant reactions for BBN. 



oy lllUOi 


IXtJclCLlOll 


oy liiuui 


l-< 00 r^Ti /~\n 


Rq 




i?8 


'^He(a,7)^Be 


Rl 


f \ 7 


Rg 


'^H(a,7)'Li 


R2 


''H(p,7)'^He 


RlO 


'Be{n,p)'Li 


R3 


''H(d, ?i)'^He 


Rii 


'Li(p, a)*He 


i?4 


2H(d,p)3H 


R12 


'^He(d, 7)''Li 


i?5 


^He(n,p)'^H 


Rl3 


^U{p, a)3He 


Rq 


3H(d,n)4He 


Rii 


^Be(n,a)^He 


R7 


3He(d,p)^He 


Rib 


^Be(d,p)2^He 



tions. Second, but not less important, is the intervening decay of the neutron, 
which alters the neutron to proton ratio at the time of effective deuterium 
production. Efficient production of deuterium marks the effective beginning 
of the nuclear phase of BBN: the nuclear path towards heavier elements, be- 
fore all ^H, ^He, ^He, is enabled by the interaction of deuterium nuclei with 
nucleons and other nuclei. The bulk of the remaining ^Li is produced via 
tritium and especially ^He radiative capture on ^He. The latter path leads to 
^Be, which eventually decays into ^Li by electron capture. The evolution of 
light nuclide abundances as a function of the temperature of the plasma is 
shown in Figure 4: they undergo the nuclear phase of departure from chemical 
equilibrium and, by the time the universe cools to few keV, they reach their 
final values. It is worth noting the dramatic change that all species undergo 
following the deuterium synthesis at 1/T>14 MeV~^. 



The most important nuclear processes for the BBN were identified in (jSmi93al ) 
(reactions Rq-Rh in Table 4). Reaction Rg provides actually only a sub- 
leading contribution to ''Li for today's preferred value of 77, being more rele- 
vant at lower 77 where ^Li production is direct rather than coming from ^Be 
synthesis. Also, if one is concerned with the traces of ^Li produced in BBN, 
the reactions R12 : ^He((i, 7)^Li, and -R13 : ^Li(p, a)^He are relevant. Due 
to their larger uncertainties, even reactions as -R14 : ^Be(?7,, a)^He and per- 
haps -R15 : ^Be((i, p)2 ^He are important in the ^Li error budget determi- 
nation. All these reactions are summarized in Table 4 and Figure 5. Both 
l eading and sub-leading reactions of some interest have been discussed e.g. in 
( Ser04br ) , which we mo stly follow here. Other compilations can be found in 
flAngggl : [Cvboi : lDes04[ ). 



These reactions are not only important to understand the nuclear physics of 
the BBN, but also to assess an error budget for its th eoretic a l predic tions. 
The first "modern" papers addressing these issues are (iKraQOl : ISmi93al ). For 
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Fig. 5. The most relevant reactions for BBN. 



example, in (jSmi93al ). the authors performed a systematic analysis of the 
nuclear network in BBN, as it had been implemented in the first publicl y 



released code (IWagGTl ). with the-at the time-new rates compiled in (ICau 
For the range of 77 allowed at the time, they identified the twelve leading 
reactions reported in Table 4. A detailed study of their rates and uncertainties 
led to a new code (jKaw92l ). implementing an updated nuclear network. 



Inferring the uncertainties on the light elements yields is conceptually a three- 
step process. First, one has to determine the uncertainties on the cross-sections 
ai{E) (as functions of energy) measured in the Lab or theoretically predicted. 
Since a first principle computation is virtually impossible for almost all pro- 
cesses, one has to rely on some form of theoretically motivated fitting for- 
mulae and interpolate between and extrapolate beyond the data. Indeed, 
the reactions in the early universe happen in a thermal plasma, so the rel- 
evant nuclear input are the temperature-dependent rates given by the convo- 
lution of the cross-sections with the Maxwell-Boltzmann distribution of nu- 
clides, i.e. {av) oc T~^/'^ dE a{E) E exTp{—E /T). These can be approxi- 
mated with analytical formulae in some simple hypotheses for the functional 
form of <j{E). Although no modern compilation relies on such approxima- 
tions, they are often used to suggest a fitting formula for the numerical inte- 
gration results. Anyway, it is the uncertainties on these rates that ultimately 
propagate onto the final errors on the nuclides. The f orr nalism /m achinery to 
treat this prob lem can be found in classical papers (as flSmi93ah ). textbooks 



f Cla83 



RolSSf) and ha s also been rep orted in several compilations of the last 



decade (lAngQQl : [Cvbo3 : lDes04l : [s"er04bl ) . so we do not repeat it here. We want to 
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remark, however, that: a) In the last decade, in particular following (jAng99l ). 
a world-wide effort in obtaining new measurements at low energy for several 
important reactions involving light nuclei has been undergone; b) in several 
cases, especially when new measurements have become available, the different 
experimental data sets do not seem to agree within the quoted uncertainties. 
In this situation, it is a tricky business to assess uncertainties in a statistically 
meaningful way, unless one has reason to believe that some of the datasets are 
affected by unaccounted systematics and decides for example to dis miss the 
older measurements. Here we follow the prescription illus trated i n flSer04bl ) 
and motivated on the basis of the arguments presented in (lDAg94f). While in 
agreement with the error estimates given in other compilations (jCyb04l : lDes04l ) 
when statistical errors dominate, we warn the reader that our procedure tends 
to produce smaller errors than other compilations when discrepancies among 
datasets exist. 



An additional technical aspect arises in the way uncertainties are accounted 
for in the BBN calculations. A lot of attention has been paid in the last two 
decades to this problem. One may a dopt Monte Carlo simu lations directly, 
with various degrees of sophis tication flKra90l : ISmi93al : iNolOOh , or use an error 
matrix approach as in (jFio98l ) and its generalization in (JCuoOJ). Each one has 
its own advantages and dis advantages, but the agreement between the two is 
typically very good (jFio98l ). In Table 5 we report the light ii uclide abundances 
for WMAP 5-year result fi^^^ = 0.02273 ±0.00062 flPunOSh (second column), 
showing in the third and fourth columns, respectively, the uncertainties due 
to l^^/i^ i^ojt) nuclear rates errors (cm)!^. The last two columns report 
the rates mostly contributing to the nuclear uncertainties, and the ir relative 
contributions in percent in a quadrature sum, according to flSer04bl l 



To illustrate the main dependence of the yields o n the nuclear rates, we report 
here the scaling relations introduce d in (|Cvb04l ) and normalized with respect 
to the predictions of PArthENoPE ( ParOSi ) around the fiducial value of Ub = 
0.02273 (WMAP 5-years). The scalings are: 



0.41 



^ = 2.53 X 10-^ R^^-^Rf-^Rf-^^R^'^-^ (^^) (^] (61) 



^He 
IT 



1.02 X 10-^i2f°-^^/2° =^S/?4-°-25/23 -20i?-o.i7^o.o8 



0.02273 



-0.59 



0.15 



Jnfi, 



(62) 



^ A recent analysis finds a larger valuejbr ^Li/H= (5.24l;Q'g7)10~^'^, due to new 
determinations of the rate Rg, see ( CvbOSi ) and references therein. This further 
exacerbates the ^Li problem, see Sections 4.4-4.5. 
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Table 5 

The light nudide abundances for WMAP 5-year result Ob/i^ = 0.02273 ± 0.00062 
(second column). The uncertainties due to Qs/i^ {(^uji,) and nuclear rates errors (an) 
are also shown in the third and fourth columns, respectively. The last two columns 
report the rates mostly contributing to the nuclear uncertai nties, an d their relative 
contributions in percent in a quadrature sum, according to ( Ser04bl ) . 



nuclide i 


central value 




(Tii 


rate 






0.2480 


+0.0002 
-0.0003 


±0.0002 


Rq 


98.5 










R2 


49 


^H/H xlO^ 


2.53 


±0.11 


±0.04 


R'^ 
R^ 


37 
14 


3He/H xlO^ 


1.02 


+0.01 
-0.02 


±0.03 


Rr 

R2 


80.7 
16.8 










i?14 


40.9 


^Li/H xlO^o 


4.7 


±0.3 


±0.4 


Rs 

Rl5 

Ri 


25.1 
16.2 
8.6 


^Li/H xlO^"^ 


1.1 


±0.1 


+ 1.7 
-1.1 


Rl3 


~ 100 



Fp = 0.2480 



0.02273 



0.39 



0.72 



7t ■, 

— = 4.7x lO-^'^i^i'^^i^pi^f 
H 



10 dI.34 dO.96 d-0.76^-0. 71^0.71 ^0.59^-0.27 



0.02273 



2.12 



0.44 



JnX), 



(63) 



(64) 



It is clear, for example, that prediction is dominated by the neutron mean 
lifetime (see Section 3.3). For all the other reactions the relative weight of the 
different cross-sections is consistent with what is reported in Table 5, modulo 
a caveat: some important reactions may have a minor impact on the uncer- 
tainty due to their better determination or vice versa. This is the case of -Ri, 
whose role does not reflect in the error budget being theoretically well under 
control (at the % level), or conversely the case of Rm,Ri5, which contribute 
ap preciably to the ^Li error budget in Table 5 due to the assumption done 
in ( Ser04bl ) that they are only known at the order of magnitude. In general, 
while improvements in the nuclear reaction rates would still sharpen the BBN 
predictions (especially for ^Li), it is fair to conclude that these uncertainties 
are at the moment negligible compared to the observational ones (see Section 
4.4). Thus an experimental campaign does not seem mandatory for BBN pur- 
poses alone. Yet, a careful assessment of the systematic uncertainties would be 
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useful in reducing the present discrepancies in data regression protocols. A re- 
alistic account especially of scale uncertainties in the already existing datasets, 
perhaps excluding unreliable measurements and correcting underestimated er- 
ror assignments would be certainly a more useful input from the experimental 
community. 



It is worth commenting here on an additional robust prediction of BBN: el- 
ements heavier than ^Li are virtually absent from the chemical composition 
of the early univers e. This was realized from the first pioneering studies on 
BBN flAlp48l : iFerSoh . and mostly explained in terms of the very high entropy 
and on the low density of the primordial plasma. They result in extremely low 
abundances A>2 elements at the weak reaction freeze-out, and an inhibition 
of three-body reaction^, respectively. To assess more quantitatively the ro- 
bustness of this result, the present authors have recently studied t he nuc lear 
processes involved in the production of A > 8 elements in BBN ( llocOTl ). It 
was confirmed that their synthesis can only start at T<60 keV (see Figure 
4), close to the freeze-out of nuclear reactions, when they have a very small 
efficiency. Heavier elements in BBN are in fact produced by an a capture over 
^Li and subsequent build-up of ^^C; their final abundance, mainly produced 
by means of an a capture over ^Li and subsequent proton capture, is to be well 
below 10~^° that of hydrogen (see Figure 4 for details). While processes till 
now neglected were identified and included, their role does not alter this basic 
conclusion. An important implication of these results follows for the first stars 
in the universe, known as Population III stars, whose formation process and 
initial mass function is thought to be critically dependent on the carbon and 
oxygen content of the environment^. While the amount of "heavy" element 
predicted in SBBN is safely below the critical level to alter the standard (but 
yet to be tested) scenario, in exotic models this might not be true. More in 
general, even when roughly reproducing ^H, ^He, ^He and ^Li yields of the 
SBBN, exotic scenarios may differ for the predictions of ^Li, ^Be or elements 
with A > 12 ('metals' in the astrophysical jargon). Examples of this will be 
provided in Section 7.1 and 9.2.3. 



3.4-1 Numerical solution of the BBN set of equations 

Despite of the fact that some simple estimate of BBN predictions can be 
made on the back of the envelope, detailed theoretical estimates that can 
be compared with experimental data require a careful numerical solution of 
BBN equations (24)- (29). This is particularly relevant when one considers 
exotic scenarios with extra parameters or where a different BBN dynam- 



^ The leading reaction for the synthesis of ^^C in stars is the 3a — > ^^C. 
^ For a review we addres s the reader to the proceedings of the last plenary confer- 
ence on 



the topic (IfS308I 1. 
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ics is considered, with the aim of constrai ning physics be yond the SMPP. 
Since the original Wagoner computer code (IWag69l : IWag73l ) much effort has 
been devoted to develop numerical tools and provide reliable numerical results 
flKaw88l : |K^^:lFio98[ iLoggSl : lEsp99l : iBurOol : lEspOOal : lOliOOl : ICvbOll : ICvb03bl : 
CocOi ICvb04l: ISer04bl : IPisOSl l In 1988, Kawano modified the BBN program 
of Wagoner (jKaw88l ) . a nd in 19 92 an updated and user friendly public ver- 
sion has been released flKaw92h . which has served reference tool for a 
wide num ber of research groups. Nuclear reaction rates were updated in 1993 
( 1Smi93al ) and radiative and Coulomb corrections to the weak rates, to which 
the ^He abundance is very sensitive, as well plasma effects, were included as 
a constant multiplicative factor to n-p weak rates. In 1993 Kernan discovered 
a relatively large time-step systematic error in the '^He abundance prediction 
of the public Kawano code, 6Yp = 0.0017, since then routinely added to its 
result. This, strictly speaking, is not fully adequate because different users use 
different step sizes and furthermore, numerical error is mach ine dep endent. 
In 1999 Lopez and Turner wrote a new nucleosynthesis code ( Lop99l ). which 
used the same nuclear rates but incorporated more accurately the radiative, 
finite nucleon mass and finite temperatur e corre ctions to the weak rates. The 
code used by Olive, Steigman and Walker flOliOOf ) agreed, over the whole range 
1 < '7io ^ 10, better than the 0.1% level, with the predicted ^He abundance 
of the Lopez and Tu rner code. In 2000, the weak rate corrections were also 
calculated in (lEsp99l ) and included in a BBN code. This code was developed 
and continuously updated over almost a decade, giving particular care also to 
the treatment of neutrino decouplin g and the nuclear reaction chain which en- 
ters the light abundances evolution (lEspOOal : ISer04br) . It was eventually made 
public in 2008 with the name PArthENoPE f|Par08[ ). Details can be found in 
( lPis08l ). where particular emphasis is given to a comparison with (|Kaw92K 
To our knowledge, only the original Kawano code or improved versions of il0 
and PArthENoPE are publicly available. 



To examine the concordance of theoretical predictions we have considered 
some recent results obtained using numerical BBN codes, either publicly avail- 
able or whose outputs have been made public by the authors as fitting formulae 
versus the baryon density or finally, simply quoted in their papers for some 
reference value of Qsh"^- These results are shown in Table 6. The list is of 
course largely incomplete, but we think it is representative enough to give the 
idea of the status of BBN theoretical accuracy. In the last column we also 
quote the theore tical un certainty due to nuclear rate (experimental) errors, as 
estimated using (jPar08l ). We see that all results are in very good agreement, in 
particular Yp agrees at the ±0.1% level in the most recent calculations. There 
is a somehow larger spread of ^Li estimates which are however all compatible 
within the larger theoretical uncertainty (Table 6). 



See the hnk http://www-thphys.physics.ox.ac.uk/users/SubirSarkar/bbn.html 



for a version of the Kawano standard code, made Linux-friendly by S. Dodelson. 
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Table 6 

A comparison of light nuclei theoretical predictions in some recent analyses. Results 
are either produced by public numerical tools, or obtained by using fitting formulae 
made available by the authors o r finall y, simply quoted in the papers. The theoret- 
ical errors are estimated using ( ParOSi ) . and account for the effects of nuclear rate 
uncertainties. Results are shown for a baryon fraction 0^/1^=0.0224. 





('Kaw92) 


('Fio98l 


fBurOO) 


fCvb04^ 
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('Par08) 
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0.2479 


0.2483 


0.2485 
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2.60 


2.60 


2.55 


2.59 


2.58 ± 0.04 


^He/H xIQS 


1.04 
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1.03 ± 0.03 


■^Li/H xlO^o 


4.53 


4.42 


4.91 


4.26 


4.50 


4.57 ±0.4 
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Coefficients 


of the fit of Eq. (65) for the Yp abundance. 
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Table 8 


















Coefficients 


of the fit of Eq. (65) for the ^H/HxlO^ abundance. 
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14.892 


-1551.6 


70488. 


-1.5390x10^ 1.3164x10^^ 




1 




6.1888 


-916.16 


56639. 


-1.6046x10^ 1.7152x10^^ 




2 




-0.60319 


118.51 


-8556.3 


267580. 


-3.0624x10^ 




3 


4.5346x10-2 


-8.7506 


624.51 


-19402. 


221200. 





We report the fit of the main nuclide abundances as function of cj;, = Vtsh"^ 
and the number of effective degree of freedom, iVefj. The fit holds for 0.015 < 
uji, < 0.029 and < A^^cfr < 7. The fitting function for all nuclei is 



EE 



(65) 



and the coefficient are reported in Tables 7-10. The fit accuracy is better than 
0.13% for Yp, than 0.3% for ^H/H and ^Li/H, and than 0.6% for ^He/H. 
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Table 9 



Coefficients 


of the fit of Eq. (65) for the '^He/HxlO^ abundance. 
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Table 10 












Coefficients 


of the fit of Eq. (65) for the ''Li/HxlO^'^ abundance. 
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-4.2339 
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4 Observational Abundances 



The abundances of primordial elements are inferred from measurements per- 
formed in a large variety of astrophysical environments. What 'precision cos- 
mology era has meant in this field is an increased number and precision of 
spectroscopic data over the past two decades. Presently, the situation is still 
quite involved due to the presence of relevant systematic errors which are 
comparable to (where not dominant over) the statistical uncertainties. Un- 
fortunately, these errors are to a large extent irreducible and intrinsic to the 
astrophysical determinations themselves, which rely on highly evolved systems 
where reprocessing, astration, and contaminations from younger systems are 
possibly present and difficult to correct for. 

Since our nearby universe is far from reflecting its primordial conditions, the 
methods proposed and developed in the last forty years to infer primordial 
yields have focused on very old and hence little evolved astrophysical regions, 
as well as on the capability to correct for the effect of the galactic evolution 
on the pristine abundances. For example, due to its very weak binding, any 
2H nucleus contained into pre-stellar nebulae is burned out during their col- 
lapse. Hence, the post-BBN deuterium evolution is expected to be a monotonic 
function of time and astrophysical deuterium measurements can be assumed 
to represent lower bounds of its primordial abundance. 
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Unfortunately, such a simple scheme cannot be applied to the more tightly 
bound ^He nucleus. In this case, in fact, in stellar interior it can be either 
produced by ^H-burning or destroyed in the hotter regions. As a consequence, 
all the '^He nuclides surviving the stellar evolution phase contribute to the 
chemical composition of the InterStellar Medium (ISM) and thus stellar and 
galactic evolution models are necessary to track back the primordial ^He abun- 
dance from the post-BBN data, at least in the regions where stellar matter is 
present. For this reason the inferred primordial ^He abundance is intrinsically 
a model-dependent quantity. 

The situation of ''He is reversed with respect to that of deuterium, since in 
this case the hydrogen burning in the successive stellar population has in- 
creased the amount of ^He as well as "metals" (nuclei with Z > 4) such as 
C, N and O. Usually, ^He is measured in old and very little evolved systems 
versus their metallicity, and extrapolating linearly to zero-metallicity. While 
this is certainly a reasonable approach, still it could lead to some systematic 
uncertainty. 

Finally, ^Li, whose bulk is also believed to be produced in the primordial 
cauldron, is a very weakly bound nuclide which has an extremely involved 
post-BBN chemical evolution. In fact, it is easily destroyed in the interiors of 
stars but can survive in the cooler outer layers of stars with shallow convec- 
tive zones, where it can be measured by means of absorption spectra. However, 
the scenario is much more involved due to the observation of enhanced lithium 
abundance in some red-giants. This suggests that ^Li formed in the interior of 
some stars may be transported by convective modes to the cooler exteriors. In 
this case these stars behave as lithium producers and thus enrich the ISM with 
lithium. Furthermore, the cosmic rays scattering on ISM nuclei can contribute 
to the total amount of ''Li and since CNO elements are necessary for spalla- 
tion processes, this would imply a correlation between post-BBN ''Li and the 
metallicity, which may help in tracking back this non-primordial component. 



4-1 Deuterium 



It is commonly believed that there are no astrophys ical sou rces of deuterium 
since it is destroyed by stellar evolution processes flEps76h and non-thermal 
production channels have been constrained to be negligible (see e.g. f lPro03h l 
Thus, any astrophysical observation can provide a lower bound for the primor- 
dial abundance. Therefore, the local ISM in the Milky Way can provide an 
order of magnitude more determinatio ns of ^H /H than high redshift Quasar 
Absorption Systems (QAS) (see Ref. (ILinOGl ) for the most recent compila- 
tion), given its easier observation al acces sibility. By using the Far-Ultraviolet 
Spectroscopic Explorer (FUSE) (IFUS99I ). a large database of Galactic ^H/H 
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measurements has been compiled. However, despite FUSE and other satelhte 
observations, providing measurements of ^H/H in almost 50 lines of sight, the 
picture of Galactic deuterium abundances remains puzzling. In particular, in- 
side the Local Bubble (<100 pc from the Sun) the deuteriur n to hyd rogen ratio 
seems roughly constant at ^H/H|p = (1.56 ± 0.04) x 10"^ (jWooOJ). However, 



beyond this bound , an unexpected scatte r of a factor of ~ 2 in ^H/H values 



is observed (|Jen99l : IWoo04l : lLin06l : lEUOTl ) as well as correlations with heavy 
element abundances, which suggest that ISM deuteri um mig ht have suffered 
stellar processing, namely astration (see for example (jRom06l )). but also that 
it may reside in dust particles which evade gas -phase observations. This is 
supported by a measurement in the lower halo (jSavOTI ) which indicates that 
the Galactic ^H abundance has been reduced by a factor of only 1.12 ± 0.13 
since its formation. As an alternative explanation, it is worth reporting the 
possible existence of a strong late infall of pr e-galactic mat erial with primor- 
dial composition (high ^H) (see for example ( PraO?! : Pro08 )) which has some 
observational evidence via the study of kinematics of high latitude gas regions. 
Finally, it is important to mention the analysis of Infrared Space Observatory 
spectra of H?, H-^H , CH4, and CHa^H in Jupiter's atmosphere which led au- 



thors of Ref. flLelOll ) to infer the value of deuterium to hydrogen ratio for the 
protosolar cloud ^H/Hp^c = (2.1 ±0.4) x 10~^, which corresponds to our galaxy 
value when its age was two thirds of the present one. 



The astrophysical environments which seem most appropriate to obtain reli- 
able measurements of the primordial deuterium fraction are the hydrogen-rich 
clouds abs orbing t he light of background QSOs at high redshifts, as recognized 
already in (jAda76l ). Conventional models of galactic nucleos ynthesi s (chemical 
evolution) predict a small contamination of pristine ^H/H fjFie96l ). However, 
a successful implementation of this method requires: (i) neutral hydrogen col- 
umn density in the range 17 < log[N(H^)/cm < 21; (ii) low metalhcity 
[M/H] to reduce the chances of deuterium astration; (iii) low internal velocity 
dispersion of the ato ms of the clouds, allowing the isotope shift of only 81.6 
km /s to be resolved ( PetOSi ). For this reason, only a handful of determinations 
have been obtained since the advent of the > 8m class telescopes in the 1990s 
in damped Lyman-a systems and Lyman limit systems: 



Q 1009-1-2956, with the absorber placed a t z 
X 10-5 (log^H/H 



(3.98t°;i) 



2.504, yielding ^H/H 



-AAOtofs) flBur98l : lKir03l l 
ii) PKS1937-1009 (I), with the absorber placed at z = 3.572, yielding ^H/H 
(3.3 ± 0.3) X 10-5 (log^H/H = -4.49 ± 0.04) flBur98h . 



iii) 



IV 



= 2.536 with metallicity 
5 (log^H/H = -4.596 ± 



HS 0105-1-1619, with the absorber placed at z 
[Si/H] ^ 0.01, y ielding ^H/H = (2.54 ± 0.23) x 10 
0.040) flOMeOlh . 

Q2206-199, with the absorber placed at ^ = 2.0762 with metallicity [Si/H] 
-2.23 yielding ^H/H = (1.65 ± 0.35) x 10-^ (log^H/H 
(jPetOlh . 



4.78tH^) 
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VI 



Vll 



Vlll 



IX 



XI 



Xll 



Xlll 



XIV 



Q0347-3819, with the absorber placed at z = 3.024855 with metaUicity 
[Si/H]= -0.95 ± 0.02 , yielding ^H/H = (3.75 ± 0.25) x I Q-^ (log ^H/H = 
-4.43 ± 0.03) fllevOlh . Note that the previous analysis in flPOdOlh yielded 
^H/ H = (2.2 4ib0.67) x 10~^ due to an incorrect velocity distribution function 
(see (ILevOll )). Similar considerations apply for the value discussed in vii) 
below) . 

Q1243+3047, with the absorber placed at z = 2.525659 with metaUicity 
[0/H]= -2.79 ±0.05, yielding ^H/H = (2.42l[J:i) x lO'^ (log^H/H = 



-4.617l°:°i) flKir03h . 



PKS1937-1009 (II), with the absorber placed at z = 3.256 with metal- 
hcity [Si/H] = -2.0 ± 0.5, yielding ^H/H = [l.Gtoio) x lO'^ (log^H/H = 
— 4.80lo;o9) dCriOJ). Since not all the ionized deuterium (D/) components 
are resolved, it is often claimed that this value of ^H/H is more dependent 
on the precise descriptio n of the kinematics of the gas and thus less robust 
against systematics (see flPetOSh l 

SDSS 1558-0031, with the absorber placed at z = 2.70262 with metaUicity 
[0/H] = -1.4 9, yielding ^R/E= (3.3ltg;^^) x lO'^ (log^H/H = -4.48 ± 
0.06) flOMeOd l 

Q0913+072 shows six well-resolved D/ Lyman series tra nsition s placed at 
^ = 2.61843 and recently observed with the ESO VLT (|Pet08l l With an 
oxygen abundance of about 1/250 of the solar value the authors of Ref. 
( Pet 08 ) deduce a value of the deuterium abundance ^H/H = x 
10-5 (log^H/H = -4.56 ± 0.04). 

Q0014+81 3, with t he abs orber placed at z = 3.32, yielding ^ H/H ~ (1 .9± 
0.5) X 10"^ f lRug96l : [Cmil . Very old and no more used after (jBur99al ). 
Q0420-388, with the absorber placed at ^ = 3.08, implying, if the deu- 
terium identification is correct, that the ^H/H ratio could have any value 
< 2 X 10"^. Whereas, if the Of/ Hf ratio is constant throughout the complex, 
then 2H/H ~ 2 x 10"^ f lCar96h . 

BR1202-0725, two Lyman- g system s placed at z = 4.383 and z = 4.672, 
yielding ^H/H < 1.5 x 10"^ f lWam96h . 

PG1 718+4807, with the absorber placed at ^ = 0.701 ^H/H G (1.8-3.1) x 
10~^ (IWeb97l : lTyt99l : iKirOll : ICri03l ). but having been criticized in past by 
Ref. f lKirOlh . 

Q0130-4021, with the absorber placed at z 2.8 with metaUi city [Si/H]< 
-2.6, yielding ^H/H < 6.7 x 10"^ (log^H/H < -4.17) flKirOOh . considered 
not very interesting and thus typically neglected. 



Recently, three different analysis of high-redshift systems appeared: 



1) in (jOMeOel; ISteOTl) quasars used are: HS 0105+1619, Q2206-199, PKS1937- 
1009 (I), Q1009+2956, Q1243+3047, SDSS 1558-0031, and the cumulative 
value reported is 

^H/H = (2.681°:^ X 10-^ 
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Fig. 6. Hammer- AitofF projection in Galactic coordinates of the positions of the 
eight quasars along which nine measurements of deuterium abundance have been 
reported (see i)-ix), in the text). 



when averaging over ^H/H determinations (jSteOTI ). or rather 

log^H/H = -4.55 ± 0.04 H/H = (2.82^°;^^) x 10-^ 

when the log's are use d; 

2) the authors of fiFieOTi ) use: HS 0105+1619, Q2206-199, PKS1937-1009 (I), 
Q1009+2956, Q1243+3047 and the cumulative value reported is 

^H/H = (2.78 ±0.29) x 10~^ 

averaging over the log[^H/H] determin ations; 

3) finally, in the very recent analysis of (jPet08l ) the following sample is ex- 
ploited: HS 0105-M619, Q0913+072, Q1009+2956, Q1243+3047, SDSS 1558- 
0031, PKS1937-1009 (I), Q2206-199. An average consistent with all the data 
is logpH/H] = -4.55 ± 0.03 or, equivalently, 

^H/H = (2.82+°:?°) X 10^^ 

The positions in the sky of the reported measurements i)-ix), are given in 
galactic coordinates in Figure 6, and in redshift space in Figure 7. We have 
performed a re-analysis of ^H/H using these results, a veraging over the values 
of logpH/H] and using the method described in Ref. ( Bar04al ) which has the 
particular advantage to give a continuous function even in the case of asym- 
metric errors. In particular, let us denote with Xj ^°^L the generic measurement 
of log[^H/H] corresponding to the i-th QSA. According to Ref. (jBar04al ) one 
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Fig. 7. The nine measurements of i)-ix) QSA's used in our analysis. The horizontal 
band represents the value of Eq. 



can define the following quantities: 



(66) 



and define a total = J^iXlifJ')^ where each contribution xf is expanded 

up to A'^ terms 

X.%) = (^)'(l-2A.(^)+5A^ (^)') . (67) 



Using this procedure we find a value of the reduced \Xmm/i^~ 1) = 
2.715, which shows the effect of some systematic effects and that one or more 
uncertainties have been underestimated. If one chooses to treat all the data 
on the same footing, one can account for this by inflating each uncertainty by 
the multiplicative factor 2.715. In this case, after repeating the procedure, the 
new minimization leads to the result 

log^H/H = -4.53 ± 0.04 H/H = (2.98^°;|) x 10"^ . (68) 



In (jPetOSf) it was argued that the determinations v) and vii) of our list are 
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less robust against systematics in the modeling of the cloud, due to the mi- 
nor number of resolved deuterium lines. To be more conservative, one can 
thus exclude these two data from the regression. In thi s case, with the choice 



for central values and symmetric error bars done in (iPetOSl ). we reproduce 
their results. However, we find that it is not irrelevant to take into account 
the asymmetric errors in the regression procedure. When using the published 
asymmetric errors in log[^H/H] and applying the same procedure as above, we 
find a significantly lower multiplicative factor 4.539)/6 = 1.897, show- 

ing indeed that the dispersion of the measurements for this dataset is more 
consistent with a statistical one. To be conservative we multiply the error bars 
by the factor 1.897 and find 

log^H/H = -4.54 ± 0.03 =^'^ H/H = (2.87+°;^^) x 10"^ , (69) 



which is the value we will be using in the following. 

To conclude, it is worth mentioning the recent proposal to use the fluctua- 
tions in the absorption of cosmic microwave background photons by neutral 
gas during the cosmic dark ages, at redshifts z ^ 7 — 200, to reveal the pri- 
mordial deuterium abundance of the Universe. This method is based on the 
strength of the cross-correlation of brightness-temperature fluctuations due to 
resonant absorption of CMB photons in the 21-cm line of neutral hydrogen, 
with those due to resonant absorption of CMB photons in the 92-cm line of 
neutral deuterium. This results to be proportional to the ratio ^H/H fixed 
during BEN. Although technically challenging, thi s mea surement could pro- 



vide the cleanest possible determination of ^H/H (jSig06l ). A difficulty which 
has been pointed out — that may prevent the viability of the method at red- 
shifts when the first UV sources turn on, z < 40 — is that when including Ly/5 
photons in the analysis, the inferred rat io ^H/H would not be constant, but 



depend sensitively on the UV spectrum (jChuOa ). 



4.2 Helium- 3 



Like deuterium, whose primordial yield is extremely sensitive to the baryon 
density parameter, 77, '^He is a crucial test of the standard BEN scenario as well. 
From the observational point of view, several environments are studied in order 
to derive its primordial abundance. Terrestrial determinations yield e.g. the 
ratio ^He/^He ~ 10~ ^ from balloon measurements or ~ 10~® from continental 
rock (jAnd93l : lRoo02l ) . These observations, which show a large spread of values, 
confirm the idea that the terrestrial helium has no cosmological nature. In 
fact, most of it is '^He produced by the radioactive decay of elements such as 
uranium and thorium. No natural radioactive decay produces ^He, hence its 
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observed terrestrial traces can be ascribed to unusual processes such as the 
testing of nuclear weapons or the infusion of extraterrestrial material. 



In the solar system, ^He is measured through the solar wind and meteorites 
flGloOOh . The most accurate value was measured in Jupiter's atmosphere by 
the Galileo Probe (IMahOSl ). These measurements support the idea of the con- 
version of the deuterium initially present in the outer parts of the Sun into ^He 
via nuclear reactions. One can infer that in the ProtoSolar Material (PSM ) out 
of which the Sun formed, ^He/^He = (1. 66^0.05 ) x 10"^ f lMahQSl : iRooOi ) and 
simultaneously '^R/R= (2.6 ±0.7) x 10^ f lMahQSh . Note that in order to trans- 
form the above mentioned ratios ^He/^He into ^He/H one can use the value 
^He/H ~ 0.1 and thus derive the important ratio ^He/^H= 0.6 ±0.2 which 
allows one to constrain electromagnetically decaying particles (see Section 9). 



The previous observation of ^He/^He is compatible with the same measure- 
ment performed in rn e teoritic gases, yielding ^He/^He = (1.5 ± 0.3) x 10~^ 



( lBla72l : lEbe74l : lGei72l : lGei93l ). Further information comes from Local Inter- 



Stellar Medium (LISM). Occasionally, LISM atoms crossing the termination 
shock region that separates the solar system from interstellar space get ion- 
ized. By counting the helium ions in this particular component of the solar 
wind, t he Ulysses spacecraft has measured a '^He/^He ratio of 2.48lQ;g2 ^ 10~^ 
( iGloQSl ). which is not inconsistent with the idea that '^He at our galaxy's lo- 
cation might have grown in the last 4.6 billion years since the birth of the 
Sun. 



Far beyond the LISM, only one spectral transition allows the detection of ^He, 
namely the 3.46 cm spin- flip transition of ^He'*', the analog of the widely used 
21-cm line of hydrogen; this is a powerful tool for the isotope identification, 
as there is no corresponding transition in ^He"*". The emission is quite weak, 
hence ^He has been observed outside the solar system only in a few HII regions 
and Planetary Nebulae (PN) in the Galaxy. The values found in PN result one 
order of magnitude larger than PSM and LISM d etermin ations (for example 
^He/H = (2 - 5) X 10"^ is measured in NGC3242 flBan02h ). confirming a net 
stellar production of ^He in at least some stars. From the expected correlation 
between metallicity of the particular galactic environment and its distance 
from the center of the galaxy, one would expe ct a gra dient in ^He abundance 
versus metallicity and/or distance. In Ref. (jBan02l ) the ^He/H abundance 
ratios are reported for the sample of simple HII regions. 



Figures (8) and (9) report the data as functions of the distance from galactic 
center and metallicity, respectively, together with the two determination from 
PSM and LISM. No significant correlation between the ^He abundance and 
location (or metallicity) in the Galaxy is revealed. The linear fits reported in 
Figures (8) and(9) as black solid lines correspond respectively to 



41 



o 

X 



1) 
X 



3,5 



2.5 



1.5 



0.5 



LISM 




1 


PSM 








8 10 12 
R^, (Ikpc) 



14 



16 



Fig. 8. The 16 purple spots represent the HIT regions studied in Ref. (jBanOd ) versus 
their distance from galactic center. Also the PSM and LISM measurements are 
reported. The black solid line stands for the linea r fit of purple data, whereas the 
red band represents the upper bound obtained in (iBanOd l. 



{^Re/B.) X 10^: 
(^He/H) X 10^ : 



2.194- 0.030 i?ga«(kpc) 
1.910 + 0.014 [0/H] , 



(70) 
(71) 



where [O /H] = log(0/H)-l og(0/H )^. Note that we adopt (O/H)0 = 4.2x10"^ 
( Asp05l ). differently from ( Ban02 ) where log(O/H)0 = 6.3 x lO"'^ was used. 
Qualitatively, these data suggest a remarkable compensation between stellar 
production and destruction of ^He. 



The failure in observing a galactic ^He dependence on metall icity, ty pically 
predicted by a chemical evolution model of the Galaxy (see (jRom03l ) for a 
review) has been referred to as th e "^He p roblem" . However, in the last years 
this subject received new insight ( lEggOGl ) by the study of 3D mixing models 
which seem to reconcile the predictions with the data. In th is scena rio, by 
assuming a more conservative approach, the authors of Ref. (jBan02l ) prefer 
to report an upper limit to the primordial abundance of ^He by using the 
observations of a peculiar galactic HII region. 



^He/H < (1.1 ±0.2) X 10" 



(72) 



This upper bound is reported as the red band in Figures (8) and (9) (ISte07r i. 



42 



o 

X 



X 



3,5 



2.5 



1.5 



0.5 




V 



-0.4 -0.3 -0.2 -0.1 



V 





[0/H] 



0.1 



0.3 



0.4 



Fig. 9. The data of Ref. (jBan02l ) versus metallicity are reported. The black sohd 
Une stands for the hnear fit of p urple data, whereas the red band represents the 
upper bound obtained in (jBanOd ). 



A different approach could be based on the observation that the ratio (^H + 
^He)/H shows a high level of stability during the galactic evolution. This is 
partially suppor ted by observations and chemical evolution mod els (see for 
example f GeiO?! )). In this case, by using the value reported in Ref. f Gei98 ) for 
PSM, namely (^H + ^He)/H = (3.6 ±0.5) x 10~^, and considering it as a good 
estimate for the primordial value one gets 



^He 
IT 



(0.7 ±0.5) X 10" 



(73) 



by using the primordial deuterium abundance discussed in the previous Sec- 
tion, whose upper bound is consistent with the one derived from Eq. (72). 



4-3 Helium- 4 



The post-BBN evolution of ^He can be simply understood in terms of nuclear 
stellar processes which, through successive generations of stars, have burned 
hydrogen into ^He and hea vier ele ments, hence increasing the ^He abundance 
above its primordial value (jSte07l ). Since the history of stellar processing can 
be tagged by measuring the metallicity {Z) of the particular astrophysical 
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environment, the primordial value of ^He mass fraction Yp can be derived by 
extrapolating the 1^— 0/H and 1^— N/H corre l ations to 0/H and N/H— > 0, 
as proposed originally in Ref.s (jPei74l : |Pei76l : |Pag86|). However, heavy ele- 
ments like oxygen are produced by short-lived massive stars whereas ^He is 
essentially synthesized in all stars, so one has to minimize model-dependent 
evolutionary corrections. The key data for inferring ^He primordial abundance 
are provided by observations of helium and hydrogen emission lines generated 
from the recombination of ioniz ed hydrogen and helium in low-metallicity ex- 
tragalactic H/j regions (ISteOTl ). Many attempts to determine Yp have been 
made, constructing these correlations for vari ous sam ples of Dwarf Irregular 
(DIrrs) and Blue Compact Galaxies (BCGs) (jlzoQSl ). These systems are the 
least chemically evolved known galaxies. Plausibly, they contain very little 
helium synthesized in stars after the BEN, minimiz ing the chemical evolution 
problems that affect e.g. the determination of ^He (jlzo04j ) . 



Uncertainties in the determination of Yp can be statistical or systematic. Sta- 
tistical uncertainties can be decreased by obtaining very high signal-to-noise 
ratio spectra of BCGs. These BCGs are undergoing intense bursts of star for- 
mation, giving birth to high-excitation supergiant H// regions and allowing an 
accurate determination of the helium abundance in the ionized gas through 
the BCG emission-line spectra. The theory of nebular emission is understood 
well enough not to introduce additional uncertainty. According to the stan- 
dard scenario, the universe was born with zero metallicity; hence, Yp can be 
determined extrapolating to Z ^ the relationship between Z and the '^He 
abundance for a sample of objects. This procedure relies on the determination 
of the individual Yp and Z values and of the slope dYp/dZ, which is assumed 
to be linear. The uncertainty affecting Yp depends directly on the uncertain- 
ties on dYp/dZ and the ensemble of the {Yp, Z) pairs. For this reason, it has 
long been thought that the best results are obtained from the analysis of 
extremely low metallicity objects like DIrrs and BCGs, since their use mini- 
mizes the uncert ainty associated with dYp/dZ. However, the authors of Ref.s 
f|Pei03al : |Pei03bl ) ^lave noted that this advantage is outweighed by the rela- 
tively higher uncertainty on the l^-values, which derives from the (unknown) 
collisional contribution to the Balmer line intensities, an uncertainty especially 
affecting these objects since collisional contribution is quite important at high 
temper atures and rapidly fades away at intermediate a nd low temperatures 
(ILurOSl ). Mostly following the analysis reported in Ref. (jPeiOTI ) we list below 



the most recent estimates of Yp. 



i) In Ref. (jlzoOj) is reported the estimate Yp = 0.2421 ± 0.0021. According 



to Ref. (iPeiOTI ) the differences with the prev ious de termination are mainly 
systematic. One is due to t he use i n Ref. (jPeiOTl ) of Hei recombination 
coefficients studied in Ref.s (jPorOSl : iPorOTI ). which yield Yp val ues about 
0.0040 higher than the previous ones. Moreover, in Ref. (jPeiOTI ) they use 
some recent Hj colhsional data 



which further increase the Yp values over 
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the older Hj collisional c orrecti ons by about 0.0025. 

The value quoted in Ref. (jOli04l ) is Yp = 0.249±0.009. The small sample size 
used and the large uncertainty affecting the parameters derived from the Hjj 
regions considered in the analysis are responsible in this case for the very 
conse rvative error estimate. Also in this case, the systematic differences with 
Ref. flPeiOTi ) are due to the Hei recombination data used by both groups 
and to the esti mation of the collisional contribution to the H Balmer lines. 
In Ref. flFukOeh . based on a reanalysis of a sample of 33 Hjj regions from Ref. 



(JIzoOJ), the authors determined a value oi Yp = 0.250 ± 0.004. In addition 
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to a different treatment of the underlying H and Her a bsorption there are 
few syst ematic effects discussed in detail in Ref. ( jPeiOTI ). 
In Ref. (jPeiOTl ). the authors present a new ^He mass fraction determina- 
tion, yielding Yp = 0.2477 ± 0.0029. This result is based on new atomic 
physics computations of the recombination coefficients of Hej and of the 
collisional excitation of the Hj Balmer lines together with observations and 
photoionization model s of metal-poor extragalactic Hji regions. 
Finally, in Ref. (jlzo07l ) is reporte d the e stimate Yp = 0.2516 ± 0.0011 when 
using the Hej emissivity of Ref. (jPor05l ) . 



All recent analyses of Yp agree on the fact that the systematic error is the 
main r esponsible for the spread of the Yp determinations. For example, in Ref. 
f lPei07l ) ;hree of the four main sources of error (Al^ > 0.01) in estimating 
Yp are reported to be of systematic nature. Different authors however, report 
different error budgets (sometimes analyzing the same objects) . Nevertheless , 
the use of new Hej recombination coefficients studied in Ref.s (jPor05l : |Por07l ) 
has sensibly ch anged the pr edictions by pushing them up a few percent as 
shown in Ref.s ( lPei07l : Ilzo07l ). For this reason we will use their determinations 
only to derive an estimate of primordial Yp. 



In particular, in Ref. (jPei07l ) the uncertainty quoted on the value of Yp seems 
to provide a more realistic estimate (0.003) of the residual indetermination 
affect ing th e ^He mass fraction, versus the more optimistic (0.001) reported in 
Ref. fllzo07l ). due to a very large sample of objects included in the analysis. For 
this reason we prefer to average the two main values as not weighte d by th eir 
uncertainties, and use the more conservative error estimate of Ref. (|Pei07l ) at 
1-0". In summary, we adopt for the ^He mass fraction 



K = 0.250 ± 0.003 



(74) 



Note that this range of values totally contai ns the one of Ref. (jlzo07l ). but 
only partially covers the determination of Ref. (jPei07l ) which points out slightly 
smaller values. Hence the effects of new Hej recombination coefficients jPor05 : 



Por07l ) on global Yp data analysis is to push slightly upward its estimate, and 



this is going to affect at a certain level the determination of cosmological 
parameters, as will be discussed in the following sections. 
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We would like to briefly point out that o ther co nstraints on Yp can be obtained 
by indirect methods. For example, in (jSalOJ) Yp was bounded from studies 
of Galactic Globular Clusters (GGC). The value they found, Yp < Yggc = 
0.250 ± (0.006)stat =t (0.019)sys, is consistent with the above estimate. Finally, 
CMB anisotropies are sensitive to the reionization history, and thus to the 
fraction of baryons in the form of ^He. Present data only allow a marginal 
detection of a non-zero Yp, and even with PLANCK the error bars from CMB 
will be l arger than the present systematic spr ead of the astrophysical determi- 
nations (ITroOSl : lHue04l : llchOa : llchOTI : iHamOSl ). On the other hand, imposing a 
self-consistent BEN prior on Yp would improve the diagnostic power of CMB 
data on other parameters, thus representing another nice synergy of CMB and 
BBN, besides the concordance test provided by 77. 



4-4 Lithium- 7 



Lithium's two stable isotopes, ^Li and ''Li, continue to puzzle astrophysicists 
and cosmologists who try to reconcile their primordial abundance as inferred 
from observations with the BBN predictions. From the astrophysical point of 
view the questions mainly concern the observation of lithium in cold interstel- 
lar gas and in all type of star s in which lithium lines are either detected or 



potentially detectable (lAspOGl ) . 



A chance to link prim ordial ^ Li with the BBN abundance was flrst proposed by 



Spite & Spite (1982) (jSpi82l ). who showed that the lithium abundance in the 
warmest metal-poor dwarfs was independent of metallicity for [Fe/H]< —1.5. 
The constant lithium abundance deflning what is commonly called the Spite 
plateau suggested that this may be the lithium abundance in pre-Galactic gas 
provided by the BBN. The very metal-poor stars in the halo of the Galaxy 
or in similarly metal-poor GGC thus represent ideal targets for probing the 
primordial abundance of lithium. Even though lithium is easily destroyed in 
the hot interiors of stars, theoretical expectations supported by observational 
data suggest that although lithium may have been depleted in many stars, t he 



overall trend is that its galactic abundance has increased with time (ISteOTI ). 



There is quite a long tale of ''Li determinations appearing in the literature, 
starting from the Spite & Spite (1982) value of rLi/H]= 2.05 ± 0.15 (by 
deflnition [^Li/H]= 12 + log;^Q('''Li/H)). For the sake of brevity we will restrict 
our analysis roughly to the determinations of the last decade. The implicit 
assumption is that, hopefully, the more recent papers have reached a better 
understanding of the systematics involved in the inference of the primordial 
^Li abundance. 



In Ref.s (lRya99l : iRyaOOl ) a study of a set of very metal-poor stars showed a 
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very small intrinsic dispersion in the ^Li abundance determinations. More- 
over, the authors found evidence for a decreasing trend in the ^Li abundance 
toward lower metallicity indicating that the primordial abundance of ^Li can 
be inferred only after allowing for nucleosynthesis processes that must have 
been at work in the early stages of the Galaxy. The primordial ^Li abun- 
dance reported is [^Li/H]= 2m±toil (^Li/H= (l.23l[J:|) x 10~^°)- Different 
studies of halo and G GC star s provide d higher lithium plateau abundance 
rLi/ H]= 2.2 4 ± 0.01 flBon97al : lBon97bf ). A similar analysis is contained in 



Ref. flBon02l ). where high resolution, high signal-to- noise ratio spectra of 12 
turn-off stars in the metal-poor globular cluster NGC 6397 were used. The 
author conclude that, within the errors, they all have the same lithium abun- 
dance rLi/H]= 2.34 ±0.06. 



In Ref. ( lMel04l ). a study of ^Li abundance in 62 halo dwarfs was performed by 
using accurate equivalent widths and a temperature scale from an improved 
infrared flux method. For 41 plateau stars (those with Teff > 6000 K) the 
^Li abundance is found to be independent of temperature and metallicity, 
with a star-to-star scatter of only 0.06 dex over a broad range of tempera- 
tures (6000 K < Teff < 6800 K) and metallicities (-3.4 <[Fe/H]< -1). Thus 
they rep ort a mean ^Li plateau abundance of [''Li/H]= 2.37 ± 0.05. In Ref. 



flChaOSh ;he authors underwent a very detailed reanalysis of available obser- 



vations; by means of a careful treatment of systematic uncertainties and of 
the error budget, they find [''Li/H]= 2.21 ± 0.09 for their full sample and 
[^Li/H]= 2.18±0.07 for an analysis restricted to unevolved (dwarf) stars only. 
They also argued that no convincing/ conclusive evidence for a correlation be- 
tween ^ Li and m etallicity can be claimed at present. More recently, the authors 



of Ref. (lAspOd ) have studied a set of 24 very high quality spectra metal-poor 
halo dwarfs and subgiants, acquired with ESOs VLT/UVES. The derived one- 
dimensional, non-Local Thermodynamical Equilibrium (non-LTE) ^Li abun- 
dances from the Lij 670.8 nm line reveal a pronounced dependence on metal- 
licity but with negligible scatter around this trend. The estimated primordial 
^Li a bundan ce is ^Li/H e (1.1- 1.5) x 10"^° ([^Li/H]= 2.095±0.055). Recently 



Ref. (jKor06l ) has reported the spectroscopic observations of stars in the metal 
poor globular cluster NGC6397 that reveal trends of atmospheric abundance 
with evolutionary stage for various elements. These element-specific trends 
are reproduced by stellar-evolution models with diffusion and turbulent mix- 
ing. They compare their observations of lithium and iron to models of stellar 
diffusion, finding evidence that both lithium and iron have settled out of the 
atmospheres of these old stars. Applying their stellar models to the data they 
infer for the unevolved abundances, [Fe/H] = 2.1 and [''Li/H] = 2.54 ± 0.10. 

The list of the last ten years estimates for ''Li abundance is then the following: 



Li/H] = 2.24 ± 0.0 1 flBon97al: lBon97bh : 



1) \'Li/ti\ = 2.2^ ±o.o L imori^/ai: moi 
ii) fLi/H] = 2.09l[5:l| flRva99URvao'oF ) 
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fLi/H] = 2.34 ± 0.06 f Bon02h : 
f Li/H] = 2.37 ± 0.05 (Mdol; 
[^Li/H] = 2.21 ± 0.09 flCha05h: 
f Li/H] = 2.095 ± 0.05 5 flAspOeh : 
f Li/H] = 2.54 ± 0.10 flKorOOh . 



All but (marginally) the latter value are inconsistent with standard BBN pre- 
dictions for the preferred range of r] singled out by CMB data, which fits 
remarkably well deuterium abundance. It is unclear how to combine the dif- 
ferent determinations in a single estimate, or if the value measured is truly 
indicative of a primordial yield. A conservative approach (similar to the one 
used for ^He) is to quote the simple (un-weighted) average and half-width of 
the above distribution of data as best estimate of the average and "systematic" 
error on ^Li/H, obtaining 



If 



2.27 ±0.23 




1.86 



+1.30 
-1.10 



X 10 



-10 



(75) 



Note that the statistical error is much smaller (of the order of 0.01), but we 
will not need it since, due to its uncertain status as tracer of the primordial 
value, ^Li is typically excluded in "conservative" BBN statistical analyses (or 
rather invoked to support particular non-standard BBN scenarios). 



4-5 Lithium-6 and "The lithium problems 



It is clear from the above discussion and from the substantial disagreement 
of almost all of the ^Li observations with the standard BBN predicted value 
(by about 0.4 dex, assuming the central value of the quoted average) that 
some piece of (astro)physics is missing. For a detailed discussio n of pos sible 
cause s we add ress the reader to the excellent review given in Ref. flAspOd ) see 
also ( Lam05al )). Here we want to remark that: i) a ~ 1.5 -i- 2 lower value of 
T] at the BBN time with respect to the best fit deduced from CMB data is 
excluded by the agreement between deuterium observations and CMB value 
of ?7, but also by the inferred upper limit of the primordial ^He abundance; ii) 
underestimated errors in the adopted nuclear reaction rates are now ex cluded : 
the lab o ratory r aeasure ments of the crucial ^He(a, 7) '^Be cross-sectio n (Nar04l : 



Bem06l : ICon07l : lBro07l ). its inferred rate from solar neutrino data (ICyb03al ) 



and the meas uremen t of the proposed alternative channel for ^Be destruction 
^Be((i, p)2a fjAng05l ) all point to the conclusion that nuclear uncertainties 
cannot explain the discrepancy between observed and predicted primordial 
^Li abundances; iii) systematic errors in the abundance analysis, although 
in principle still possible, seem very unlikely. The introduction of 3D model 
atmospheres, even accounting for non-LTE, has not resulted in a significant 
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upward revision of the lithium abundance obtained from more primitive ID 
atmospheres. 



The most hkely causes of the "problem" are thus: a) either some modification 
to the BBN scenario; b) or, perhaps more likely, that the lithium abundance 
of very metal-poor stars is not the one of the primordial gas. We will illus- 
trate later some scenarios of the type a). Here, however, we want to point 
out that explanations of the type b) are probably not trivial and might in- 
volve both reprocessing in situ (the observed stars) and earlier lithium syn- 
thesis/depletion in the pre-galactic, young universe environment. Indeed, even 
assuming that s ome di ffusion and turbulent mixing mechanism like the one 
pointed out in (IKorOGl ) can explain the ''Li problem, still an issue remains 
with ^Li. The presence of the fragile ^Li isotope, which is produced during the 
BBN at the level of ^Li/H ~ 10~^^ — 10~^^, has been recently confirmed in a 
few metal-poor halo stars, with some hint of a pl ateau vs. metallicit^f^ with 
abundance as high as ^Li/H ~ 6 x 10~^^ flAspOGi ). These data (at least par- 
tiall y) confirm the first observations report ed in literature already f ifteen y ears 
ago flSmi93bl : ISmiQSl : lHob97l : ICavQQl : iNisOoh . In particular, in Ref. flAspOeh ^Li 
is detected in 9 of 24 analyzed stars at the > 2a significance level. How- 
ever, it is worth str essing that these observations have been questioned in 
Ref. (jCayOTI : ICayOSi ). According to these papers the convective asymmetry 
generates an excess absorption in the red wing of the ''Li absorption feature 
that could mimic the presence of ^Li at a level comparable with published 
values. This would mean that only an upper limit on ^Li/''Li can be derived 
at present. 



Both ^Li, ^Li can be produced by fusion {a + a — >Li) and spallation (p+CNO 
— > LiBeB) reactions by ordinary cosmic ray primaries impinging on nucleons 
and nuclei in the intergalactic medium (see e.g. (I VanOOl ) ) . Additionally, obser- 
vations are performed not in "inert" gas clouds, but in stellar atmospheres, 
where thermonuclear burning depletion (of particularly fragile nuclei) is a 
crucial effect. It follows that primordial production mechanisms and later as- 
trophysical effects might be competing. Viable astrophysical candidates for 
the acceleration of cosmic rays, able to account for the ^Li observe d at low 
metallicity include: the massive black hole in the Galactic center (jPraOSi ). 
radio-l oud AG Ns (one of which c ould have been present in our Galaxy in the 
past) (INatOGl ). a nd Pop III stars (IRolOa : IRoIOGI ) which may also explain the 
depletion of ^Li (jPiaOGi ). but have been rece ntly ch allenged on the light of 
further constraints adopted in the model of flEvoOSl ) . The shocks developed 
during structure formation (jinoOSi ) . which however conflict with astrophysi- 



^ However, taking into account predictions for ^Li destruction during the pre-main 
sequence evolution tilts the plateau suggesting a ^Li increase with metallicity. Basi- 
cally, fairly uncertain stellar pre-main-sequence destruction of ^Li could be respon- 
sible for an apparent plateau 



RicOd : iRicOSi : lAspOa l. 
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cal constraints (jProOTI ). are not powerful enough, since in the early times of 
Galaxy formation, the masses of the assembling dark h aloes w ere still quite 
small and the corresponding virial velocities insufficient (jPraOSi ) . 



Perhaps, the most convincing explanation proposed until now is that "^Li may 
be produc ed in s itu from stellar flares within the first billion years of the 
star's life ( TatO?! ). In particular, the anomalously high ^He/^He ratio found 



in solar fiares (~0.5) and the kinematically much more favorable process 
^He('^He,p)^Li compared to aa — > ''Li fusion reactions provide a physical 
mechanism for producing large quantities of ^Li without overproducing ^Li. 
The main issue with this or other scenarios might be the difficulty to test 
them. 

A great help in solving this i ssue m ight come from detecting lithium in a 



different environment. In Ref. (jPro04l ). it was proposed to independently test 
the pre-Galactic Li abundance by looking at high velocity gas clouds falling 
onto our Galaxy, with metallicities as low as 10% of the solar one. If these 
low-metallicity clouds have a mostly pre-Galactic composition, with a small 
contamination from the Galaxy, they might allow probing the lithium abun- 
dance at least free of the possibility of thermonuclear depletion in situ. Another 
proposal to detect the cosmological reco mbinati on of lithium via its effect on 



the mi crowav e background anisotropies (jSta02l ) has been proved not to be 



viable (jSwiOSi ). 



We will not use ^Li to derive constraints in this review. However, using the 
observed ^Li as an upper limit to its primordial value turns out to be a pow- 
erful constraint in some regions of parameter space for exotic models. While 
reporting them in the following, we warn the reader that their robustness relies 
on the assumption that no destruction or major reprocessing of the ''Li and 
^Li observed in the halo stars has happened, which might be overly optimistic. 



5 Standard BBN theoretical predictions versus data 



The goal of a theoretical analysis of BBN is to obtain a reliable estimate of 
the model parameters, once the experimental data on primordial abundances 
are known. In this Section we will consider only the case of the standard BBN, 
where the only two free parameters are the value of the baryon energy density 
parameter Qsh'^ (or equivalently the baryon to photon number density, rj) and 
possibly, a non-standard value for the relativistic energy content during the 
BBN. The latter, after annihilation can be parameterized in terms of the 
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Fig. 10. Values of the primordial abundances as a function of ?/io, calculated for 
AA'eff = 0. The hatched blue bands represent the experimental determination with 
1 — 0" statistical errors o n Yp, ^H, and ^Li, while the red band is the upper bound 
obtained in Ref. (jBannd ). Note that for high value of rjio all ^Li comes from ^Be ra- 
dioactive decay via electron capture. The vertical green band represents the WMAP 
5-year result O^/i^ = 0.02273 ± 0.00062 dPunOSl V 



effective number of neutrinos we liave recalled in Section 3.2, 



(76) 



Similar analyses have been recently presented by various groups, which might 
be slightly different depending on t he adopted values of and/or ^H/H ex- 



Cvb04: 


Cvb05: 


Han02: 


Man07: 



In the minimal scenario the parameters reduces to the baryon density only, 
since AA^eff is assumed to vanish. Figure 10 shows the dependence on 7710 of 
the final value of the primordial yields, calculated using PArthENoPE, along 
with the experimental values of the abundances and their corresponding un- 
certainties, as discussed in the previous Section. 
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Fig. 11. Likelihood functions for ^H/H (narrow) and Yp (broad). 
To get confidence intervals for 77, one constructs a likelihood function 

£(r/)(xexp(-x'(^)/2) , (77) 

with 

x\v) = T.[Mv)-xt']wM[xAv)-xf'] . (78) 

The proportionality constant can be obtained by requiring normalization to 
unity, and Wiji^r]) denotes the inverse covariance matrix, 

where CTjj and cxj exp represent the nuclear rate uncertainties and experimental 
uncertainties of nuclide abundance Xi, respectively (we use the nuclear rate 
uncertainties as in Ref. ( Ser04bl )). while by (rfj^other denote the propagated 
squared error matrix due to all other input parameter uncertainties (r„, Gn, 
etc.). We use the following values for the experimental yields of and '^He 
(see previous Section): 



^H/H = (2.87toil) X 10-^ Yp = 0.250 ± 0.003 . (80) 

We first consider abundance alone, to illustrate the role of deuterium as 
an excellent baryometer. In this case the best fit values found are flsh^ = 
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0.021 ± 0.001 (r/io = 5.7 ± 0.3) at 68% C.L. , and ^Bh'^ = 0.021 ± 0.002 at 
95% C.L. . A similar analysis can be performed using "^He. In this case we get 
neh^ = 0.028tg:l]J^ (r/io = r.QtHjB at 68% C.L. , and ^Bh^ = 0.028^^?^ at 
95% C.L. . Figure 11 shows, as from our discussion in the previous Section, that 
the determination oiflsh^ is mainly dominated by deuterium. In any case, the 
result is compat ible at 2-a with the WMAP 5-year result Qsh"^ = 0.02273 ± 



0.00062 (jDunOSi). The slight disagreement might have some impact on the 



determination from CMB anisotr qpies o f the primordial scalar perturbation 



spectral index Tig , dbS noticed in (jPet08l ). where the BBN determination of 
flBh^ from deuterium is used as a prior in the analysis of the five year data of 
WMAP. 

In Table 11 we report the values of relevant abundances for some different 



baryon densities, evaluated using PArthENoPE fjPar08l ). Notice the very low 
prediction for ®Li (see discussion in Section 4) and that, for these values of 
baryon density, almost all ^Li is produced by ''Be via its eventual electron 
capture process. 

If one relaxes the hypothesis of a standard number of relativistic degrees of 
freedom, it is possible to obtain bounds on the largest (or smallest) amount 
of radiation present at the BBN epoch, in the form of decoupled relativistic 
particles, or non-standard features of active neutrinos (but see our previous 
discussion in Section 3.2). Figure 12 displays the contour plots 68% and 95% 
C.L. of the total likelihood function, in the plane {flBh'^,Ncs). After marginal- 
ization one gets fi^/i^ = 0.021 ± 0.001 and N^s = ^-^StoH at 68% C.L., and 
n^h'^ = 0.021 ± 0.002 and N^s = 3.18^^^ at 95% C.L. . Hence the global 
analysis results to be compatible with iVeff = 3.046 and fis/i^ = 0.02273 found 
by WMAP at 1-a level (see Figure 12). 



^ Note that r/io reported in this section is conventionahy defined, according to 
PArthENoPE , as 7710 = 273.491^^/1^ which is sHghtly different from the definition of 
Eq. (9), but simply connected to it. 

Table 11 

The theoretical values of the nuclear abundances for some value of VtBh?- 



^Bh^ 


0.017 


0.019 


0.021 


0.023 


0.028 




0.245 


0.246 


0.247 


0.248 


0.250 


2H/H xlO^ 


4.00 


3.36 


2.87 


2.48 


1.79 


3He/H xlO^ 


1.22 


1.14 


1.07 


1.01 


0.903 


^Li/H xlO^° 


2.53 


3.22 


3.99 


4.83 


7.08 


■^Be/H xlO^o 


2.15 


2.89 


3.69 


4.56 


6.88 


^Li/H xlO^^ 


1.72 


1.45 


1.25 


1.08 


0.791 
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Fig. 12. Contours at 68 and 95 % C.L. of the total likelihood function for deuterium 
and "^He in the plane {Qsh'^ ,^es)- The bands show the 95% C.L. regions from 
deuterium (almost vertical) and Helium-4 (horizontal). The red cross corresponds 
to the standard N^g and Q.Bh? = 0.02273 as indicated by WMAP 5-year results. 

6 BBN and Neutrino physics 

We have already stressed how large is the impact of neutrino physics and neu- 
trino properties on BBN. In fact, the discovery of neutrino masses via oscilla- 
tions, combined with the stringent bounds on the effective electron neutrino 
mass via tritium decay experiments, has had a profound impact on the phe- 
nomenology of active neutrinos in the early universe. At the moment we know 
that: (at least two) neutrinos are massive; all the masses are small {mi, < eV, 
possibly much smaller); individual lepton numbers are violated (with mixing 
angles much larger than in the quark sector), although it is not known if the 
overall lepton number, L, is conserved. We do not know yet if the neutrino 
mass term in the Lagrangian is of the Dirac (~ [Ai + h.c], flavor indexes 
omitted) or Majorana (~ [M. vIvr + /i.e.], where v1 is the charge conjugate 
field) type, but in either case new physics is required. In the Dirac case, one is 
forced to introduce the yet undetected right-handed neutrino fields, vr. In the 
Majorana case, one assumes (differently from the SMPP) that the lepton num- 
ber L is violated and introduces a Majorana mass operator, which is allowed 
for neutrinos, being the only neutral fermions in the SMPP. The important 
news for BBN are that even the incomplete knowledge of the mass matrix 
M. that we have at present is enough to conclude that the phenomenology 
of active neutrinos in the BBN is greatly simplified. A plethora of cases once 
popular in the literature are now excluded. Among the ones which were of 
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remarkable interest only a decade ago we can mention: (i) a lower-than-three 
effective number of neutrinos due to a "massive z/^-" (improper, not being a 
mass eigenstate) ; (ii) a decaying 'Vr"; (iii) the thermalization of right-handed 
neutrinos (for the Dirac mass case), which is inhibited by the smallness of the 
neutrino masses by which they are coupled to the active states. We do not 
treat these issues further and ad dress the interested reader to these histori- 
cal topics to the review (jDol02bl ). In the following, we focus on the bounds 
on electromagnetic interactions of neutrinos in Section 6.1, while Section 6.2 
treats other exotic interactions. The topic of neutrino asymmetry is briefly 
reviewed in Section 6.3, while Section 6.4 treats the impact on BBN of sterile 
neutrino states. 



6.1 Bounds on electromagnetic interactions of neutrinos 



Dropping flavor indexes, the most general structure of effective neutrino elec- 
tromagnetic interactions is 



^int = -ey u-f^uA^'-a^u-ff.-f^u dxF^'^ - \ a^p (/i + e 75) z/F"^ (81) 



where F"^ is the electromagnetic field tensor, = [7^,7/3], and the form 
factors {e,^, a^, /i, e}, which are functions of the transferred squared momentum 
g^, in the limit ^ correspond to the electric charge, anapole moment, 
magnetic and electric dipole moment, respectively. 



In principle, Dirac neutrinos may have a very small electric charge e^. BBN 
bounds may be derived by requiring both that right-handed partners are not 
populated and that neutrinos are not kept in equilibrium too long after the 
weak freeze-out, which would alter the photon-neutrino temperature relation. 
However, for the range of masses presently allowed, the BBN bounds are never 
competitive with other astrophysical or laboratory constraints, as for instance 
the red giant bound of < 2 x 10^^^ (lRaf99l ) . Actually, the ind irect b ound 
comi ng frorn the neutrality of matter is stronger [e^ < 10~^^, (jFoo90l ) : see 
also f DavOOl ) and reference therein for details), so we ignore in the following 
a possible neutrino charge. 
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The possibility of a neutrino charge radiua ^° I (which can be negative), 



, ^ 6 / de,{q^) 



m) 



has been considered in the hterature. After a long debate, it has been finall 
establ i shed tha t (r^) is a well-defined (gauge- independent) quantity (IBerO 
Ber02l : lBer04al ). The Standard Model expectations are in the range (r^) ^ 
1 4nb. Differently from the case of the magnetic moment, the charge radius 
does not couple neutrinos to on-shell photons, so stellar cooling arguments are 
not very sensitive to (r^). Yet, for Dirac neutrinos the channel e~e^ — >■ i'ri^r is 
still effective, provided that new physics (NP) violates the cancelation between 
vector and axial contribution that otherwise applies in the SMPP. If this 
cancela tion does not take pl ace, th e corresponding bounds from SN 1987 A 
( lGri89l ) and nucleosynthesis (lGri87l ) read 



I ) 1 ^ 2 (7) nb, from SN1987A (BEN) 



(83) 



For Majorana neutrinos, the previous bounds do not apply. In this case, how- 
ever, even in the z/T--sector where BEN may have a sensitivity comparable to or 
better than laboratory experiments, even a change of one order of magnitude 
above the SMPP level in the channel e~e~^ —>■ VrVr due to new physics would 
only bring changes at the 0.1% level in 1^, so that only labora tory bounds are 



meaningful. For a further discussion of this point, see (jHir03l ). 



Another consequence of the existence of neutrino masses is that a neutrino 
magnetic moment is naturally present (and it is usually expressed in terms 
of Eohr magnetons, /Ub), although for Majorana particles only off-diagonal 
elements in flavor space are non-vanishing. There are two possible processes of 
interest for EEN: (i) the thermalization of Dirac neutrinos via e.g. the process 
^ e; (ii) a radiative decay of the kind v]^ — > z/^, 7^ ^^e latter being 
possible also for Majorana neutrinos. In absence of primordial magnetic fields, 
the EEN bound on the diagonal elements coming from the thermalization of 
right-handed neutrinos is not as restrictive as the one coming from red gia nt 



cooling argument (from plasmon decay 7* — > vv, yU < 3 x 10 ^"^ flRaf99l ) ) 
The radiative decay rate for a transition z — > j is 




In complete analogy, one can define an anapole radius (r^) from a^. For Majorana 
neutrinos, symmetries require some of the e.m. couplings to vanish. Since the astro- 
physical/cosmological bounds do not typically distinguish between (r^) and (r^), or 
electric and dipole magnetic moments, we shall quote bounds on (r^) and /i in this 
loose sense. 
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and typical bounds lead to a lifetime definitively too long to affect the BBN 
cos mology . In any case, the very cosmological bounds on the neutrino lifetime 
in (jMirOTI ) exclude any effect of the radiative neutrino decay at the BBN 
epoch. However, in presence of very strong primordial magnetic fi elds, th e 
BBN bound via spin-precession may be as str ingent as yU < 



10' 



Although very model-dependent (see (lEnq95l ) . in particular Eq.s (50,51)), this 
is to our knowledge the only bound probing the level of intensity expected for 
the dipole moment in the SMPP enlarged with a Dirac neutrino mass term, 
which is = 3eG'j.m^/(8y27r2) = 32 x lO'^o {mJeV)i^B- 



6.2 Bounds on other exotic interactions of neutrinos 



Besides anomalous electromagnetic interactions, neutrin os might undergo non- 
standard-interactions (NSI) with electrons. Ref. (IManOGl ) considered low-energy 
four-fermions interactions of the kind 



''-NSI 



^5) 



with the NSI parameters constrained by laboratory measurements to 

be at most of OiV). It was found that, for NSI parameters within the ranges 
allowed by present laboratory data, non-standard neutrino-electron interac- 
tions do not essentially modify the density of relic neutrinos nor the bounds 
on neutrino properties from cosmological observables. Qualitatively, this de- 
pends on the fact that a large modification of the neutrino spectra would only 
be achieved if the decoupling temperature were brought below the electron 
mass. The presence of neutrino-electron NSI within laboratory bounds may 
enhance the entropy transfer from electron-positron pairs into neutrinos, up 
to a value of A^eff = 3.12 (and AYp ~ 6 x 10"^), which are almost three times 
the corrections due to non-thermal distortions that appear for standard weak 
interactions, but still probably too small to be detectable in the near future, 
even for PLANCK. 

Another scenario is the one where neutrinos couple to a scalar or pseudoscalar 
particle with a Yukawa-type interaction of the kind 

L,^ = \d^,c\>d^^-\rn\<\?-<\>Y^v'^\,,vi , (86) 



where j are flavor indexes and, in case of pseudoscalar coupling, A^j — > 75A£j. 
We shall denote these couplings simply as A if we ignore flavor effects and only 
refer to constrai nts within a factor of 0(1). A famous case of this kind is the 
Majoron model (jChiSlI ; iGelSlI : iGeoSlI ) where the Majoron is the Goldstone 
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boson associated to the breaking of the lepton number symmetry (and thus 
T^tp ~^ 0). Although laboratory bounds rule out the original model as an 
explanation of the small neutrino masses, still it represents a prototype of "se- 
cret neutrino interactions", where neutrinos interact with a sector precluded 
to other standard model particles, and may thus have stronger interactions 
among themselves at low energies, than predicted by the SMPP. In the early 
universe, a large enough A would allow t o populate thermally the species 0. 



For -C 1 MeV, it was found in flCha94l ) that A < 10 , if one considers one 



additional boson (AA^^efr = 3/7) to be incompatible with the observations. If 
this is instead considered viable, the same process uu ^ cfxf) may be respon- 
sible for a "neutr inoless universe" well after the BBN epoch, provided that 



<C (lBea04j ) . 



Relatively less attention has been paid to the case where the associated (pseudo) 
boson is massive (although a particle of this kind might have other "cosmo- 



logi cal virt ues" , as providing a warm dark matter candidate, see e.g. f lLatOTl )) 



In fICuoOSi ) the authors assumed MeV-scale (p particles produced at early 
epochs via additional couplings with other SMPP particles and later decay- 
ing into neutrinos in out-of-equilibrium conditions (with a rate r(</) — > up) = 
3A^m0/(87r)); when this happens before the photon last scattering epoch, the 
produced neutrino burst directly influences the CMB anisotropy spectrum, as 
well as the late LSS formation. They show that current cosmological obser- 
vations of light element abundances. Cosmic Microwave Background (CMB) 
anisotropies, and Large Scale Structures (LSS) are compatible with very large 
deviations from the standard picture. They also calculate the bounds on non- 
thermal distortions which can be expected from future observations, finding 
that the present situation is likely to persist with future CMB and LSS data 
alone. On the other hand, the degeneracy affecting CMB and LSS data could 
be removed by additional constraints from primordial nucleosynthesis or in- 
dependent neutrino mass scale measurements. 



In (jHanOJ), the particle was considered to be the inflaton, only coupled 
to neutrinos and with a mass >> 1 MeV, to determine via BBN and 
other cosmological observation the constraint on the lowest possible reheating 
temperature Tr/^. In particular, the author derives constraints from partial 
thermalization as well as neutrino spectral distortions. Barring fine-tuning, 
the res ulting bound is Trh ^ 4 MeV. A factor ~ 2 lower bound was found in 
fllchOSi ) if no (j) boson is included but oscillations are taken into account. It is 



worth noticing that in scenarios with late-time entropy production a constraint 
arises anyway due to the incomplete background neutrino thermalizati on, even 



if the n eutri nos do not couple to (f) directly, as already noted in ([Kaw99 



KawOOl ) . In (jSer04al ). light scalar particles annihilating into neutrinos were 
considered to analyze the effect on BBN of MeV-scale dark matter particles, 
invoked to explain the exc ess of 5 11 keV photons from positron annihilation 



from the Galactic Center ((BoeOJ). If such particles only couple to neutrinos. 
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they need to be heavier than ~ 1 MeV to be consistent with the Yp constraints. 
If they have an additional couphng to e"*" e~ at the level required to explain 
the Galactic Center positrons, the bounds may be more stringent (but then 
depend on the details of the model). 

We have seen that a Dirac mass term could in principle be responsible for the 
production of right-handed neutrinos in the primordial plasma. While this 
possibility is excluded by the smallness of neutrino masses, it is still possible 
to populate uji via direct right-handed currents mediated by Wr bosons, of 
the kind i^rH^rI^r (or analogous coupling with right-handed charged leptons). 
These are possible in some extensions of the standard electro-weak model. If 
one assumes that the right-handed interaction has the same form as the left- 
handed one but with heavier intermediate bosons, one can obtain from BBN a 
lower hmit on their mass of the order of ^ 75mw ( Ste79 : OliOd : BarOSb), 



which depends however on the exact particle spectrum of the physics beyond 
the SMPP up to ~ 75 mw 



6.3 Neutrino asymmetry 



The origin of the most fundamental parameter in BBN, the baryon asymmetry 
= {ns — nR)/n^ (or simply rj at late times), is not known. While SMPP 
and SMC contain all the ingredients required to generate it dynamically from 
an initially symmetric universe (B, C, and CP violating interactions, depar- 
ture from thermal equilibrium) (ISakGTl ). the amount of CP violation and the 
strength of the electro-weak phase transition are insufficient to account for an 
asymmetry as large as ?7 ~ 6 x 10~^°. The usual theoretical attitude towards 
the cosmic lepton asymmetry rji is that sphaleron effects before/at electroweak 
symmetry breaking equilibrate the cosmic lepton and baryon asymmetries to 
within a factor of order unity (the relatio ns in the limit of ultra-relativistic 
SMPP particles can be found in flHarQOh l If this is the case, for all phe- 
nomenological purposes r/^ is vanishingly s mall. Sphalero n effects are a crucial 
ingr edient in most baryogenesis scenarios (jDin04j : lRio99l ) , including leptogen- 
esis ( BucOSi : Dav08). Yet, no experimental evidence for or against these effects 
exists, and (even barring the — phenomenologically viable — alternative that 
r]B and its leptonic counterpart r]L are simply "initial cosmological conditions") 
models have been envisioned where the lepto n asyrnmetry i s large, as for exara - 
ple via Affleck-Dine mechanism or Q-balls flCas99l : lMar99l : iMcDOOal : lKaw02f ). 
Since charge neutrality implies that the electron density matches the proton 
one, we do know that a large lepton asymmetry could only reside in the neu- 
trino sector. This asymmetry can be parameterized in terms of the chemical 
potentials of the different flavor species, /Xi^^, or better in terms of the degen- 
eracy parameter = n^jTy^ which is constant in absence of entropy releases. 
For neutrinos distributed as a FD with temperature Tj^^, the asymmetry in 
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each flavor is given by 



Without further input, the quantities ^£ are not determined within the Stan- 
dard Model, and should be constrained observationally. Over the years, BBN 
with a lepton asymrnetry has been studied by many authors and in differ- 
ent scenarios (|Wag67l: IfyS: lKan92l: lKoh96l: IWhiOol: lE spOOal: lEs pOObl: IEsdOI 



Han02l : IOri0ilsti02ilBar0'3JjCuo04lKne04ISer04bfe Pop OSbl i lSimoiah 



There are several effects of 7^ on BBN. The most important one (at least 
for relatively small ^e) is a shift of the beta equilibrium between protons and 
neutrons, which is however insensitive to and ^r- The leading flavor-blind 
effect amounts to a mere modiflcation of the radiation density entering the 
Hubble expansion rate equation by the amount 



eS 



E 



7 U/ 7 U 



4n 



(88) 



Moreover, for sufficiently large ^ the n e utrino decoupling temperature is 



higher than in the standard case (lFre83l : lKan92l ) . so that in principle one 



could get a non-standard Ty{T) evolution. A non-zero slight ly modifi es the 
partial neutrino reheating following the e^e~ annihilation, too flEspOObl l Both 



effects are however typically negligible for the ranges of ^ presently allowed. 
The greater sensitivity to ^ than ^^^^ made the constraints on the latter 
quantities looser, allowing on the other hand a richer phenomenology within 
a quasi-standard scenario. 

Again, the new knowledge on neutrino mixing parameters has rescued the sim- 
plicity of the standard cosmological scenario. A few years ago it was realized 
that the measured neutrino oscillation parameters imply that neutrinos reach 
approximate chemical equilibrium before the BBN epoch. This is due to the 
effects of the background medium on the evolution of the neutrino matrix den- 
sity. In presence of neutrino asymmetry the medium term in the Hamiltonian 
becomes 

H^ = ^mg{V,,V^,Vr)±V2GF{Q-Q) , (89) 



with the -|- sign for i/, the — sign for u. In particular neutrino self-interactions 
synch ronize the neutr i no oscil lations and drive all the potentials to the same 
value f Dol02a ; Won02 ; Aba02l ) . Assuming the standard value for N^q., from the 
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Yp range follows the bound = —0.008 ± O.Olsf"!. In Figure 13 we show the 
predictions for the primordial light-element abundances as a functi on of th e 
neutrino degeneracy parameter ^, taken to be equal for all flavors (jSerOSbl ). 
The gray band i s the 1 a predicted range, including both the uncertainty on 
n of Ref ■ (ISpeOSi ) and the nuclear reactions and uncertainties adopted in Ref. 
flSer04bh . 



The bound relaxes by a factor of 2-^3 (depending on other priors used) if 
addit ional degrees of freedom are present in the plasma, i.e. iVefj is allowed to 
vary (jBar03al : ICuo04l : ISimOSal ). At present, the BBN is by far the best cosmic 
"leptometer" available, and is virtually the only one sensitive to the sign of 
(^e- Most of its sensitivity derives however from the sensitivity of Yp to the 
variation of the weak n — p rates, so in order to improve these bounds the 
systematic error in the determination of primordial helium remains the major 
obstacle. Yet, even lacking further progress in this direction, in the near future 
the BBN role will be still important in combination with other cosmological 
observables to break degeneraci es amorig different pararn eters, as in the case 
of the PLANCK CMB mission flHam08l : 1Pop08al : IPopOSbl l 



6.4 Sterile Neutrinos and BBN 



Sterile neutrinos are, by definition, standard model gauge group singlet fermions. 
Their only coupling to SM particles arises via their mass and mixing param- 
eters with active neutrinos and, provided their mass is not too high and their 



mixing sufficiently small, they are long-lived particlea I. Here, we shall only 
refer to the case where there is only one additional neutrino mass eigenstate 
z/4, with mass 7714. Even in this case, the resulting 4x4 neutrino mass matrix 
U is described by 4 masses, 6 mixing angles and 3 CP-violating phases (and 
possibly, other 3 phases, not entering oscillations, if neutrinos are Majorana 
particles). Namely, one mass, three mixing angles and two more phases with 
respect to the 3x3 matrix for active Dirac neutrinos. Since not even the 
mass pattern of active neutrinos (or their complete mixing matrix) is known, 
it comes with no surprise that sterile neutrinos can manifest quite a rich 
phenomenology. This is especially true in cosmology, due to the relevance of 
medium effects. Basically, sterile neutrinos with a typical mixing element with 



F or som e regions in parameter space, the equalization may not be complete, 
see (|Pas08l l. Also, at le ast one w ay around the equalization of the chemical potentials 



has been proposed in (|Dol04bl ): an hypothetical neutrino-Majoron coupling of the 
order g ~ 10~^ can suppress neutrino flavor oscillations in the early universe, in 
contrast to the usual weak interaction case. 

For a typical mixing element with the active sector of the order sin 9s one 
expects a lifetime for decay into a neutrino and a photo n of th e order ~ 
2048 TT^ / {9 aG'p sin^ 26s 1^^4)7 "^4 being the sterile mass scale 
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Fig. 13. Light-element abundances as a function of the neutrino degeneracy pa- 
rameter. The top panel shows the primordial ^He mass fraction Yp, whereas the 
other panels show the ^H, '^He, and '^Li number fractions relative to hydrogen. The 
gray 1 a error bands include the u ncertainty of the WMAP determination of the 
baryon abundances of Ref. (ISpeOah as well as the uncertainties from the nuclear 
cross-sections of Ref.( Ser04bl ). Updated, from (SerOsS). 



62 



the active sector of the order sin 6s can be populated via incoherent scattering 
with a rate which, under "normal" conditions, writes 



sin^ 29s Ta , (90) 



Ta ~ G\T^ being the active neutrino scattering rate. However, in the presence 
of matter with a potential V , the effective mixing angle can be efficiently 
enhanced whenever the resonance condition 6m^^/2p ^ is fulfilled, giving 
rise to potentially large effects even when a small vacuum mixing would make 
the sterile neutrino undetectable in laboratory experiments. BBN is sensitive 
to sterile neutrinos through the following three effects: (i) the partial and 
total population of a sterile state induces A'^cfr > 3 and thus affects the Hubble 
expansion rate; (ii) if are produced only after the decoupling of the active 
neutrinos from the cosmological plasma, they lead in general to a depletion of 
Ue and Ue, thus affecting the weak n — p rates; (iii) the depletion can he u — u 
asymmetric, again affecting in particular, the weak rates. 



The inter est in the physics of s terile neutrinos in BBN has a long history 



see e. 



flKirSSi : iBargol : lEnqQOl : iKaiQOl : iBarQll : lEnq92al : ICli92l : IShi93l : lLis99 



Aba03al )). Due to the otherwise large parameter space, it has roughly fol- 



lowed the appeal that, from time to time, sterile neutrinos have had in ex- 
plaining anomalies in the neutrino phe nomenology. A n incomplete account in- 
cludes m4 ~17 keV in the b eta decay (jFra95l : IWie96l ). the KARMEN anomaly 
(m4 ~33.9 M e V) (lArm95n . and in the last few year s, the L SND anomaly 
(m4 ~ leV) flAguOll ). The recent Miniboone results flAguOTl ). although not 
ruling out the possibility of more exotic physics, strongly disfavor or rule out 
the simplest sterile neutrino models to explain the LSND signal. At the mo- 
ment, there is no clear theoretical or experimental argument suggesting the 
existence of (sufficiently light) sterile neutrinos , yet t heir rich physics contin- 
ues to attract a lot of interest (for a review, see (jCirOSi )). However, this implies 
that in absence of theoretical or experimental prejudice, one has to scan over 
a large parameter space. 



Most of the old literature referred to mixing between a sterile and an active 
state, neglecting mixing among active neutrinos. This allowed for several sim- 
plifications, but the results are clearly unphysical given the fact that we know 
that neutrinos do mix, and the mixing is large. A density matrix formalism 
is necessary, which differs from the one we introduced in Section 3.2 in sev- 
eral points. The vacuum Hamiltonian Hq includes now four mass eigenstates 
and a 4x4 mixing matrix. The refractive term Hi writes now in flavor basis 
(considering only diagonal elements) 



ffi = diag(V;,l^^,K,0) . 



(91) 
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The matter potentials Ve for each flavor is (with r]iy = rjy^ + vj^^ + vj^^ ) 



Vn 



Ve = ±V2GFn^ 

= ±V2GFn^ 

Vr = ±V2GFn^ 
where + applies to i^, — to 



Vn . 



8V2G 



FP I Pv^+Ve _|_ Pe-+e 



Ml 



8V2G 



fP 



SV^GfP 



3m: 



2 Pur+^T 



(92) 
(93) 
(94) 



As initial conditions one assumes usually thermal populations for the active 
neutrinos and a vanishing one for the sterile n eutrinos (this assumption is 
relaxed in some papers, as (jKir04l : iKirOTl : iKirOGi ). where a partial filling of the 
initial sterile state has been considered). In general, assuming that the active- 
sterile mixing angles are small (to be consistent with the phenomenology in 
the laboratory) is the only reasonable simplification. Also, as long as T -C m^, 
in the limit ^23 = t^/^ and there is a — r symmetry which further 

simplifie s the structure of U. A quite thorough analysis has been performed in 
(jDol04a] ). at least for the range Sml^ < leV^ which gives rise to the majority 
of phenomenologically distinct cases. We summarize here the main features, 
while addressing to the original literature for details. 

• If the terms in the square brackets are very small or have the natural value of 
the baryon asymmetry, rji ~ 10^^, they are dynamically negligible compared 
to the terms of order O^Gf/M^t z) at high temperatures and to the vacuum 
term Hq at low temperatures. This is also the situation considered for the 
standard decoupling in Section 3.2. In this case, the matter potential is 
always negative, and the existence of resonance conditions depends only on 
the mass-square differences 5m\. If 5m\ > 0, Vz, the sterile-active mixing is 
never resonant, and the analysis simplifies considerably. If however Sml^ < 
for some value of i, the system may undergo one, two or three resonances, 
and many sub-cases are possible. 

• If a large neutrino asymmetry is present, the square brackets terms may 
be large enough to change the impact of sterile n eutrino s on BBN, typ- 
ically weakening the constraints, as first noted in (jFoo95[ ). The reason is 
that the dominance of the flavor-diagonal medium term compared to the 
off-diagonal term due to mixing suppresses the active-sterile oscillations, 
producing sterile neutrinos less efficiently than in a symmetric background. 
Apart for larger value of the potential at a given temperature, another pe- 
culiarity of this case is that its sign is opposite for u and u: a pattern of 
resonances appears ind ependen tly of the sign of Smf, and each one only in 



the ly (or u) sector. In (jChuOGi ). some attention has been paid to strategies 
to lift the cosmological bound on sterile neutrinos invoked to explain the 
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excluded by exp 




sin 26lsnd 



Fig. 14. The allowed LSND region at 99% C.L. (yellow/light shaded area) com- 
pared to the cosmological bounds from BBN and LSS in the presence of primordial 
asymmetries. The darker shaded area is already excluded at 99% C.L. by other ex- 
periments. The regions below and to the left of the thin lines are allowed by BBN 
because they correspond to 1^ < 0.258. The regions below and to the left of the 
thick lin es are allowed by LSS because they correspond to Q^h'^ < 0.8 10~^. Prom 
(jChuOfil V 



LSND anomaly with a moderate asymmetry (say, rj ~ 10~^). The weak- 
ening of the BBN bounds for a growing asymmetry is represented by the 
shift from the sohd purple line to the dashed ones in Figure 14. Apart for 
the asymmetry, the distortion of the neu trino momentum distributions is 
negligible in the cases studied in (jChu06l ). It is well-known, however, that 
significant deviations from a pure FD distribution could occur during the 
evolution. Typically they c an be relevant either for relatively small mass 
splittings, 5m|j < 10~^eV^ (iKirQSl ) or, for eV scale masses, via matter res- 
onances post weak decoupling, which might leave both active and sterile 
neutrinos with a highly nonthermal spectrum for some choices of neutrino 
parameter s and energies. This however requires larger initial asymmetries, 
7]^ > 0.01 (lAba05l : ISmiOd ). 



An interesting aspect implicitly omitted above is that, even if the neutrino 
asymmetry is vanishingly small before the onset of the active-sterile oscilla- 
tions, a large u — u asymmetry (up to 77 ~ 0.1) might be dynamically generated 
in the active neutrino sector (and compensated by an opposite one in the sterile 
sector) via a resonant matter effect. This scenario requires quite special choices 
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of the parameters: the mixing must be sufficiently small, sin^ 26s ^ 10 
\5mf\ < f or some i, and \5m^ \ sin^2^s < few x 10 ~^eV^. For a review we 



address to (jDol02bl ) and reference therein (see also (|Dol04aJ) for some more 
details). Perhaps more general is the following consideration, which is often 
overlooked in the literature: once a sterile neutrino is introduced, at least two 
additional CP- violating phases are naturally introduced in the neutrino mixing 
matrix. The vacuum oscillation probability between an active and the sterile 
state are naturally CP-violating, even when ^13 —>■ and no CP violation hap- 
pens in the active-active oscillations. One might generate a flavor- dependent 
lepton- asymmetry as large as 0(sin^ 2^^) — which in a large part of the allowed 
parameter space is much bigger than t]b — without resonances. 

Finally, some BBN constraint also arises for more massive sterile neutrino 
states, in the keV-MeV range, by requiring that the energy density stored in 
the sterile states p opulated via mixing does not exceed e.g. AA^efr = 1- As re- 



viewed in ( lAbaOll ). the bounds are not very competitive compared with others. 
For significantly heavier states, the most interesting bounds may come from 
cascades and dissociation of light elements from the sterile neutrino decays, 
which we address in Section 9. 



7 Inhomogeneous nucleosynthesis 



In the standard scenario of nucleosynthesis all constituents are homogeneously 
and isotropically distributed, in accordance with the hypothesis of a FLRW 
universe. Anisotropic models, studied as early as in the 60's, or models with 
adiabatic fluctuations in the radiation energy density may affect nucleosyn- 
thesis essentially via a variation in the expansion rate. The emergence of a 
concordance cosmology model, supported by CMB and LSS data essentially 
confirm observationally that, apart from small initial adiabatic fluctuations 
(~ 10~^) as seed of structure formation, no signiflcant departure from the 
homogeneity or isotropy is required. The logical possibility that large fluc- 
tuations existed at the horizon scales at BBN epoch (size of the order of 
10~^ deg. in the CMB) is not well motivated either in the favored inflation- 
ary scenario to generate the perturbations. Thus, this li ne of re search has 
faded away in the last decade: we address the reader to (lMal93l ) for a his- 
torical overview. However, it is perfectly consistent with the present cosmo- 
logical scenario to speculate on a varying baryon to photon and neutron to 
proton ratio on small scales. These isothermal (or isocurvature) fluctuations 
may alter BBN in a non-trivial way, and these scenarios are known as in- 
homogeneous BBN (IBBN) models. Till the 90's, the interest in IBBN was 
due to two aspects: (i) several theories can lead to inhomogeneous distribu- 
tions of neutrons and protons at the time of nucleosynthesis: a flrst order 
quark-hadron phase transition ( WitSJ; lKur88al : lSuni90|; 



Maim , the CP vio- 
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lating interaction of particles with the bubble enucleated in the electroweak 



phase transition ( Ful94 : Hec95 : Kai99 : MegOSl ), the ph ase transition involv 



ing infl ation-g enerated isocurvature fluctuations (jDol93l ) or kaon condensation 
phase ( Nel90l ) : (ii) there was some hope that IBBN scenarios with — 1 — 
thus consistent with a flat, matter dominated universe, or at least an open uni- 
verse without non-baryonic dark matter — were phenomenologically viable (see 
e.g. ( Salsd )). Both motivations have lost observational or theoretical support 
in the last decade. The agreement between the value of Qsh"^ from standard 
BBN deuterium abundance and CMB anisotropies is indeed quite a compelling 
result in favor of the simplest homogeneous scenario. Theoretically, due to the 
high mass of the Higgs, the electroweak phase transition is a smooth cross-over 
for the SM particle content, and even the QCD phase transition — although less 
well established — appears to be a crossover or weak flrst order one, probably 
i nsuffic ient to produce significant departures from the standard BBN scenario 
fISchOSh . 



In the following, we limit ourselves to introduce a nd briefi y describe some 
recent calculations in IBBN, addressing the reader to (jMal93l ) for a throughout 
overview of the topic. Since IBBN has to reproduce very closely the SBBN 
yields while requiring additional parameters, most of its phenomenological 
interest has declined, too. Today, perhaps the only way to discriminate IBBN 
vs. homogeneous BBN lies in the different predictions for the intermediate 
(CNO) or heavy elements. 



7.1 Baryon inhomogeneous models 



An important point to be considered is the sensitivity of the predictions to the 
form of the inhomogeneity itself. Many studies simply have a 2-phas e mode l 
with a fraction / at high density and 1 — / at low, while in the paper (jBar83l ) 
a more realistic distribution (log normal) is considered, which allows a bound 
to be placed on the variance. 



If there are large fluctuations in the nucleon density, the differential transport 
of neutrons and protons can create neutron-rich regions where heavy elements 
can be formed. Neutrons diffuse by scattering on electrons and protons, pro- 
tons scatter on neutrons and Coulomb scatter on electrons, but the mean free 
path of protons is about 10^ times smaller than that of neutrons. The diffusion 
of other species is negligible with respect to neutron scattering due to their 
larger masses. Neutron diffus i on, how ever, was not considered in the earliest 
codes of IBBN ( 1Zel75l : lEps75l : lBar83l ) . where regions of different nucleon den- 
sity were treated as separate homogeneous BBN models. The mass fractions 
from each model were then averaged , with a weight given by the correspond- 
ing size. A later generation of codes (lApp87l : lAlc87l : lKaj90l ) introduced in the 
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calculation nucleon diffusion, but only at early times and high temperatures, 
before the starting of nucleosynthesis. This led to the neutron enrichment of 
the low-density region but, once the original protons were consumed, neutrons 
could form ^He only when other protons were produced by neutron decay. The 
main consequences of this situation on the light element abundances were: a) 
since four neutrons (two of which decaying in two protons) were needed for 
each "^He nucleus, the final yield of ^He was reduced; b) nucleosynthesis time 
scale were tuned by neutron decay rate, extending the process to cooler tem- 
peratures and allowing to survive in larger quantities; c) the high neutron 
density could help the production of heavier elements through neutron-rich 
isotopes. However, once neutron diffusion during nucleosynthesis is taken into 
account, all the three previous effects are weakened, since when neutrons are 
rapidly consumed in the high density region, where nucleosynthesis begins 
first, the excess neutrons in the low density region diffuse back. Kurki-Suonio et 
al. made the significant step forward of treating nucleon diffusion both before 
and duri ng nucleosynthesis, with planar sym metric baryon inh omogeneities 
( Kur88bl ) or cylindrical and spherical models (Kur89; KurQOa ). In order to 
decrease the nur nber of zones needed to obtain a high accuracy, nonuniform 
grids were used f lKurQObI : iMatQOl : lMat96l : lOriQTI : iLarOSl ). The diffusion equa- 
tion is sufficient for describing the motion of particles in an IBBN model if 
the background fluid is stationary, as in the case of neutrons, which are much 
more mobile than the ions and electrons they scatter on. The evolution of 
ions at low temperatures is more complicated, due to momentum transfers in 
collisions with other ion components, which move with comparable fluid veloc- 
ities. In this case, the common diffusion approximation has to be relaxed and 
one needs to tak e into account dissipative processes through hydrodynamic 
equations flAlcQOl : Ijed94al : lKei02l ). 



In an inhomogeneous code with treatment of neutron diffusion, the region 
considered is divided into several zones, s, where the time evolution of the 
number densit y of the specie z, rij obeys the following evolution (jLar05 : 
Kai99l : iMatQoh 



'~'"'i,s J, J- I -p ^ k,s l,s -p ^ i,s -' J 

nB,s i ki^ij i ij^ki 



The first three terms are usual, corresponding to reactions which create or de- 
stroy nuclides and to the expansion of the universe, while the last one is due to 
diffusion of isotope i between zones. The parameters which appear in Eq. (95) 
are the inhomogeneity distance scale, r, which measures the physical distance 
between inhomogeneity regions at the starting temperature, T ~ lOMeV, the 
stretching function, ^(r), which implements the non-uniform grid, marking the 
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zone edges, and the neutron diff usion co efficient, Dn, which is a function of 
proton density and temperature ( App87 ) P^ . For small distance scales, r < 1 



light-hour, the inhomogeneities are smeared out by neutron diffusion before 
nucleosynthesis starts, and the IBBN results approach standard BBN results. 
The constant parameter p changes with geometry (for example, for the spher- 
ical symmetry p = 2). Other important parameters are the density contrast, 
R, which is the ratio between the high and low densities, taken as high as 10^, 
and the volume factor, /„, that is the fraction of space occupied by the high 
density region. Note that the higher R, the larger the number of zones needed 
for a sufficient accuracy of the calculation. 

The plots in Figure 15 are co ntour m aps of the ^He mass fraction, and ^H/H 
and ^Li/H, taken from Ref. (ILarOSl ) . and present the characteristic features 



of an IBBN prediction on light element abundances. For small values of rj 
neutron diffusion homogenizes neutrons very quickly and protons as well, since 
this happens before weak interactions go out of equilibrium. This means that 
the final abundances are the same as a homogeneous model. The shift in 
the contour lines towards low values of t] for distance scales starting from 
~ 2 — 3 X 10^ cm is due to the fact that, for these values of r^, neutron diffusion 
and homogenization take a time of the order of the weak interaction freeze- 
out. This implies that protons start to be not efficiently homogenized and 
nucleosynthesis occurs before in the high density shells with a larger proton 
density and has the characteristics of an earlier nucleosynthesis (larger ^He 
and ^Li and less ^H). So, to recover the same values of the abundances of a 
model with lower rj one needs lower values of rj. Up to ~ 10^ cm the proton 
number density is unchanged except for a slight increase due to neutron decay, 
and this explains the almost vertical contour lines in this range. The depletion 
of neutrons in the high density shells at temperatures just after nucleon freeze- 
out leads to neutron back-diffusion and a more efficient nucleosynthesis (and 
a magnification of the production of ''Be) in the high density region. Due to 
this overproduction of ''Be, which gives ^Li after decay, the contour lines in 
the lowest plot in Figure 15 have a larger shift to lower rj. When starts to 
be larger than 10^ cm, neutron back-diffusion does not affect all the shells 
and nucleosynthesis is concentrated only in some of them. This leads to a 
decrease of the final '^He abundance accompanied by an increase in deuterium 
production, which corresponds to a shift of the contour lines towards high rj. 
Finally, contours turn back to low values of t], since for very large diffusion 
cannot homogenize neutrons before nucleosynthesis, a large neutron density 
remains in the high density region, giving rise again to ^He overproduction 
and a suppressed yield. In this region of r,, results are equivalent to the 
average of two separate homogeneous BBN models, one of high density (with 
high ^He and ^Li and low ^H) and one of low density. The basic shapes of 



■^Note that, as remarked in (lKur90bh . in their Eq. (21) the factor vr/lG IS missms; 
from the numerical value. 
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Fig. 16. Concordance betwee n the observ ational constraints on ^He, ^H, and ^Li 
and the model of Ref. ( LarOSi ) . Prom Ref. ( LarOSi ) . 



the contour lines of Figure 15 are common to all IBBN models: for different 
geometries and values of the parameters there will be regions in where 
neutron homogenization and diffusion occur at times between weak freeze-out 
and nucleosynthesis or after nucleosynthesis. 



Figure 16 shows the concordance betwee n the o bservational constraints on 
^He, ^H, and ^Li and the m odel of Ref. (ILarOSl ): upper plot is for the ^Li 
constraints from Rya n et al. while lower plot is for the ^Li data of 

Mele ndez fc Ramirez (jMslQJLhi both plots the concordance region between 
'^He (jlzoOJ) and data (IKirOSl ) is shown in yellow. While upper plot have 
a concordance region for ^Li only for a depletion factor ranging from 2.8 to 
5.9, the lower plot does not need any depletion in ''Li if rj < 5 x 1 0'^ (see Ref. 
(LarQGa) for a similar analysis using r„ = 878.5 ± 0.7stat ± O.Sgyst f Ser05al )). 
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One interesting consequence of IBBN results is the larger range in the deple- 
tion factor one can obtain for ^Li with respect to the analogous prediction of 
homogeneous BBN. Figure 16 shows, moreover, that IBBN allows for larger 
values of t] than SBBN (up to 77 ~ 10~^), requiring at the same time a depletion 
factor for ^Li to obtain concordance with the observational limits. 

Since lithium is produced quite late in nucleosynthesis, its yield is partic- 
ularly sensitive to the late-time transport phenomena such as hydrodynamic 
ion diffusion. In thi s respect, the results of codes wh ich tak e into account these 



phenomena (jKei02l ). which are not considered in (ILarOSl ) . might explain the 
observed depletion of lithium. This is a consequence of late separation of ele- 
ments due to Thomson drag: Thomson scattering of electrons on background 
photons makes the diffusion of ions inefficient which must drag electrons with 
them to keep charge neutrality. On the other hand, protons and helium ions, 
for instance, are allowed to diffuse in opposite directions. While protons diffuse 
out, helium and lithium get concentrated in the high density regions, leading 
to enhanced destruction of ^Li, ^H, and ^He. 

Heav y elemen t productio n in the framewo rk of IBBN was first investigated in 



Ref. f|Jed94bl ) and then in (IMatOSbl : iMatOTI ) . in the approximation of neglecting 



baryon diffusion. The authors claim that there is a parameter region, for the 
volume fraction / and density contrast R, in which heavy elements can be 
produced enough to affect the observation, while keeping the light element 
abundances consistent with observation. The results show that BBN proceeds 
through both the p-process and the r-process, with the transition between the 
two due to the Coulomb barriers of proton-rich nuclei and the amounts of 
neutrons when heavy elements begin to be synthesized. 



7.2 Matter-antimatter inhomogeneities 



A possible scenario which gives rise to an inhomogeneous baryon- t o-photon 



ratio is antimatter BBN (ABBN) flSte76l : iRehQSl : iRehOll : iKurOOal : iKurOOb 



SihO if ). Different baryogenes i s models can giv e rise to matter -antimatter do- 
mains fiGioQSl : iDoiga IDoIQSI : iKhlOOal : iKhlOObI : iMatOi IPolOSh . In the ABBN 
scenario, the antimatter regions have radius r^, while R is the antimatter- 
matter ratio in the universe. Antimatter regions should be small enough to 
be completely annihilated well before recombination, in order to satisfy CMB 
constraints. Their size determines the time when most of the annihilation takes 
place, before or after significant amounts of ^He are produced by nucleosynthe- 
sis. The first case is reali zed for typ ical radii between 10^ and 10'' m (comoving 



distance at T=l KeV) flGio02l lP. Thanks to the different diffusion scale of 



Smaller antimatter regions would annihilate before neutrino decoupling without 
any effect on BBN. 
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Fig. 17. Light element yields in ABBN as a function of r] and iVcg for ta = 10^'^ 
m and, from top to bottom, R = 10~^, 10~^'^, 10~^'^. Thick sohd lines, thin solid 
lines, a nd dashed lines are for Y^, Zo^i^qH/H, and Zof^^QLi/H, respectively. From Ref. 
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Fig. 18. Combin ed allow ed region in (7?, Np ff). The dashed line is the CMB+SNIa 
constraint from ( HanOd ) . From Ref. ( Gio02l ) . 



protons and neutrons, the latter can more easily move to antimatter regions 
and annihilate. This would produce a reduced neutron to proton ratio with 
respect to the standard case, which can be compensated by a larger expansion 
rate at BBN, provided by more relativistic degrees of freedom, Ngg > 3, and 
results in the same Yp. Correspondingly, the speed-up of expansion shortens 
the time interval available for nucleosynthesis, and implies smaller yields for 
all other light nuclides, yet this can be compensated by the increase of the 
reaction rates due to a higher value of 77. The net result is thus, a shift of the 
agreement of theory vs. data towards larger 77 for large AA'^efr, see Figure 13, 
but with the bonus that N^s is no longer constrained, as shown in Figure 18. 



8 Constraints on fundamental interactions 



8.1 Extra Dimensions and BBN 



^.1.1 A short journey to extra dimensions 



The idea of introducing extra (spatial) dimensions to generalize the 4-dimensio- 
nal theory of fundamental interactions and unify different forces is quite old. 
As early as 1919, thus shortly after the birth of General Relativity, Kaluza 
considered a 5-dimensional version of Ein stein th eory which described gravity 
and electromagnetism in a unique setting (jKal2ll ). Shortly after, in 1926 Oskar 
Klein stressed the role of having a compact fifth dimension in orde r to ev ade 
constraints from observations of large a ccessibl e extra dimensions ( lKle26l ). A 
nice review on Kaluza-Klein theories is (IBaiSTl l. 
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Higher-dimensional theories had perhaps their golden age starting from the 
late 70 's, after the discovery of the remarkable properties which superstring 
and supergravity theories have for particular spacetime dimensionalities. Quite 
recently, theories with one or more extra dimensions with a fundamental scale 
of TeV~^ have been advocated as possible way to address the long-standing 
proble m of h i erarchy between the elec troweak and the much higher Planck 



scale ( lArk98l : lAntQSl : lRan99al : iRanQQbl ). In these scenarios, the fundamental 



gravity scale is lowered down to the TeV range, and the observed Planck 
mass emerges as an effective scale at low energies, smaller than the Kaluza- 
Klein (KK) excitation mass scale. This is due to the dilution of gravitational 
interactions in the large (millimeter-sized) extra dimensions (flat scenarios), 
or the particular configuration of the gravitational field which provides a static 
solution to Einstein's equations (warped extra dimensions). A general feature 
of these theories is to assume that ordinary matter is confined to standard 3-1-1- 
dimensional spacetime, a brane embedded in a (4-|-d)-dimensional manifold, 
while gravity can propagate in the whole higher dimensional spacetime. 

Interestingly, large (experimentally accessible) extra dimension models can 
be t ested using co llider physics, as for example at LHC, for a review see 
e.g. flGiu07l : ISunOSh . On the other hand, they may have a large impact on 



the cosmological evolution of our universe. The issue of understanding the 
phenomenological implications of "brane-cosmology" has been addressed by 
several scholars in the last ten years, mainly aimed at discussing how these 
scenarios can be constrained by cosmological observables, CMB and BBN 
among others. In the following, after a brief summary of the aspects of the 
extra dimension models which are relevant for our discussion, we (mainly) 
focus on the constraints which can be obtained by exploiting BBN. 

We start by introducing the D = A+d dimensional Einstein and matter action, 
which can be written as 



S = I d'xd^yV^-^R+y^L^ , (96) 



where the first term corresponds to the Einstein-Hilbert action, R being the 
4 + (i-dimensional scalar curvature for the 4 + d metric g and Md playing the 
role of the D-dimensional reduced Planck mass, while the second term contains 
the matter Lagrangian, with the SMPP fields localized on the 3-1-1- dimensional 
brane y = 0. In the case of compact extra dimensions (we will consider the 
specific case of a (i-dimensional torus of radius 6) and for a factorized metric, 
i.e. if the 4-dimensional part does not depend upon the d extra coordinates, 
the action can be reduced to a 4-dimensional action at low energy, smaller than 
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the inverse compactification radius S, by integrating over the y coordinates 



S 



-R 



-gLr, 



(97) 



from which we read the expression of the Planck mass Mp = 
GeV in terms of Md and 6, 



1/2 



1.2-10 



19 



(27r5)-^ = Md{V^Md/Mp] 



2/d 



(9J 



For Mp, ~ TeV, the simplest case of one extra factorized dimension is excluded 
as it leads to a value of 5 which is too large, of the order of the scale of the 
solar system, while the scenario is viable for > 2, as in this case 5 < 1 mm. 

More generally, one can allow for an expli cit depe n dence o f the 4-dimensional 
metric on the extra coordinates. As in f BinOOal : ChuOO ). and usually con- 
sidered in almost the whole literature, we consider a (non-factorized) d = 1 
model, where the extra dimension (the bulk) is compactified on the line seg- 
ment S^/Z2- The metric (preserving 3-dimensional rotation and translation 
invariance) can then be written as 



-n 



(r, y)dT'^ + a^(r, y)dx^ + 6^(r, y)dy'^ 



(99) 



The Z2 symmetry identifies the points y and — so one can restrict to < y < 
1/2. Two three-branes are placed at y = (our "visible" brane) and ?/ = 1/2 (a 
hidden brane, which absorbs the gravitational flux lines of the visible brane). 
The metric is obtained as usual from Einstein's equations (upper case latin 
indexes A, i? = 0, 1, 2, 3, 5 run over the 5-dimensional spacetime), 

Gab = M-^Tab , (100) 



where the stress-energy tensor is the sum of contributions of ordinary matter 
on the visible brane, bulk matter and fields living on the hidden brane, 

rpAB rpAB 



Each term corresponds to a perfect fluid, parameterized as usual in terms 
of the energy density p and pressure P and specified by the equation of 
state P = P{p), with furthermore T°j^ = T°f^ = 0, so that there is no 
flow of matter on the branes along the fifth dimension, and finally T^^fj^ = 

diag{—pbulk, Pbulk, Pbulk, Pbulk, Pbulk,b)- 
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Before considering the cosmological scenarios corresponding to this framework, 
one should look for the static solutions of Eq. (100), the analogous of (empty 
space) Minkowski spacetime. Apart for the case T^^ = which leads to a 
(factorizable) trivial spacetime, Randall and Sundrum (RS) ( Ran99a : Ran99bl ) 



found a new solution by considering a pure bulk and brane cosmological con- 
stant terms. Choosing T^"^^ = Adiag{—1, 1, 1, 1, 1), a static solution is in fact 
obtained if A < and the two brane tensions PvisMd = ^visMd fine-tuned 
to the values 



Kis = -AMd = ±^-QAMl , (102) 



so that the effective cosmologi cal const ant in the 3-dimensional space exactly 
cancels. In this case one finds flRan99ah h{T,y) = bo = const and 



a{T,y) = exp ±bo\y\ 



EE] 

6M| J 



exp {±bQ\y\m) 



(103) 



Choosing a negative value for Ayis, so that the solution corresponds to Eq. 
(103) taken with the positive sign, leads to a nice solution of the hierarchy 
problem, since the fundamental mass scale on the invisible brane M5 ~ Mp is 
red — shifted on our visible universe by the conformal factor exp {—bom/2), 
which can explain the large relative ratio of Planck and electroweak scales for 
a moderate value of mbo ~ 10^. In other words, the non-trivial dependence of 
the metric upon y implies that the KK zero-mode of the graviton wavefunc- 
tion is peaked around the invisible brane and has an overlap with the visible 
brane suppressed by the exponential "warp" factor exp(— m6o/2). Yet the KK 
tower mass gap is potentially as low as the TeV scale, and thus these graviton 
excitations can lead to testable effects at high energy colliders as LHC, see 



e.g. (jGiuOTI ) and references therein. 



On the other hand RS also observed that choosing our brane with a positive 
cosmological constant Ayis > 0, though does not solve the hierarchy problem, 
nevertheless it has the nice properties of allowing for a non compact fifth- 
dimension, as one can take the limit y — > 00 maintaining consistency with 
short- distance force experiments. This scenario is also the one which is more 
interesting from the cosmological point of view, as we will discuss soon. 



8.1.2 Brane cosmology and BBN 

Adding matter on the branes will lead to an evolving universe analogous to 
the standard FLRW model. Depending on the choice of the corresponding 
reference static solution one starts with, the prediction for the Friedmann-like 
equation governing the visible scale factor, i.e. the value of a in the vicinity of 
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our brane y = can be significantly different, leading to testable predictions 
for cosmological observables such as BBN, CMB and structure formation. 



The e quation governing the evolution of a(r, 0) = ^oiT) has been worked 



out in (jBinOOal ). If p and P denote energy density and pressure on our visible 
universe, thus dropping the index vis in the following, one obtains the standard 
conservation equation, 

p + 3(p + P)— = , (104) 
ao 



which leads to the usual power behavior for p ~ ao as well as the 



evolution equation for ao, 
ao , f ao\^ 1 



p{p + 3P)-- 



ao Vao/ SeMf^'"^ ' Sb^Ml 

where the time derivative is with respect to t, with dt = n{T, 0)dT and a flat 
metric in the ordinary 3-dimensional space has been assumed for simplicity. 
This expression shows two remarkable properties, namely that it is indepen- 
dent of the energy density and pressure of the second brane, a manifestation of 
the local nature of Einstein theory and, furthermore, that the energy density 
enters quadratically rather than linearly as in conventional cosmology. 

If the dynamics is dominated by the brane energy density, so that one can 
neglect the last term in Eq. (105), using Eq. (104) the second order equation 
(105) can be put in the form 

±^ay) = ^_±(^p^a^) (106) 



which gives 

= —^p^ + 4 • (107) 
36Mr at ^ ' 



The second term in the r.h.s of this expression dep ends upo n the free integra- 



tion constant C and behaves as a radiation term (IBinOObl ). thus its popular 



name of "dark radiation", though the sign of C can be also negative. 

This result for the Hubble parameter strongly differs from the usual Friedmann 
law, unless one consider the very special case of a radiation dominated phase 
driven by a positive C, while for a negligible value of C one gets 
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thus a slower expansion rate compared to the standard result ao{t) ~ ^2/(3+3«))_ 
If we assume that Eq. (107) can be applied to the present universe, writing 
the energy density as a fraction of the critical density ~ 1, 

^2 H^Mj, C 6^ C , , 



with the radius of the fifth dimension, 5 (see Eq. (98)), which thus should 
be of the order of the present Hubble radius, Hq^ QfJ. This is ruled out by 
observations (since we could then observe five- dimensional gravity directly), as 
well as the possibility that dark radiation provides the dominant contribution 
to the expansion today. 

Strong bounds on this model also come if we assume that it describes the 
evolution of the universe at the earliest stage we can probe in a quantita- 
tive manner, namely during BBN. A rough constraint can be obtained by 
considering the different behavior of the Hubble expansion parameter, which 
changes the neutron to proton ratio freezing temperature, Td, as well as the 
time-temperature relationship in the ter nperatu r e range from To down to the 
deuterium formation at ~ 0.1 MeV ( BinOOal : ChuOd ). The value of n/p at 



Td is given by the standard relation n/p= exp(— Am/To), withP^ 

Gin^j^ . (110) 



On the other hand, the time-temperature relationship during a radiation dom- 
inated epoch is given by Eq. (108) with w = 1/3, t ~ T~^. Inserting numerical 
values and using the expression of the relativistic degrees of freedom, g^, during 
BBN one gets 

TD~7.5(^^j MeV , (111) 



and 



t ~ 2.8 



4 / M5 



/ T 



VTeV/ VMeV 



(112) 



Thus, assuming that all neutrons at T/v are eventually burned into ^He nuclei. 

Derivation of this result is presented in a slightly different form in (jBinnOal ^. where 
the second order equation for the scale factor is used. 

We consider the case C = 0. If dark radiation dominates at BBN, it would be 
difficult to reconcile later evolution with, say, CMB and structure formation data. 
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we have 



n 



2-K 



exp 



Am 



7.5MeV VTeV 



t(Tj^)-t{Tjy) 



(113) 



which imphes that a correct mass fraction Yp 
a too large compactification radius. 



0.25 requires M5 ~ 8 TeV and 



It should be noted that considering more than one (flat) extra dimension and 
an empty bulk could be potentially in better agreement with both BBN and 
the requirement of sub-mm extra dimensions, but a careful analysis of this 
scenario has not been considered in the literature in details, evaluating the 
whole network of light nuclei produced during BBN, in particular and ^Li. 
We also mention that it has been pointed out that in general, the presence of 
the KK tower of gravitons may lead to overclosure of the universe, unless the 
highest temperature ever achieved was of the order of MeV, thus with a severe 
impact on the whole BBN scenario, as well as on the stan dard inflationary 
picture for early production of perturbations (IHanOll : iFaiOll ). 



Interestingly, a much more promising brane cosmology can be obtained by 
exploiting the RS model. Assuming that the stress-energy tensor T^^^ cor- 
responds t o a cosm ological constant, pbuik '- 



bulk 



bulk,5; 



it has been 

shown in (jBinOObI ) that one can integrate the (0,0) component of Einstein's 
equations and obtain the generalized Friedmann equation in the vicinity of 
the visible brane, 



1 

^Pbulk + 



36MI 



Pvis ~^ 



C_ 



(114) 



This result holds independently of the metric outside and in particular of the 
time evolution of the scale factor b. If one assumes that the energy density pms 
can be decomposed as the sum of the contribution of ordinary matter p and 
a (positive) cosmological constant A^s, and the latter is flne-tuned as in the 



RS model, see Eg. (102), one recovers a standard cosmology (jCli99l : ICsa99 
ChuOd : iBinOObh . 



A. 



18M^^ ' SGMf 



1 



P 



C_ 



(115) 



if one identifles 

A,: 



87rG 



N 



6M| 



(116) 
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and the limit A^j^ >> p is assumej^. Of course, choosing the original RS 
proposal, which provides a beautiful solution to the hierarchy problem but 
requires a negative tension on the visible brane, one would obtain a negative 
sign relative to the conventional Friedmann equation, so that the visible brane 
behaves as an anti-gravity world (jShiOOl ). In particular, this implies that, as 
soon as the universe becomes matter dominated, i t would collapse on a scale 
of the order of the matter-radiation equality time (ICsa99l ). 



From Eq. (115) we see that adjusting the value of A^j^ as in Eq. (116) for each 
given M5, the Hubble rate depends on two new parameters, the fundamen- 
tal scale M5 which controls the quadratic term in p and the dark radiation 
constant C. Bounds on both these parameters can be obtained from BBN, 
and have been discussed by sev eral authors u s ing the simple argument on_^He 



mass fraction described above flCsa99l : ICli99l : iBinOObI : iMaaOOl : iLanOll : iMizOl 



Bar02al : iFlaOGi ). namely requiring that both terms be sufficiently small that 



an acceptable value for Yp be produced. A more careful analysis based on a 
full numerical integration of the BBN dynamics and conside ring th e predic - 
tions for a nd ""Li as well, has been instead performed in (jlch02al : lBra02l ). 
In particular (jlch02al ) only consider the dark radiation term and its effect on 
both BBN and CMB, neglecting the effect of linear perturbations of dark ra - 
diation during photon decoupling and recombination. The result of fllch02al ) 
(and of (jBra02l ) in the limit of large M5) can be easily translated in terms of 
the well-known bound on the effective number of neutrino species (see Section 
3), 



C 



Stt 



;il7) 



The allowed range for C can be easily obtained from the result on A'efr of 
Section 3. More interestingly, if the value of M5 is sufficiently low, the effect 
of the term can be non- negligible, but it can be compensated by a large 
and negative value of C. A degeneracy is thus expected in the M5 — C plane 
which qualitatively is of the form a/M^ + AA'^cfr = const, see Figure 19. 

We have assumed in the previous discussion that the size of the extra di- 
mensional sp ace is s tabilized by some dynamical mechanism, as for example 



discussed in ( lGol99l ). but indeed it might have some dynamics during cosmo- 
logical epochs. In general, for a homogeneous and isotropic model, the evolu- 
tion of the extra dimensions is controlled by a single scalar field, the radion, 
with a canonical kinetic term which interacts with ordinary matter, see e.g. 



( iLiOd ) . These interaction terms lead to a time dependent behavior of the Higgs 



vacuum expectation value, v, which in turn affects fermion masses. The BBN 



The standard behavior is indeed recovered only in the stron ger hmit p « 
A^j^/Mp if the bulk spacetime is not exactly anti-de Sitter ( ShiOOl ) . 
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Fig. 19. The BBN- allowed range in the M5 — C plane. Prom (|Bra02l ) 



is extremely sensitive to such variations, since changing v produces a different 
Fermi coupUng constant, Gp, and neutron proton mass difference, both enter- 
ing the determination of the freezing temperature, Tp, as well as shifting the 
pion mass which influences the nucleon potential and the deuterium binding 
energy. We will discuss this issue in details in the next Section. 



Finally, BBN can be also used to put constraints on theories which consider 
the possibility of bulk neutrinos. As we have mentioned, in the extra dimension 
scenarios the SMPP particles are assumed to be localized on the visible brane, 
while gravitational interactions can propagate in the bulk. If the fundamental 
scale is of the order of TeV, this leads to a serious problem in understanding 
the smallness of neutrino masses, which is usually thought to be produced 
via a see-saw mechanism. For sub-eV neutrino masses, this scheme requires 
the existence of some new mass scale of order 10^^ — 10^^ GeV, much larger 
than the extra dimension scale. Furthermore, operators as LHLH/Md {H 
and L are the Higgs and left-handed fermion doublets, respectively) could be 
induced in the low energy Lagrangian which lead to an unacceptable neutrino 
mass. Thi s problem was realized quite early on and severa l solutions have been 
proposed flPieQQl : lMoh99l : iMohOOl : lArkOll : ICalOlal : ICalOlbl ) , which postulate the 
existence of one or more gauge singlet neutrinos in the bulk which couple to 
the lepton doublet on the brane. Their corresponding KK modes give rise to 
an infinite tower of sterile neutrinos labeled by an integer n, which mix with 
active neutrinos and have masses typically of the order of m„ = n 10~^ eV 
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for an extra dimension size of mm. This mixing has a relevant effect on both 
solar and atmospheric neutrino phenomenology since for each mode with mass 
larger than the three active neutrino with Dirac mass /x, (i = 1, 2, 3) the re will 
be a corresponding vacuum mixing angle 9 ~ fii/rrin (see e.g. (ILamOll ) for a 
review and references therein). 



Sizeable effects are also expected in cosmology ( BarOd : Aba03b; Goh02 ). In 
particular, KK modes of bulk neutrinos can be produced in the early universe 
before BBN by incoherent scatterings or coherent oscillations, thus contribut- 
ing to the total radiation content parameterized by A'^cff. Bounds on this pa- 
rameter from BBN amounts to require that the whole tower of modes should 
be equivalent to no more than approximately one active neutrino species. Fur- 
thermore, the photoproduction of and ^Li by decays of modes after BBN 
may potentially spoil the whole standard nuclei production scenario. These 
effects limit the possible values for the extra dimension length scal e S as fun c- 
tion of the Dirac active neutrino mass, see e.g. Fig ure 20 (f rom flAba03bh l 
For illustration, for a neutrino mass of order 0.1 eV (jAba03bl ) find from BBN 
O.OlMeV < 6-^ < lO^MeV. 



8.2 Variation of fundamental constants 



8.2.1 Introductory remarks 

Already in 1937 P.A.M. Dirac first introduced the idea that the fundamen- 
tal constants of physics may be in deed y ariable parameters characterizing 



the particular state of the universe (jDir37l ). Physicists have long scrutinized 



this possibility. On one hand, a strong effort has been devoted to embed this 
paradigm into a definite theoretical framework. For example, theories with ex- 
tra dimensions such as Kaluza-Klein or string theories, naturally predict that 
4-dimensional constants may vary (in time and space), since they represent 
effective values in the low energy limit, and are sensitive to the size and struc- 
tures of extra dimensions. Any variation of these invisible dimensions, for ex- 
ample over cosmological times, would lead to varying 4-dimensional constants. 
On the other hand, measuring variations of fundamental constants has been 
pursued at the experimental level by several groups and techniques, ranging 
from short time laboratory based measurements, to astronomical or geological 
scale studies and, finally, to cosmological time variation searches. All these in- 
vestigations provided quite strong constraints on the possible time evolution 
of e.g. the fine structure couphng, a, or the Newton gravitational constant, 
G^. Both theoretical and experimenta l aspect s of this intriguing research issue 



are beautifully covered in the review (IUza03l ) 



In the following we will summarize the impact that fundamental constant time 
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Fig. 20. Constraints on the extra dimension length s cale R versus active neutrino 
masses fii for bulk neutrino scenarios. Prom (jAbaOSbl l. 



variations have on BBN, by describing the main effects on the hght nuclei 
abundances and the bounds which therefore is possible to obtain using BBN 
as a probe. Unless otherwise stated, we will limit to constrain the departure of 
these parameters at the BBN epoch compared to present values, independent 
of a specific theor etical model for their ev olution (this is instead the approach 
considered e.g. in (jBar02bl : ISan02l : IShaOSi ) . ) As a general remark, we would like 
to stress that all results depend quite strongly upon the general assumption 
which are made to perform the analysis, the priors in presently fashionable 
Bayesian language. In fact. 



1) there are too many constants which enter the BBN physics (the fine struc- 
ture coupling, the Newton constant, the strong interaction scale, Aqcd, the 
Yukawa couplings, the Higgs vacuum expectation value), so that if they are 
all considered as free independent parameters to be fixed by data, one lacks 
predictive power; 

2) there is no unique theoretical framework which allows for an unambiguous 
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determination of their relative evolutionary history. 



For these reasons, there are two typical strategies which have been exploited. 
One may assume that only a single fundamental constant (or a subset of them) 
is promoted to the role of a free parameter to be constrained, keeping all the 
others as fixed. Several analyses of variation of the fine structure coupling, 
which we describe in the next Section are based on this assumption. On the 
other hand one can consider a specific theoretical model, which reduces the 
number of independent parameter one starts with, and this typically allows 
for tighter constraints. Examples of this approach has been considered in the 
last decade in details, based on dilaton inspired theories, or rather on unified 
gauge theories, where all gauge coupling of the SMPP get unified at some high 
mass scale Mqut- 

To conclude this short introduction, let us stress that checking for a time (and 
space) variation of fundamental constants is only meaningful for adimensional 
quantities, such as a. This is due to the fact that measurement of dimen- 
sional parameters is strongly intertwined with both the system of units and 
the particular measurement technique which is employed, so that an absolute 
determination of, say, the time evolution of the speed of light, c, is meaningless. 

This can be illustrated with a simple example. Suppose we want to measure 
the value of c by using a light clock device, a source S and a mirror placed at 
a distance D away from S, which reflects the light ray back to S. The value of 
c is then computed as the ratio of the distance 2D over the total time elapsed 
from emission to light collection, which is expressed in terms of an adimen- 
sional number, once units of length L and time T are specified. For example, 
we can choose an atomic clock to specify T, using the hyperfine splitting tran- 
sition rate of caesium-133 atoms. In this case the particular combination of 
fundamental parameters mlc^a'^/mph (the typical hyperfine frequency) is kept 
fixed by definition, rrip being the proton mass. If we choose the length unit L 
as the distance traveled by light in k units of time T, it is rather trivial that 
there is no possibility to detect the time variation of the speed of light, as also 
c in this case is kept fixed by definition. On the other hand we may use the 
standard prototype platinum-iridium bar as the value of L, which depends on 
the interatomic spacing of the material, and thus ultimately on the value of 
the Bohr radius = h/nieCa. In this case if we find that at different epochs 
the value of c has changed in these units, this amounts to say that the ratio 



is time dependent. Therefore, either a or the adimensional ratio rrip/mf, or 
both, change with time. 



c rrip 1 




L/T rrie a 
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As this simple example shows, evidences for time varying fundamental con- 
stants are in all cases evidences for particular combination of adimensional 
quantities, which depends upon the particular choice of units which is adopted. 
In the following, when time changes of dimensional parameters are considered, 
a ratio of two independent and homogenous constants will be always implicitly 
understood, as for example the ratio Agc/j/mg, with rriq some quark mass, or 



8.2.2 Varying the fine structure constant 

There are several direct measurements in the laboratory on the varia tion of o; 
over relatively short time period, using different techniques, see e.g. ( 1Uza03l ). 
The geological limit from the O klo natural rea c tor is about < 10~^ 

over a period of few billion years flDamQa : iFujOOl : IOli02l ). Astrophysical obser- 
vations of high red-shift quasar absorption lines provid e the only evidence for 
a possible time variation of a flWeb99l : iMurOSi : lMur04j ) . 



5a/ a 



-0.57 ±0.11) ■ 10" 



(119) 



A result compatible with zero is inste ad reported by (jSriOJ; ICha04f l. which 
however has been criticized in (jMurOTl ). 



Over longer time scales, the value of a can be constrained by cosmological 
observables, such as CMB and BBN. Indeed, CMB anisotropies are a very 
good probe of a since its value is imprinted in the ionisation history of the 
universe. The main effects are a change in the redshift of recombination due to 
a change of hydrogen energy levels, the modification of the Thomson scattering 
cross-section whi ch is proportio n al to o; ^, and at sub leading level a change of 



He abundance (lAveOd : lAveOll : lMar02l ). If the value of a is increased, the 



last scattering surface moves towards larger redshifts, which correspond to a 
shift of first Doppler peak towards larger /'s. Moreover, this shift produces a 
larger Integrated Sachs- Wolfe effect (i.e. more power around the first peak), 
while the high multipole diffusion damping is decreased by a larger a, thus 
increasing the power on very small scales. With present CMB data, including 
those from the first year release from WMAP Collaboration, one finds the 
bound —0.06 < 5a /a < 0.02 at 95% C.L., while a Fisher matrix analysis 
shows that future experiments such as PLANCK should be able t o const rain 
variations of a during CMB formation with an accuracy of 0.3% (jRocOJ). 



The value of a enters the physics of primordial nucleosynthesis in several ways 
and, remarkably, BBN represents the earliest reliable probe of possible varia- 
tion of the fine structure coupling over cosmological times, though it suffers 
from being mod el dependent with resp ect to CMB analysis. This issue was 
first studied in (lKol86l : iBarSTI : ICam95l ). where the focus was on the abun- 
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danc e of '^H e, while a detailed analysis has been presented mainly in (IBerQQl ) 
and (jNol02l ). which we follow for our discussion. 



During the early stages of BBN, at the n/p ratio freeze out temperature, T^, 
the fine structure coupling affects the weak n-p rates in two ways. First of 
all, it changes the neutron-proton mass differe nce Am, which can be only 
phenomenologically parameterized as in (jGas82l ). where the authors find 



Am(MeV) ~ 2.05 - 0.76 [l + ^ 



(120) 



This result is obtained by studying the the behavior of the nucleon masses 
versus a, which determines the electromagnetic quark masses and binding 
energy. The fact that this parametrization, though quite reasonable, is not 
deduced from first principles in a QCD-based calculation, is the main source 
of possible uncertainties and renders all predictions model-dependent. Fur- 
thermore, weak rates also depend upon a when QED radiative and thermal 
corrections are included. However, these corrections are at the level of few 



percent, see e.g. ( lEsp99l ). so the effect is very small for moderate variation of 
a, i.e. 8a I a << 1. Since the ^He mass fraction is mainly sensitive to the n/p 
ratio at freeze out, while it depends more weakly on the whole set of nuclear 
reaction rates, this implies that at first approximation the whole dependence 
of Yp upon a is through Am, which also fixes the decoupling temperature . 
since 

2 

^ 1 + exp(Am/Tz)(Am)) ' ^^"^^"^ 



The dependence of other nuclei on a is more involved, and one should scruti- 
nize the way the fine structure constant appears in the set of nuclear reaction 
rates. The leading effect is due to the change of Coulomb barrier. Charged 
particle reactions at low energies take place via tunneling through a repul- 
sive potential, and this produces most of the energy dependence in the cross- 
section. 



a{E) = 




(122) 



where rj = aZiZ2y fi/2E, with Zi and Z2 the charges of the incoming nuclei 
and /i the reduced mass. If one neglects th e depen dence on a of the astrophys- 
ical S'-factor and of reduced mass, as in (jBer99l ). then only the Sommerfeld 



^® Notice that the Fermi coupling constant Gp does not depend upon the values of 
gauge couplings in the SMPP, but on the Higgs vacuum expectation value only. 
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parameter rj would change linearly with a, and o ne can o btain the thermally 
averaged rates versus 6a/ a, see Tables 1 and 2 of (lBer99l ). as well as the frac- 
tional variation of light nuclei abundances. Not surprisingly, the most dramatic 
changes are those of ^Li, due to the strong Coulomb barrier in its production. 



The analysis of (jBer99l ) has been improved in (jNol02l ) . see Figure 21, by in- 
cluding several corrections previously neglected, which accounts for the a de- 
pendence of the S'-factor, namely the normalization of initial-state Coulomb 
penetrabilities, final state charged particle Coulomb interactions, linear de- 
pendence on a of photon coupling to nuclear currents in radiative processes, 
photon energy and barrier penetration in external direct captures and, finally, 
the electromagnetic contribution to the nuclear masses. All these corrections 
can amount to a third of the total dependence of ^Li on Sa/a, while they are 
smaller for ^He and ^H. 



These results can be used to get bounds on 6a/ a by comparing data with 
the results of a fully numerical computation of light nucle i versu s a, us- 
ing a modified version of public BBN codes. For example (lAveOll ) reports 
at 95% C.L. 6a/ a = —0.007 ± 0.009, using the experimental val ues of ' ^He, 
Yp = 0.244±0.002, and of deuterium, 2H/H=(3.0±0.4)■10-^ while flNol02h for 
the same experimental values obtain 6a/ a = — O.OOTj^^n m?- If we u se the ^H/H 
and Yp of Section 4 and the response matrix formalism of (jDenOTl ) we find the 
updated bound -0.007 < 6a/a < 0.017 at 95% C.L.. It should be noticed that 
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a varying value of a does not significantly improve the global fit of data when 
^Li is a lso included in the analysis. If one uses log]^Q(^Li/H) = — 9.911qj3 
flRvaOOh . combined with the deuterium experimental result as above, one 
would find a reasonable agreement for 6a/ a ^ 0.23 (the standard value being 
excluded at more than 3-0"), but including ^He gives a minimum very far 
from this value but with a very low likelihood flNol02h . Stated differently, it 
seems difficult that the problem of the lower observed '^Li with respect to theo- 
retical expectation could be alleviated by allowing a different value of a at the 
BBN epoch, unless one also includes iVgff as a free parameter in the analysis. 
In this case, a larger value of a, which suppresses the theoretical value of ''Li, 
and a slower expansion rate at BBN, N^s < 3 (which is degenerate with a and 
can compensate for the otherwise too larg e amou nt of ^He for positive 6a/ a) 
can reasonably fit all nuclear abundances fllch04t ). 



8.2.3 The role of Higgs vacuum expectation value, fermion masses and Aqcd 

In the SMPP the vacuum expectation value, v, of the Higgs field provides the 
mass term to vector bosons and Z via the Higgs-Englert-Brout mecha- 
nism, as well as to all fermions via Yukawa couplings. Any variation of the 
weak scale, v, or better to say of the ratio of the weak to strong scales, v/Aqcd, 
and weak to gravitational scale over cosmological times can therefore, have a 
dramatic impact since it induces a time dependence of all massive particles 
after spontaneous symmetry breaking of the electroweak symmetry. 



The main effects on BBN are the change of the Fermi coupling constant, the 
neutron-proton m ass difference, the elec t ron ma ss and finally, the binding en- 
ergy of deuterium fjPixSSl : ISch93l : Ilch02bl : IYoo03| ) . The change oiGp = 1/ \/2v'^ 
re-scales all weak reactions which keep neutron and proton in equilibrium, and 
thus shifts the decoupling temperature, T^. Similarly, these rates are also af- 
fected by the variation of the neutron-proton mass difference, due to the u 
and d quark mass contribution to Am. This contribution is proportional to 
the difference between the two corresponding Yukawa couplings and can be 
estimated once the electromagnetic contribution to Am is singled out. Since 
the latter can be calculated with relative ly little uncertainty, being the Born 
term of the Cottingham formula (jGas82l ). one obtains, see Eq. (120), 



Am(MeV) = 1.29 + 2.053— 



(123) 



with the assumption that in this case only v is considered as a free parameter, 
while both a and the Yukawa couplings are fixed to their standard values. 
This est imate i s consistent with recent lattice QCD calculation with dynamical 
quarks (jBeaOTl ). Finally, the electron mass is shifted linearly with v, 6me/me = 
6v/v. The changes in Gp, Am and me fixes the new value of Am/To, which is 
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a key parameter for Yp, as we stressed already several timea I. For example, 
if V is increased, the lower value of Gp and the larger mg (which reduces the 
available phase space) implies a less efficient n-p rates, so that their freeze out 
occurs earlier (more ''He), but the larger effect is due to the increased Am, 
which reduces the n/p ratio, and therefore the final value of Yp. 



Finally, changing v alters the deuteron binding energy BA I through the 
change in the pion mass m^. Since the pion is a Goldstone boson, its mass 
scales as a geometric mean between weak and strong scales as found by Gell- 



Mann, Oakes and Renner fiGelGSl ). ~ (m„ + mdjAqcD oc v. It is worth 
noticing that a change of Bd can be induced by a variation of f , as we are 
presently discussing, or a time evolving values of Yukawa couplings, yi, with 
i = u,d, or finally by a change of the strong interaction scale, Aqcd, which 
we consider later on, 



Sm^ _ 1 / 6Aqcd _^ ^ ^Vu+^W 
2\ Aqcd V Vu + Vd , 



(124) 



Therefore, the effects which we now describe of a modified value of Bd will be 
relevant when any of these parameters is assumed to be time-dependant. 

The value of Bjj fixes both the initial conditions for BBN, given by the Nuclear 
Statistical Equilibrium, as well as the cross-sections that burn deuterium into 
heavier nuclei. As B^ increases, deuterium becomes more stable, and BBN 
would start at a higher temperature epoch T^. Interaction rates of the BBN 
network occur more rapidly at higher temperature (though the cross-sections 
have changed) leading to a more efficient processin g. One thus expe ct in 
this decr eased ^ H/H ratio and an increased Yp (IKneOSbl : lDmi04l ). see 

also Table I in (jDenOTI ). 



Static properties of nuclei and, in particular, binding energies depend strongly 
on m^r, which sets the length scale of attractive nuclear forces, and gives 
the dominant contribution to two- and three-body potential via the one and 
two pion-exchange terms. Though accurate calcu lations have been per formed 
which reproduce several experimental properties (jPudQTl : iPieOll : iFlaOTI ). a nu- 
merical determination of the functional dependence of binding energies upon 



The change in electron mass also produces a different value for the — e~ ener gy 
density, and thus of the Hubble parameter, but the effect is sub-leading ( YooOd ). 
'^^ Regardless of its expression in terms of fund amental parameters, the effect of 
changing Bo has been studied for example in (jDmiOJ). Combining a full BBN 
data analysis, which also includes ^Li abundance with the WMAP prior on r], they 
find 6Bd/Bd = -0.019 ± 0.005, and point out that this 4-(T shift of the D binding 
energy might reconcile the whole BBN picture with data, in particular the low value 
of observed ^Li. 
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rriTr is still lacking. On the other hand, this dependence has been studied ex- 
tensively in the framework of low energ y effective theory tha t respec t the 
approximate SU{2)l ® SU{2)r of QCD flBul97l : lBea02l : iBeaOSi : lEpeOSh . De- 
pending on the values of the coefficients (fixed by data) which weight the 
s-wave four-nucleon operator expansion of the lagrangian density, the result 
seems to suggest that deuterium remains loosely bound for a wide range of 
mjr, or rather that the value of Bd strongly changes with m^r (IBeaOSl ). Despite 
the large uncertainties, the results are compatible with a linear dependence 
for small variations around the current value, 



BoiMeV) = 2.22 1 



5m, 



(125) 



with -10 < r < -6 flBeaOa lEpeOST ). 



One further caveat is represented by the fact that the role of strange quark 
mass in nuclear quantities such as binding ene rgy is not fully understood, being 
so close to the non-perturbative scale, Agcn fjFla02l : lDen07l ). The sensitivity of 
Bo to rris has been estimated in (IFlaOSl ) . by inspecting the role of a meson mass 
to nuclear potential. They find SBd/Bd ~ —176ms/ms- Some authors have 
also pointed out that large pion mass variation could also render deuterium an 
unbound system, and discussed the values of \ yhich \y ould lead to stable 
di-neutron (or di-proton) systems (jFla02l : lDen03l : |KneOJ). This would have 
dramatic consequences on the standard BBN picture, leading to a stable or 
long-lived ^Be, and by-passing the A = 5 bottleneck through formation of 
^He. Both effects would produce a large enhancement of lithium or metallicity 
production in the early universe. 



Implementing all changes described so far in a BBN numerical code, and 
assuming w as a free parameter, one can typically constrain 5v/v at the level 
of percent. In particular, one can estimate 6Yp ~ Sv/v{%). Using both ^H/H 
with an error of 30% and a rather low value for ¥„ = 0.238±0.005, which seems 



presently disfavored, Yoo and Sherrer find —0.7% < 6v/v < 2.0% flYooOSl ) . 
Interestingly, a looser bound (10%) is obtained from (pre-WMAP) data on 
CMB, by considering the effect of a varying electron mass in both the Thomson 
scattering cross-sectio n and h ydrogen binding energy. We have updated the 
bound from BBN of (lYooOSl ) . using the data on ^H/H and Yp discussed in 
Section 4 and find -0.005 < 5v/v < 0.012 at 95 % C.L. 



The scenario with a varying strong interaction scale, Aqcd, can be described 
in quite a complete analogy. In this case, assuming that the weak scale is 
kept fixed, while v/Aqcd is time dependent, the effects are due to the change 
of neutron-proton mass difference, since the electromagnetic contribution is 
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weighted by the strong scale, 



Am(MeV) = 1.29-0.76 



6A. 



QCD 



A, 



QCD 



(126) 



and the shift in d euterium and heavier nuclei binding energy (jFla02l : lKne03bl ) . 
In particular, in (IKneOSbl ) results are obtained by exploiting a simplified BBN 
network up to nuclei with A=3. This introduces some systematics in the result, 
but it has the benefit of reducing the dependence on Aqcd only through the 
two parameters Am and Bd- The result is a bound on Aqcd at the level of 
10%, if one corrects for such a systematic effect. 



The effect of binding energy on BBN is also considered in (jFla02l ) . where 
the authors quote \SAqcd/ Aqcd\ < 0.06 as a conservative bound. Interest- 
ingly, they also point out that even when both the strong scale and quark 
masses are modified by the same amount, so that all binding energies re- 
main unchanged (and so does the reference temperature, T^r), nevertheless 
the freezing temperature is changed, since it also involves the Planck mass 
scale which is assumed to be fixed. In this case one obtains a bound of the 



order of {6Aqcd/Mp)/{Aqcd/Mp) < 0.1 flFla02D . 



8.2.4 Correlated variation of fundamental constants in unified scenarios 

A summary of bounds on fundamental parameters considered in the previous 
Sections is presented in Table 12. As we mentioned already, they refer to a sin- 
gle parameter analysis, i.e. assuming that all but a single fundamental constant 
are held fixed. On general grounds, one could expect that this assumption is 
rather ad hoc. For example, in models with extra dimensions, or based on em- 
bedding of the SMPP into a larger gauge symmetry group, it is quite natural 
that more than a single coupling or fundamental scale, if not all, would be 
time dependent during the (homogeneous) expanding history of the universe. 
When using BBN to constrain such a variation (and CMB anisotropics to a 
lesser extent as they are not sensitive to strong interactions) , one immediately 
faces the problem of several degeneracies among these parameters. This re- 
sults in a much less predictive power, unless a specific theoretical framework 
is assumed in the analysis, which reduces the number of free parameters of 
the theory, and relates possible variations of different fundamental constants. 

The purpose of this Section is to describe the impact on the BBN of simulta- 
neous variation of all involved fundamental parameters and discuss the con- 
straints on these variations which can be obtained when some particular model 
is assumed (GUT theories, string-theory inspired dilaton scenario). 

The first steps in this programme are basically the following: to identify those 
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Table 12 

A summary of BBN constraints on fundamental parameters discussed in Sections 

8.2.2 and 8.2.3. 

Data Range Ref. 

baja = 0.244 ± 0.002 

2H/H = (3.0 ± 0.4) X 10-5 -0.016 ^ 0.002 (95 % C.L.) (AveOl) 

baja = 0.244 ± 0.002 

^H/H = (3.0 ± 0.4) X 10-5 -0.024 ^ 0.003 (95 % C.L.) (Nol02) 

(5Q/a = 0.250 ± 0.003 this paper 

^H/H = (2.87lg:[i|) X 10"^ -0-015 ^ 0.014 (95 % C.L.) (using fPenOT)^ 

bvjv = 0.238 ± 0.005 

^H/H = (3.010:5) X 10-5 -0.007 ^ 0.02 (Yoo03) 

(Jv/z; Yp = 0.250 ± 0.003 

^H/H = (2.87l|]:g^^) X 10"^ -0-005 ^ 0.012 (95 % C.L.) this paper 

Mqcd/Aqcd Yp = 0.238 ± 0.005 

2h/H = (2.6 ±0.4) X 10-5 ~-0.lH-~0.1 rKne03b'l 

bAqcD/AQCD ^H/H = (1 ^ 10) X 10-5 ~ -0.06 ^ ~ 0.06 (Fla02) 



combination of fundamental parameters which influence the BBN dynamics 
and quantify how light nuclei abundances change when this parameters are 
left as a free input. In general, this task is very involved, as large variation 
of, say, the strong scale parameter, Aqcd, or quark masses can change the 
standard picture of BBN in quite a dramatic way (e.g. by rendering deuteron 
an unbound state or predicting bound (pp) or (nn) states which may add new 
paths to nucleosynthesis). On the other hand, it is easier to study the prob- 
lem perturbatively, in the neighborhood of the standard values that all these 
couplings, mass scales, etc., take today as measured in laboratory. A com- 
prehensive analysis of this sort was lacking until quite recently, and h as been 
mainly pursued in a very detailed study of Dent, S tern and Wetterich (jDenOTI ) 
(see also a semi-analytical tour de force of flLanOGl ) in the chiral limit of QCD). 
We will follow their approach in the following, and choose to work in the units 
for which the strong interaction scale is held fixed (see our discussion in Sec- 
tion 8.2.1). One can then identify the following set of fundamental constants 
playing a role at BBN epoch, which we collectively denote by 99^, k = 1, ...7: 
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the Newton constant, GajEII, the fine structure couphng, a, the Higgs vacuum 
expectation value, v, the electron mass, m j I . the light quark mass differ- 
ence, = nid — mu, the averaged light quark mass, Mg = {rrid + mu) /2 oc m\ 



( iDenOTi ). Finally, the baryon density parameter is also included. Notice that in 
this analysis variation of the strange quark mass, itIs/Aqcd, is not accounted 
for. 

The leading linear dependence of nuclear abundances Xj {i =^H, '^H, ^He, ^Li, 
''Li) upon small changes of these parameters is then encoded in the response 
matrix R defined as 

Rik = ^jr^ . 127 



This matrix can be evalu ated by numerically integrating the BBN equations. 



In particular, in (jDenOTl ) this has been performed in two steps, by first vary- 
ing a set of nuclear physics parameter Vj which includes binding energies of 
light nuclei up to '^Be, nucleon mass, neutron-proton mass difference, neutron 
lifetime and a subset of the (fk, («, f^e, Gn and 77) and computing the matrix 

a, = ^?^ , (128) 



Xj dvj 



and then relating the variation of the Vj to small changes of the ipk-, 

f.. = ^|^ . (129) 

The response matrix R is therefore given by the matrix product R = CF. 
This decomposition turns out to be useful as the determination of C is rather 
robust, while computing F requires some theoretical assumption, as for ex- 
ample the dependence of binding ene rgies on the pion mass. In Table 13 the 



entries for the R matrix are reported (jDenOTl ). The last row is the logarithmic 



variation for small changes of the baryon fraction around the reference value 
VLsh"^ = 0.022. We have checked that the result is very weakly dependent on 
the reference point in the range 0.015 ^ 0.025, the variation being of the order 
of 5% for and even smaller for ^He. 



The results of this Section have some overlap with the analysis presented later of 
scenarios when one varies the gravitational action, including a time-evolving Newton 
constant. 

Differently than in our previous discussion, me is assumed an independent param- 
eter with respect to v. This is equivalent to consider the lepton Yukawa couplings 
as varying parameters. 
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Table 13 

The response matrix R for ^H, ^He, ^He, ^Li and ''Li (from ( DenOTi )). The reference 
value for the baryon fraction is Q^/i^ = 0.022. 



^He "^He 



^Li 



^Li 



Gn 


0.94 


0.33 


a 


3.6 


0.95 


V 


1.6 


0.60 


me 


0.46 


0.21 




-2.9 


-1.1 




17 


5.0 




-1.6 


-0.57 



0.36 


1.4 


-0.72 


1.9 


6.6 


-11 


2.9 


5.5 


1.7 


0.40 


0.97 


-0.17 


-5.1 


-9.7 


-2.9 


-2.7 


-6 


-61 


0.04 


-1.5 


2.1 



Using these results one can construct the function (we do not use ^Li in 
this analysis) 

X — 2 ' 2 ' 

0-2 H 



where the theoretical abundances are computed using the linear expansion 
in terms of R (but we have retained an exact dependence on the baryon 
fraction). In this way, one gets an estimate of the kind of constraints which 
can be obtained on all (pj, and more importantly, it serves to illustrate the 
degeneracies which appear between the several pairs of fundamental constants, 
at least as long as we consider small variations with respect to the standard 
results, so that the linear expansion used in Eq. (130) is legitimate. 

As an example, we have considered the experimental values and errors of 
Section 4. The results are shown in the figures of Tables 14 - 15, where the 
various curves correspond to a Fisher matrix analysis around the reference 
model with fi^/i^ ~ 0.0209 and all other parameters at their standard values 
(this is indeed, locally, the minimum of the total x^)- When one of the two 
selected parameters is the baryon density, the plots correspond to the case 
where only one fundamental constant is varied while all others are fixed to the 
standard value. In this case, one can read from the contours the typical order 
of magnitude of the constraint which can be obtained in this single-parameter 
analysis, which is at the few percent level for a, v, Aq and Mg, while it is almost 
one order of magnitude larger for Gn and me. More interestingly, the other 
bidimensional contours show the degeneracies among pairs of parameters. In 
particular, notice the strong correlations of Gn with a and nie and between 
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Table 14 

The bidimensional 68 and 95 % C.L. contours illustrating the correlation of funda- 
mental parameters in a BBN analysis (^H/H and Yp as in Section 4). 





-0.1-0.05 0.05 0.1 -0.1-0.05 0.05 0.1 -0.1-0.05 0.05 0.1 

ha la 5a /a 8v /v 



the pairs v — a, rrie — a, Aq — a, v — rrie, and v — Ag. 

The difficulty of treating simultaneous variations of more than a single fun- 
damental parameter is somehow alleviated if one works in the framework of a 
definite theoretical model beyond the SMPP, which imposes further relation- 
ships among (some) of them. A typical example is given by Grand Unified 
Theories (GUT) which assume that the three SMPP gauge couplings get uni- 
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Table 15 

Figures of Table 14 continued. 
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fied at some high mass scale, M^, where the dynamics is dictated by a larger 
(simple) symmetry group. T his sceriario and its irnplications for BBN ha s 
been worked out in details in (jCal02al : lLan02l : ICal02bl : lDen03l : iDenOTl : ICocOTI ). 
The general argument is that in GUT theories unification of couplings implies 
that the various fundamental constants are likely to vary simultaneously and, 
for example, a change of the fine structure coupling implies a much larger 
variation of the strong interaction scale, Aqcd, at low energies. 

If we denote by ac the (common) value of the three SU{3)c x SU{2)l x U{1)y 
couplings, aj, at Mq (typically Mq ~ 10^^ GeV), then their running is given 
by Renormalization Group equations. 



hi , (Mg 



(131) 



where the coefficients hi depend on the particle multiplets. If ac undergoes 
time variation at a cosmic time above the unification scale, then one can traces 
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the correlated variation of th e SMPP couplin gs at low energies. Neglecting the 
threshold effect one obtains (jCal02al : lLan02l ) 



SA 



QCD 



A, 



QCD 



(34^40)^ 

a 



(132) 



Also Yukawa couplings and the Higgs vacuum expectation value v (which is for 
example, tied to the supersymmetric breaking scale in most supersymmetric 
models) are expected to vary and to be related to time evolution of a. However, 
the size of the effect is model dependent, or simply unknown, so it is usually 
parameterized in terms of an unknown constant, St/i/t/i = c5a/a, and similarly 
for V. In general, in addition to the time evolution of ac, one cou ld also c onsider 
other possibilities, where the GUT mass scale Mq also varies flCal02bl l 



The BBN constraints in these scenarios on the variation of a (and all o ther re- 



Lan02 



Coc07 



lated p arameters) have been first analyzed in a qualitative manner in ( 
Deii03) and then worked out using a full BBN numerical study in ( 
DenOTl ). with two main results. On one hand, the constraints become typically 
tighter, |(5a/Q!| < 10~^ — 10~^. Moreover, in some cases the theoretical value 
of ^Li is lowered and can be rendered compatible with the experimental re- 
sult. This is the case, for example, when the weak scale, v, is determined by 
dimensional transmutation, so t hat variation of the largest top quark Yukawa 
coupling induces changes of v fjCam95l : ICocOTI ). or assuming that the varia- 
tion of fundame ntal constants are triggered by an evolving dilaton scalar field 
(Ilch02bl : ICocOtI i 



8.2.5 Varying the Newton constant and scalar-tensor theories of gravity 

The effect of a varying gravitational constant, Gn, on the BBN is through 
the expansion law, H oc v^GJv, which determines both the neutron-proton 
density ratio at freeze-out (Yp), and the efficiency of burning abundance 
when nucleosynthesis starts. We have already seen in the previous Section 
the typical constraints which can be put on SGn/Gn, of the order of few 
percent. In particular, using the values of Yp and ^H/H already discussed, we 
find -0.036 < 5Gn/Gn < 0.086, at 95 % C.L., which is dominated by the 
effect on ^He. In this case, the limit can be read off directly from the bound 
on Ncs, since a change of G^v is equivalent to the change 



- ^ AM 

g7 " ■ 



(133) 



Recently, the effect of a varying G^ has been considered by several authors, 
with similar results, depending on the experimental va lues fo r light n uclei 
adopted in their analysis (for earlier studies see e.g. (jYan79l : lAcc90h l In 
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( 1Cop04l ) it is stressed the sensitivity of to the value of the Newton con- 



stant. In the pessimistic case of a large systematic error on Y„, deuter ium can 



provide a 20 % bound. On the other hand, the results of (jCyb05l ). which 
adopt Yp = 0.249 ± 0.009, suggests that Yp still gives the strongest pos- 
sible constraint. If one makes the assumption of a monotonic behavior for 
the time dependence of the gravitational constant. Gat ~ t~^, the bound on 
6Gn can also be cast in a constrai nt on the exponent x and thus on the 
present value of Gn/Gn- In fICvbOSh it is found -0.0029 < x < 0.0032, and 
-2.4 X 10-^Vr"^ < Gn/Gn < 2.1 x 1Q-^^ jt-\ Our estimate on 5Gn/Gn 
translates into -0.6 x IQ-^^ji-^ < Gn/Gn < 1-6 x 10" 
13.7 Gyr for the lifetime of the Universe. 



-13-yJ, 1 ^gi]2g 



An interesting class of models which allows for a time depende nt grav i tationa l 
coupling is represented by scalar-tensor theories of gravity (jJor49l : iBraGll ). 
where the latter is mediated, in addition to the usual spin-2 gravity field, by a 
spin-0 scalar ip which couples universally to matter fields. Its dynamics define 
the evolution of an effective Newton constant for ordinary matter, to which 
(p is "universally" coupled, as one can read by the action of the model in the 
"Einstein" frame, 



S 



d'^x 
167rG, 



-g [R - 2g^,d^d^^ - V{^)] + F-^i^)g^,; ^ , (134) 



where G^, is the bare gravitational coupling, V{(f) the scalar field potential, 
and F{(p) a positive function, which enters the coupling of ordinary matter 
fields \1/ to the metric. Strong constraints on the present (de noted b y the index 
0) values of the (generally) if post-Newtonian parameters ( WilOll ). 
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(135) 



where 



dlogF 



-1/2 



dip 



da 
dip 



(136) 



are set by solar system experiments, as the shift of the Mercury perihelion 
or the Shapiro delay of radio signals from the Cas sini spacecraft as it passes 
behind the Sun, 70 - 1 = (2.1 ±2.3) x 10"^ flBer03al f^. The value of ao should 
be thus, very small, while /3o can still be large. 



This bound is also usually presented in terms of a lower limit on the Brans-Dicke 
parameter loq > 40000, the constant which weights the scalar field kinetic term in 
the Jordan frame and in the minimal model with Vi^p) = 0, F{p) = p>. 
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Fig. 22. Bounds on the post-newtonian parameters oq and (5 from ^He in 
scalar-tensor theories in the case of a massless dilaton with quadratic couphn g, 
F{ip) = exp(-/3(/92)_ rpi^g Yg^i^g Qf baryon density is Q.Bh^ = 0.0224. From (jCocOd ). 



The implications of models described by the action in Eq. (134) have been 
discussed by several aut hors with different choices of the arbitrary functions 



Fi^) and Vi^Y see e. g. flCas92d: ICas92bl: ISerOi ISerOd: lEto97l:ISan97l : 
AlvOll: ICheOll : ISerOi lKne03al : lNag04l : IChOSl : IPetOSi : ICocOGI : IPeFoi : 
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LarOGb 



NakOd ). In the case of a standard Brans-Dicke theory the value of primordial 
^He produced during BBN can be evaluated semi-analytically as a function 
of the Brans-Dicke parameter uj (which weights the scalar field kinetic term 
in the Jordan frame) ( AlvOll ). in case one consider the particular solution 
if = const during the radiation dominated epoch, giving the bound uq > 
100. The fact that the BBN constraints on the model depend strongly on 
the particular scalar-tensor theory considered in the ai ialysis , thus providing 
quite different bounds on ljq, have been stressed in (jSer02l ) . while a quite 



genera l numerical code fo r BBN in scalar-tensor models is described in ( 1Cli05 



CocOGI ). In particular, in fjCocOGi ) it is studied in details the case of a massless 
dilaton with a quadratic coupling, V{ip) = 0, F{ip) = exp(— /3(y9^), including 
the mass threshold effects when the universe cools down. The latter is due to 
the variation of the trace of the energy momentum tensor of ordinary matter- 
radiation, p — 3P, whenever a single specie becomes non-relativistic, which 
changes t he sour ce term in the Klein Gordon equation for (f. The results of the 
study of (ICocOGI ) show that the strongest bound comes from ^He and indeed 
it is stronger than from solar system experiments, ao < 10~^ for j3 > 0.3, 
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Fig. 23. Bounds on the parameters a and A which give acceptable results for ^H/H, 



Yp and ^Li/H for a scalar-tensor theory defined by a 
initial conditions (pre-BBN) for ip are (fin = —1.3, (pi 



and V (ifi) = A ^^. The 



0. From (jLar06bl ). 



see Figure 22. Finally, we mention that particular choices of the scalar field 
potential V{ip) (and of the initial conditions for ip before the BBN) can solve 
the lithium problem, yet leading to values of ^H/H and Yp compatible with th e 



experimental values. One example is illustrated in Figure 23, from fiLarOGbf ). 
where the potential is chosen to be V{(p) = A^cyj^, and F{ip) = exp(— 
We see that there is a region in the parameter space which corresponds to 
primordial abundances of all light nuclei in agreement with experimental data. 



8.2.6 A varying Cosmological Constant: Quintessence models 



Observ ations of ty pe la sup er novae (lAst05l : lRie04j ). structure formation (ICol05 



TegOGi ) and CMB (iDunOSi ) provide quite a strong evidence for an accelerated 
expansion of the Universe at recent times. The simplest explanation of this re- 
sult is to invoke a new component of the total energy-momentum tensor in the 
form of a cosmological constant A (ACDM model), which is a nice fit of data, 
yet it gives rise to two related fine tuning problems. On one hand, the value 
of A happens to be extremely small with respect to the typical energy scale 
of any fundamental physics model, as the Planck mass (122 orders of mag- 
nitude larger), or supersymmetry breaking scale, or finally, the electroweak 
mass scale (54 orders of magnitude off). Moreover, it seems really a coinci- 
dence that the energy density stored in the form of A and matter appears to 
be of the same order of magnitude just today, Qa ~ 0.75 and Qb+dm ~ 0.25. 
If the cosmological constant had been larger it would have changed the whole 
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structure formation history, while a smaller value would have been irrelevant 
for observations. 



To solve (or to alleviate) these two related problems, many attempts have been 
made in the last two decades and based on the idea that the "dark energy" 
whose present value is given by Q\, is dynamically changing and is due to 
the evolution of a s calar f i eld Q ( the qu i ntessence, k-essence , etc. fi e ld) see 



e.g. flWetSSi: iRatSSi : ICalQSl : ICopQSl : iFerQSl : lLid99l : ISte99l : IZlaQOl : iBraOOl : IChiOO 



ArmOOl : lArmOll ). In particular, depending on the choice of the field potential, 
typically modeled as an exponential, V{Q) = Vq exp(— AQ/Mp) or an inverse 
power, V{Q) = \A^'^°-/Q°-, the field Q evolves according to an attractor-like 
solution of the equation of motions. This means that for a wide range of 
initial conditions at early times, the evolution of Q rapidly converges towards 
these solutions (tracker solutions), which lead naturally to a crossover from 
the radiation dominated epoch to a dark energy dominated era at late times. 



If the Q-field provides a non-negligible contribution to the total energy density 
during radiation dominated regime, it follows that its early dynamics can be 
constrained by nucleosynthesis, and later on by CMB power spectrum. A 
first rough estimate can be obtained by requiring that at the neutron-proton 
decoupling temperature To ~ MeV, the scalar field energy density Qq should 
be small enough not to disturb the even tual am ount of frozen neutrons, i.e. 
of the final ^He mass fraction, see e.g. (IFerQSl ) . This bound can be cast in 
terms of the largest acceptable value for deviations of the effective number of 
neutrino from its standard value at that particular temperature, as in general 
the Q energy density would not scale simply as radiation 



7AiVefr/4 



10.75 + 7AiVeff/4 



< 0.09 



(137) 



where we have used our bound at 95 % CL. on Ncs of Section 5. 



More detailed analysis have been presented in (jYah02l : iKneOSall. In (jKne03af l. 
the authors consider as a model the particular potential of ( AlbOOl ). 



V{Q) = (Q - Qo)' + A exp(-Ag) 



(138) 



and obtain the bound on the quintessence contribution to the total energy 
density during BEN fig < 0.12, A > 5.7 at 99 % C.L.. In flYah02h the analysis 
is performed for bot h the originally proposed inverse power potential of Ratra 
and Peebles ( iRatSSi ). 



(139) 
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Fig. 24. Allowed values of the Q field potential parameter a and initial Pq/pb 
at z = 10^^ (B stands for "background", i.e. ordinary radiation and matter). 
The plots refer to power law (a) and SUGRA corrected potentials (b). Models 
in which the tracker solution is obtained by the BBN epoch are those in the band 
denoted by "Tracker Solution at BBN". Values of a to the right of the lines la- 
beled wq = —0.6, —0.4, —0.2 are excluded by requiring that the present equation of 
state be sufficiently negative. The BBN constraint for a maximum energy density 
of the quintessence field o f 0.1 % and 5.6 % are shown as dotted and solid lines, 
respectively. Prom ( Yah02l ) . 



and a modified version of it based on the hypothesis that the quintessence 
field be a part of supergravity models. In this case, for a flat Kahler potential, 
the potential receives an extra factor exp(3(5^/2Mp). Their results are shown 
in Figure 24, for the following choices of Yp and ^H/H 



0.226 < Yp < 0.247 , 
2.9 X 10^^ < ^H/H < 4.0 X 10"^ . (140) 

The plots show the bound on the initial ratio of the Q field to background (B) 
(ordinary radiation and matter) energy density at initial redshift z = 10^^, 
assuming equipartition of Q energy at that epoch, = V{Q), versus the 

potential parameter a. The regions denoted by "Tracker Solution at BBN" 
denotes models in which the evolution of Q has already reached the tracker 
solution. The main effect of BBN is to exclude a large family of possible kinetic- 
dominated solutions with a Q-energy density exceeding that of relativistic 
species prior or during nucleosynthesis. 



103 



8.3 Miscellanea 



8. 3. 1 Testing Friedmann equation 

We have already stressed several times that changing the expansion history 
in the early universe has a big impact on the BBN predictions for light nuclei 
abundance, in particular ^He. Several modifications to the standard Friedmann 
equation have been already discussed, arising in a variety of contexts such as 
non-standard theory of gravity, time evolving Newton constant, or brane- world 
models. One may also try to perform somehow a blind test of the validity of 
Friedmann expansion law without a particular theoretical framework behind, 
and using a suitable parametrization of possible deviations from the standard 



behavior. This has been considered for example in (jCar02l ). where the Hubble 
factor is expressed in terms of two parameters. 



where Hi and the exponent a should be constrained by data, once we fix 
the bar yon density parameter, rj. A similar analysis has been performed in 



( lMas03l ) . with a slight different definition of H{T). Using Eq. (141) one can 
compute the values of the three leading quantities which enter BBN dynam- 
ics: the freeze-out temperature. To, the time elapsed from To and the onset of 
BBN at Ttv, and the value of the expansion rate at T/y. Interestingly, Hi and 
a show a large degeneracy. The same Helium mass fraction and ^H abundance 
can be obtained increasing at the same time the Hubble rate at every tempera- 
ture (i.e. Hi), which raises the freezing temperature and thus leads to a larger 
initial neutron to proton density ratio, and the value of a (for a definite Hi a 
larger a means a relatively lower expansion rate at Tjv). Furthermore, for each 
nuclide there are two branches in the Hi-a plane for which Yp or ^H/H are the 
same, see Figure 25: a lower branch, where also the standard result, a = 2, 
Hi ~ MeV, is lying, and an upper branch, for higher values of Hi, which how- 
ever cannot simultaneously fit both ^He mass fraction and deuterium. The 
whole compatibility region is thus characterized by a single (almost) linear 
behavior, \og{Hi) oc a, with 1.5 < a < 3. This bound can also be used to 
constrain possible non- universal coupling of gravity to the three neutrino gen- 



erations as considered in (IMasOSl ) . where the study is also extended to the case 



of a degenerate BBN, or the p ossibility that matter and antimatter may have 



different couplings to gravity (jMas04l ) . A different gravitational coupling to 



bo sons and fermions, Gn,b and Gn,f respectively, has been instead considered 



in flBar04br ). with the result 0.33 < Gn,b/Gn,f < 1-10 at 2 - cr. 



Similar analyses have been also applied in the framework of metric-affine grav- 
ity models, where the non-Riemaniann effects are encoded in a fictitious fiuid 
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Fig. 25. Contours of constant helium (top) and deuterium (bottom) in the a — Hi 
plane. Chosen values are Yp = 0.24 and ^H/H=3 • 10~^. S quares and filled circles 
correspond to r] = 10~^ and r] = 10"^'', respectively. From ( Car02l ) . 



with equation of state P = p ( iKraOSi ). or which account for the quantum 
gravity corrections for matter fields computed in a lo op quantum gravity ap- 
proach, with a non-canonical equation of state (IBojOSl ). Finally, one can bound 
the value s of para meters w hich measure departure from Lorentz invariance 
( iLamOSbl ). see e.g. (jKosOJ) for a review of the Lorentz- and CPT- violating 
extension of the SMPP. 



105 



The sensitivity of BBN to any modification of the Friedmann equation can 
be also exploited to test a specific prediction of general relativity, which has 
no analog in Newton theory, namely the fact that pressure contributes to the 
acceleration of the scale factor, 



a 
a 



4tiG 



N 



(P + 3P) 



(142) 



Indeed the Friedmann equation knows about the pressure contribution in Eq. 
(142) via the Bianchi identity. If one allows for devi ation f rom this distinc- 
tive feature of Einstein theory and introduces, as in (jRapOSl ) . a free constant 
parameter, which modifies the spatial component of Einstein's equations. 



a 
a 



N 



(p + ^xP) 



(143) 



one obtains during a radiation dominated expansion (neglecting the spatial 
curvature term). 



1 + X 8ttG 



N 



(144) 



Bo unds on x can then be obtained from the BBN. Not surprisingly, the results 
of (IRapOSl ) show an excellent agreement with the standard expectation, x = 1- 
Notice that variation of this parameter is completely degenerate with a change 
of the Newton constant during the BBN, or a non-standard effective number 
of neutrinos Ncs. 



8.3.2 Primordial Black Holes and BBN 



Primordial black holes may form in the ea rly universe in presen ce of sub- 
horizon density perturbations of order unity (jZelGTl : lHaw7ll : ICar75l ) . Assessing 
cosmological constraints upon their density over some mass ranges can provide 
important in formatio n on the primordial density fluctuations. For a recent re- 
view see e.g. (IKhlOSbT) . These c onstra i nts haye been d i scussed by many author s 



view see e.g. iiKmuaPi) . i nese c onstra i nts nave Deen a i scussed py many autnor ! 
in the 1970s flHaw?! ICar76l : IZel76l : IZelTj- IVai78al : lMiv78l : IVai78bl: IhinSOl) 



while updated analysis using BBN flKoh99h and CMB data flRic07l : iTasOSh 



have been presented quite recently. Primordial black holes evaporating by the 
BBN epoch {t < 10'^ s) have a mass lying in the range M < 10^° g, s ince there 



is a simple relation between mass and lifetime, Tb/i, see e.g. (jKoh99l ) 



M ~ 10^ f —\ 
Vsec/ 



1/3 



(145) 
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Fig. 26. The BBN bounds on the primordial black hole fraction versus their mass. 
Regions to the rig ht and above the curves are excluded by primordial nucleosynthe- 
sis. From ()Koh99l l. 



They emit various particles such as neutrinos/antineutrinos, photons, quark- 
gluon jets which produce hadrons through fragmentation processes. All these 
(high energy) particles interact with species already present in the thermal 
bath, and induce several effects which can change the standard picture of 
light nuclei production. For example, production of high energy neutrinos and 
antineutrinos change the weak interaction freeze-out temperature, and thus 
the neutron to proton ratio. Similarly, large hadron injection after the weak 
process decoupling temperature may revival chemical equilibrium between 
nucleons, thus leading to quite different amounts of ^He and deuterium with 
respect to the standard case. 



In Figure 26 we report the results of fjKoh99l ) for /5(M), the primordial black 
hole initial fraction to the total energy density as function of their mass. The 
bound for 10^ g < M < 10^° g, /9(M) < 10~^° is strongly constrained by 
the upper bound on Yp (they use Yp < 0.252 in their analysis), since in this 
mass range the effect of black hole evaporation is to delay the freeze out of 
n-p chemical equilibrium. On the other hand, for larger masses and lifetimes. 



107 



the extra produced neutrons cannot be burn into ^He and thus contribute 
to deuterium abundance, which provides a stronger constraint in this region, 
(3{M) < 10-22. 



8.3.3 Mirror world 

The theory of a hidden mirror world, an exact duphcate of our visible world, 
is based on the product of two identical gauge symmetry group GxG', with G 
the SMPP group, SU (3)c x SU {2)^x11 (l)y , in the minimal case. The two sec- 
tors communicate through gravity and, possibly, by other interactions (kinetic 
mixing of photons and mirror photons^, neutrino mixing, common gauge 
symmetry of flavor) and are exchanged by the action of a discrete symmetry, 
the mirror parity, which implies that both particle sectors are described by 
the same action, and are characterized by the same particle content. For a re - 



view and a comprehensive list of references we address the reader to flBer03bl ) . 
The mirror world would influence the whole evolution history of the universe 
through its contribution to the expansion rate, and can represent a natural 
candidate for dark matter, as well as providing a mechanism for baryogenesis. 
Of course, if mirror particles populate the primordial plasma with the same 
densities of ordinary particles, they would contribute to the effective neutrino 
number for a too large value Nes ~ 6.14 at the onset of BBN, which is ex- 
cluded by light nuclei abundances. Thus, the mirror particle density should be 
reduced, and characterized by a plasma temperature, T', lower than the ordi- 
nary photon temperature, in order to be compatible with the bound on AN^s. 
This can achieved if the inflationary reheating temperature is different in the 
two sectors, and all possible interactions between the two worlds are weak 
enough to ensure that they do not come into mutual thermal equilibrium. 

The ratio T'/T can be computed by using entropy conservation. Use of the 
definition of Section 2 leads to 

T{t) \.,) \g..(T)) -''\g..(T)) ' 



with s {s') the entropy density of ordinary (mirror) species. During the radi- 
ation dominated epoch we have 



H{T) 



3 30 



g.{T)T' + g,{T')T'' 



1/2 



2^ BBN and CMB bounds on the mixing of phot on with a hidden light abelian gauge 
boson have been recently considered in ( JaeOSi ) . 
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Fig. 27. The primordial mirror ^He mass fraction as a function of x (see Eq. 146). The 
dashed curve is the approximate result of Eg . (149). The solid curves are obtained 
via exact numerical calculation. Prom ( BerOll ). 



3 30 



nl/2 



(147) 



where the last approximate equality holds as long as the value of x is not too 
small. On the other hand the parameter x can be re-expressed in terms of 
AAi'cff. Since photons, electron/positron pairs and active neutrinos correspond 
to g^: = 10.75 we have 



AN, 



cS 



6.Ux^ 



(148) 



and thus the BBN bound AA^cfr < 0.4 gives x < 0.51. 



It is interesting to sketch how nucleosynthesis proceeds in the mirror world, in 
particular to compute the final yield of mirror ^He. As usual one can use the 
limit that all (mirror) neutrons get bound in ^He and compu te the decoupling 
and deuterium formation time. In this case we have (IBerOll ) 



2 exp 



1 + exp [Am/TDil + a;-4) V6] 



(149) 



where it has been used the fact that the weak interaction decoupling tempera- 
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ture in the mirror world is larger than Td by a factor (l+x""*)^/^ and that, since 
T'p^ ~ T/v unless the mirror baryon density parameter t]' is very different than 
, the t ime of nucleosynthesis is t{T' = T^) = t(TAr)/(l+x~'^)^/^. As noticed in 



BerOll ). the estimate (149) is not valid for small x and rj' = 10~^° — 10~^ ~ 77, 
since in this case deuterium production can become ineffective, thus inhibit- 
ing the whole BBN r eaction network. The result of versus x as numerically 



computed in (iBerOll ) is shown in Figure 27. Notice that for large rj' > 10~ 
the value of is in good agreement with the estimate provided by Eq. (149). 
Indeed, this large rj' /rj regime is particularly interesting, since it corresponds 
to a sizeable contribution of mirror baryons to the present energy density of 
the universe, 

Tr=^- ■ (150) 
VLb T) 



In this case mirror baryons might constitute the dark matter or one of its 
components I. In view of the bound on x, this requires a large value for the 
ratio rj' /rj > 10, so that the mirror world would contain a considerably bigger 
fraction of ^He than the visible world. Further details on the th ermody namics 
of the early Universe with mirror dark matter are available in (jCiaOSi ). 



9 Massive Particles &; BBN 



The existence in the primordial plasma of one or more species of massive par- 
ticles {mx ^ rrie) besides baryons is at very least a likely possibihty, given the 
numerous pieces of evidence in favor of the existence of dark matter (DM). A 
wide variety of observations suggests that most of the matter in the universe is 
not in a visible form (with Qdm — ^^b), and several of them also imply that 
DM is non-baryonic and cold, i.e. made by particles with a non-relativistic 
momentum distribution. Direct and indirect evidence, starting more t han sey - 
enty years ago with the orbital velocities of galaxi es with in clusters ( lZwi33l ). 
now includes also t he rotational speeds of galaxies (IBorOOl). gra vitational lens- 
ing ( Dah07| : CloOG ). the cosmic microwave background (Dun08), th e large scale 
structure (iTegOa ) . and the light element abundances themselves flOliOOl ). 



At first sight, it might appear that DM has no impact on BBN since the 
expansion of the universe at the BBN epoch is dominated by radiation, and the 
presence of additional massive particles at the BBN ep och is thus dyn amically 
irrelevant. We recall that, for example, the analysis of (jCar02l : iMasOSi ) without 



The growth of perturbations and the power spectru m in presence of a sizea ble mir - 
ror baryon density has been studied in (jBerOll : llgnOi ^ and in more detail in 
and reference therein. 



110 



assuming strong priors on rj show that the behavior of the Hubble parameter 
should be very close to the radiation dominated regime during the BBN, see 
Section 8.3. Once the WMAP prior on rj is used, this resuh would narrow 
even further. Thus, even from an observational point of view, any effect due 
to massive particles at the BBN must be of non-gravitational nature. 



Additional couplings for the DM particles are far from being a remote possibil- 
ity: there are instead strong motivations which point to a connection between 
dark matter and the electroweak scale. A stable particle with an electroweak 
scale mass and couplings would naturally be produced in the thermal bath 
of the e arly universe in an amount similar to the observed matter abundance 
( lLee77l ). From a particle physics perspective, the hierarchy problem appears to 
require new physics at or above the electroweak scale. Furthermore, stringent 
constraints from electroweak precision measurements (and from the stabil- 
ity of the proton) indicate that these new particles respect symmetries which 
limit their interactions. Such symmetries can also lead to the stability of one or 
more of the new particles, such as the lightest exotic particle in the spectrum of 
R-parity conserving supersymmetry, K-parity in universal extra dimensional 
mod els, or T -parity in little Higgs models. For a review of DM models, see 
e.g. (lBer04bl ). 



Dark matter cannot bring electric charge by definition (otherwise it would be 
visible) and — although th e reason is less obvious — it appears that it cannot be 
strongly interacting (see (jMacOTl ) and reference therein for a recent overview 
of the stringent bounds on such scenarios). However, DM might be produced 
as a decay product of charged (or strongly interacting) progenitors, and there 
are also viable particle physics scenarios where this can tak e place , the most 
popular of which being the so-called super- WIMP scenario flFenOSl ). Charged 
or strongly interacting metastable particles at the BBN epoch may form bound 
systems with baryons, altering the nuclear reactions pattern and thus the 
yields of light elements. The analysis of this scenario, which goes under the 
name of "catalyzed nucleosynthesis", is addressed in Section 9.2. On the other 
hand, particles which decay (or annihilate) into Standard Model products, 
even if only weakly interacting, may affect BBN via the cascades they induced 
in the plasma: the injection of secondaries triggers directly or indirectly non- 
thermal reactions, altering again standard BBN predictions. Of course, if the 
lifetime satisfies rx -C 1 s, the particles have decayed well before the BBN and 
no meaningful bound can be derived. 



In general, annihilations have a similar effect as decays, the main difference 
being the time (or redshift) distribution of the injection: for decays the rate 
per unit volume is ux/tx, while for annihilations is given by nxnx{av). If 
we parameterize (av) oc an{T/mx)^ (n = 0, 1 correspond then to s-wave and 
p-wave annihilation respectively), in the early universe when homogeneity is 
a good approximation, the former gives an injection rate scaling as (1 + z)^, 
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or (1 + 2;)^+"' for annihilations. In the following, unless otherwise stated, we 
will generically refer to "decays" to indicate any mode of injection of energetic 
particles in the plasma, thus implicitly including the annihilating case. Some 
specific features of the annihilation injection mode will be discussed at the 
end of Section 9.1.2. 



9.1 Cascade Nucleosynthesis 



The phenomenology of non-thermal BBN bounds strongly depends on the 
branching ratios of the secondaries the X particle decays into: the hadronic, 
electromagnetic (photons, e^), neutrino or inert (exotic invisible) ones. The 
last case requires two or more (meta) stable particles, which does not happen 
very naturally in most realistic models of physics beyond the SMPP. Perhaps, 
one exception is provid ed by t he case of the axino decaying into a gravitino 



and an axion, see e.g. (lChu93l ). In general, the only constraint in this case 
comes from the requirement that the universe at T ~ 0.1 MeV is radiation- 
dominated, which implies 



We do not tr eat thi s scenario and its constraints further. For more details, see 
for example fISchSSl l 



It is worth noting that even a decay into neutrinos X Y+u has effects on the 
BBN. The effects are two-fold: i) energetic neutrinos can create charged lep- 



tons by annihilation onto the thermal neutrino background I: although with a 
suppressed efficiency, sufficiently high energy neutrinos can also produce pions 
via z/ + z/j/i ^ 7r+ + vr^ which affect the p — n equilibrium, ii) electromagnetic 
(and possibly hadronic) showers are induced by 3 or 4 body decay channels 
via a virtual or real weak boson propagator, like X ^ Y + u + + e~; the 
previous final state is always considered as kinematically allowed, when a tree 
level X — > y + z/ is included in the analysis, since — given the binding ener- 
gies of nuclei — to induce any change to BBN the phase space available must 
be ^ 1 MeV anyway. The importance of these sub-leading channels for phe- 
nomenological constrai nts from astrop hysical arguments has also been emp ha- 



sized in other contexts (jKac07l : lMac08l ). We address the reader to (IKanOTbl ) for 
a more through analysis. In the remaining of this Section we describe electro- 
magnetic cascades (Section 9.1.1) and hadronic ones (Section 9.1.2). It is worth 
noting that, for a very broad range of lifetimes involved (0.1 s < rx < 10^^ s). 



In principle, they can also upscatter in the plasma, but decays are only relevant 
when T <C 1 MeV, when almost no charged particle populate the plasma. 
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the thermalization of the secondaries takes place in a time which is neghgible 
with respect to the Hubble time, and therefore redshifting of particles can be 
safely neglected. 



9.1.1 Development of the electromagnetic cascade 



There are some features of the electromagnetic cascades in the primordial 
plasma which significantly simplify the treatment with respect to hadronic 
ones. When the injected e~^,e~,'j are energetic enough, the cascade develops 
very rapidly by a combination of two processes: the pair-production on the 
CMB thermal distribution (7 + 7cmb e"*" + e~) and the inverse Compton 
scattering of the non-t hermal electrons and positrons (e^ + 7cmb — >■ -|- 7) 
off the CMB photons flAha85l : IZdzSOl : ISveQOl : |Pro95l : iKawQSh . There is a crit- 
ical energy, Ec{T), above which the non-thermal spectrum is quickly cutoff 
by pair-production, lea ving virtually no photon available to photo-dissociate 



the few available nuclei I. Although the typical energy for pair-producing 
on the bulk of the CMB distribution is ml/ {Equb) — ml/2.7T, since reac- 
tions on the energetic tail of the distribution are also important . Eg tu rns 
out to be smaller: Nume rical calculations estimate Ec — ml/22T ( iKawQSl ) or 
Eq — mg/23.6T (jPro95l ) . By equating Ec to the and '^He binding energies 
of 2.2 MeV and 19.8 MeV respectively, one infers the keV-scale characteris- 
tic temperatures below which a small but non-negligible fraction of 7's (at 
the percent level) is available to photodisintegrate the light nuclei. This is 
why meta-stable particles with a too short lifetime cannot significantly affect 
thermal BBN yields. Meaningful bounds on electromagnetic cascades are only 
achieved if tv > 10^ s for deuterium dissociation and tx > 10'' s for the more 
tightly bound ^He (in a radiation-dominated universe t oc oc E^). 



Another interesting feature is that the e.m. cascade develops so rapidly that 
most of the effect depends only on the total amount of injected energy Eq and 
the time of injection, rather than the nature and energy of the primary. It is 
customary to define a parameter representing the energy stored in meta-stable 
particles before the cascade begins (hence the zero subscripts), for example 
(x = 'n^x^xfl/n^fl, which together with the decay time tx characterizes al- 
most completely the process. A quasi-universal shape of non-thermal photons 
is reached very quickly. Numerical simulations have found a spectrum below 



'No e+,e~ remain instead available for electro-disintegrations, since the cross- 
section with the abundant CMB photons is effective to quickly cool them down to 
thermal energies. 
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Fig. 28. Number of destroyed nuclei per TeV of electromagnetically interacting 
energy injected in the plasma at temperature T. From ( Jed06bl ) . 

Ec well approximated by fjPro95l : iKawQSl ) 



Kn 



7_ 

Ex 



for E^f < Ex 



for Ex <E^ < E, 



(152) 



c 



where Kq = Eq/IExI'^ + \n{Ec/ Ex)]) is a normalization constant such that 
the total energy in 7-rays below Ec eq uals th e total energy Eq injected. One 
has Ex ~ 0.0264i?c according to Ref. (jPro95l ). or Ex ~ O.OSEc according to 
( 1Zdz89l ). A more accurate calculation of the evolution of this spectrum should 
include interactions of these "break-out" photons via photon-photon scattering 
7 + 7cMB -^7 + 7, (mainly redistributing the energy of energetic 7-rays right 
below energy Ec), Bethe-Heitler pair production 7-|-p(^He) — > p('^He)-|-e--|-e"'', 
Compton scattering off thermal electrons 7 + e" — > 7 -|- e" (with the produced 
energetic e~ in turn inducing inverse Compton scattering and thus further low- 
energy 7's) and, of course, nuclear photodisintegration, which directly affects 
BBN. 



In Figure 28 we show the amounts of destroyed nuclei per TeV of injected elec- 
tromagnetic energy, Eq, at a temperature T. The sharp rises at characteristic 
temperatures are due to the pair production threshold effect for energetic 7, 
which depend on the binding energy. In Figure 29 we show a similar plot for 
the nuclei produced as a result of the dissociations. It is evident that, although 
a minor amount of starts to be produced from spallation of ^He already 
at T ~ 1 keV, a major secondary production of ^H, and ^He happens at 
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Fig. 29. Number of produced nuclei per TeV of electromagnetically interacting en- 
ergy injected in t he plasm a at temperature T, due to photodisintegration and fusion 
reactions. From ( JedOGbl ). 



T < 0.5 keV, when ^He photodisintegration becomes relevant. In this range, 
not only does ^H, and ^He production over-compensate their photodisin- 
tegration (which takes place at a slower rate) but a significant synthesis of 
^Li is induced by the non-thermally produced '^H and "^He, with a subleading 
contribution from direct photodisintegration of ^Li and ^Be. For a compila- 
tion of the reactions involved in the casca de nucleo synthesis calculations, we 
address the reader to the excellent review f|jed06bl l Also, some relevant new 
nu clear as trophysics data and their implications for BBN have been presented 
in flKusOSh . 



In Figure 30 we show a typical result for the excluded regions in the plan e 
Cx-Tx (where refers to a time t -C tx), in this case taken from (ICyb03al ). 
In the left panel only deuterium constraints are considered, in the right panel 
constraints from all elements are included. The light blue regions are excluded 
due to underproduction of deuterium: in the left shoulder at relatively small 
Tx this is due to the fact that direct photodestruction of dominates. For 
r > 10^ s "^He destruction is important, and leads typically to over-production 
of ^H, unless the injected energy is so high that in turn even this secondary deu- 
terium is destroyed (upper- right corner of the plot). In this parameter range, 
typically even ^He and ^Li are destroyed to a level inconsistent with obser- 
vations, but these additional constraints basically overlap with the deuterium 
ones. On the other hand, if one considers ^Li production from non-thermally 
produced A = 3 nuclei, one can also disfavor the red region shown in the 
right panel due to over-production of ^Li, which is however less robust due to 
the likely reprocessing of this fragile isotope in the observed stellar systems. 
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Fig. 30. Regions of the parameter space Cx-tx excluded by: left panel: overpro- 
duction (dark blue) or under-production (light blue); right panel: including also 
the constraints from overproduction of ^Li ( dark red ), shown on the top of the ^Li 



(yellow) and ^He (green) constraints. Prom ( Cvb03al ) . 



Note that if the metastable particle is the progenitor of the DM candidate, it 
must fulfill the condition rrixUx > AmpUB, i.e. a typical value for (x fulfills 
^ 2 X 10~^ GeV. Larger values imply a very large mass difference between 
the X particle and the DM candidate, smaller values refer to metastable par- 
ticles which account at most for a sub-leading fraction of the DM today. So, 
a way to summarize previous constraints is to say that BBN excludes electro- 
magnetically decaying particles as progenitors of DM when their lifetime is 
longer than tx > few x 10^ s. 



9.1.2 Including hadronic channels 

Generally, in the decay of the X particles one expects both hadronic and 
electromagnetic channels. The treatment of cascades in presence of hadrons is 
significantly more involved, since it introduces more particles, more timescales 
and more processes to take into account. The flow-chart in Figure 31 summa- 
rizes the decay scheme and physical processes involved. Note that a hadronic 
branching ratio unavoidably leads to secondary electromagnetic showers, which 
affect BBN along the lines described in Section 9.1.1. To have a flrst quali- 
tative understanding of the effects of hadronic particles, it is worth recalling 
that different hadronic species interact in the plasma via a few reactions, which 
assume different relevance at different times: 
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Fig. 31. Flow-chart of the decay effects on BBN. From (jKawOSal ). 



• mesons, mostly tt^ and kaons (with rest frame lifetime in the 10~^ s range), 
have only an effect at times t ~ 1 — 10 s, when ordinary weak interactions are 
not efficient anymore, but they still have time to interact before decaying. 
They mostly act by enhancing the n/p ratio and thus the final value of Yp. 

• antinucleons, by the preferred tendency to annihilate onto protons, have a 
similar final effect of increasing n/p. Compared with mesons, they also have 
an additional peculiar effect at later times (t ~10^ s) by annihilating onto 
Yp and leaving ^H, ^H, and ^He among secondaries. 

• nucleons: at early times, nucleons thermalize via electromagnetic processes: 
magnetic moment scattering off for neutrons. Coulomb stopping off 
and Thomson scattering off thermal photons for protons. However, at late 
times other energy loss mechanisms start to dominate for high energy nucle- 
ons, namely nucleon-nucleon collision and nuclear spallation reactions. Due 
to the different electric charge, these nuclear processes are already dominant 
for neutrons at t > 200 s, while for protons only at t > 10^ s. When they 
are effective, a cascade nucleosynthesis can take place: each nucleon-nucleon 
scattering will produce another energetic nucleon (a single 100 GeV nucleon 
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Fig. 32. Light element nuclei yields for a particle of mass mx=l TeV decaying as 
X — > qq as a f unction of the temperature T of the cascade injection in the plasma. 
From (|jed06bh . 



can produce several tens of 10 MeV nucleons) and their effect of spallation 
over '^He will produce many ^H, ^H, and ^He nuclei. The total effect will be 
more efficient than for antinucleons: a single 100 GeV antimatter particle 
has a much shorter mean free path before annihilating with a ^He nucleus 
and produce the mentioned secondaries, and nucleon induced production is 
therefore much more efficient. Non-thermal nucleon injection lead to an in- 
creased Yp abundance at t <200 s, increased abundance at 200 s< t <10^ 
s, or decreased ^Li abundance at t falO^s. Spallation of ^He to produce ^H, 
■^H, and '^He may have as a secondary effect the synthesis of ®Li. 

To illustrate this point, in Figure 32, we report the light element yields pro- 
duced by the decay X — > gg of a single metastable particle with mass mx=l 
TeV, as a function of the photon temperature at which the injection takes 
place. It is interesting that initially, after hadronization of the quark-ant iquark 
state, on average only 1.56 neutrons result; all others are created at T <90 
keV by the thermalization of injected neutrons and protons due to inelas- 
tic nucleon-nucleon scattering and ^He spallation. Similarly, all the ^H, ^H, 
and "^He nuclei are due to ^He spallation processes and n-p nonthermal fu- 
sion reactions (for ^H) induced by the thermalization of the injected energetic 
nucleons. 



For a more quantitative analysis, a numerical treatment is required. Tech- 
nically, the problem of solving cascade nucleosynthesis is highly non-trivial, 
especially in presence of significant hadronic branching ratios. One major dif- 
ficulty is that one cannot study the non-thermal effects independently of the 
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earlier standard BBN stage, as it is possible (at least as first approximation) 
in the case of electromagnetic cascades. Although already in the 80 's several 
authors have estimated the effects of injecting a non-thermal population of 
particles in the plasma at the BBN epoch, only recently several authors have 
followed self- consistently the reactions taking place during their thermaliza- 
tion and coupled them to the simul taneously ongoing thermal nuclear network 
f lKawOSbl : Ijed04bl : ICvbOd : Ijed06bl l The latest treatments typically adopt a 



Monte Carlo technique to calculate the interaction probability of the particle 
shower, and codes such as PYTHIA are employed to determine the energy 
distribution of the shower particles given the initial branching ratios in the 
decay. Nucleon energy losses must be taken into account as well. For a more 
detailed overview of the scheme and techniques of the numerical treatment, 
see e.g. flKaw05al : [jed06bl l 



The resulting constraints can be presented similarly to the ones for e.m. de- 
cays, still it is crucial to specify not only the injected energy, but also the 
hadronic branching rat io, Bt- For example, in Figure 33 we report the BBN 
constraints derived in f|jed06bl ) for the decay of particle of mass mx=l TeV 



and hadronic branching ratio, 5/^=3. 33x10"^. The quantity VLxh? is the con- 
tribution that neutral particles would have given to the total energy density to- 
day, would they have not decayed, and is proportional to the parameter C,x in- 
troduced in the previous Section {VLx = {Mxn\h?) / pcr = {n^h'^ I Pcr)Cx — 
3.9 X lO^Cx/GeV). Colored regions are excluded and correspond to constraints 
imposed by upper limit for Yp (orange area), upper limit on (blue), upper 
limit on ^He/^H(red), and lower limit on ''Li (light blue). The yellow region 
violates the less conservative bound from ^Li/''Li. References for observational 
abundances used to get these constraints, which partly differ fror n the on es 
we have compiled in Section 4, can be found in the original paper f|jed06bl ). 



As a consequence of the different mechanisms dominating at different times, 
the most stringent constraints are given by: the overproduction of ^He at early 
times (r^ <10^s), overproduction of for 10^ s< < 10^ s, overproduction 
of ^Li for 10'^ s< < 10^ s, and an overproduction of the ^He /^H ratio for 
Tx >10^s. Figure 34 illustrates how the constraints depend on the branching 
ratio, -B/i, this time for a particle mass, m^=100 GeV, with the allowed region 
below the lines. In the limit — ^ 0, one recovers bounds due to the e.m. 
channels, while even for Bh ^ 1% metastable particles with a relic abundance 
comparable to today's DM one are excluded if tx > 10^ s. Of course, one can 
turn the argument around and explore the possibility that decays may explain 
some discrepancies existing between the SBBN model and observation. In 
particular, both "lithium problems" mentioned in Section 4.5 may be solved 
by hadronically de caying particles with lifetime ~ 2000 s (see e.g. Figure 2 



in (jJedOGal : iJedOSa] )). It is worth mentioning that in some of these scenarios 
"low-energy" observables (like the properties of gravitino dark matter) would 
be linked to very high energy physical scales, as the reheating temperature 
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Fig. 33. Conservative BBN constraints on the abundance of relic decaying neutral 
particles as a function of their lifetime; particle mass is m-^=l TeV and the hadronic 
branching ratio i?/j=3.33x 10~^. ^^^^^ contribution neutral particles would 

have given to the total relic density today, would they have not decayed. Colored 
regions are excluded and correspond to constraints imposed by observations (see 
text). From (|jed06bl ). 




Fig. 34. Constraints on the relic abundance of neutral decaying particles as a func- 
tion of their decay time, for a particle mass mj(.=100GeV and varying branching 
ratio, B^. Numbers in the picture refer to the solid line above and stand for the 
corresponding log^g ^h- Dotted lines are the actual constraints if one considers the 
less conservative ^Li/^Li. From ( Jed06bl ) . 
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Fig. 35. The abundances of ^Li as produced in scenarios with hadronic cascades in- 
jected by neutrahno annihilation, for different channel s, as labeled. The annihilation 
rate is taken to be {av)= 10" 



cm /s. From (Ijed04bh . 



and the leptogenesis scale, see for example (ICerOGi : lKan07al : ISteOSi ). 



In case of an annihilation mode, the major peculiarity is that in order to 
have the same injection rate of a decay at an epoch Zinj, a very large anni- 
hilation rate at epochs z > Zinj is required. Typically, in order to produce 
significant modifications of the light element yields, often an unphysically 
large (av) is required. Stronger bounds follow then from non-BBN consid- 
er at ions^jonless a relatively light particle {mx < few GeV) is considered (see 
e.g. f McDOObI )). A possible exception is provided by the ^Li nuclide. It has 
been shown that, for typical electroweak-scale WIMP masses, the primor- 
dial ^Li yield from DM annihilation exceeds the SBBN abundance for any 
{av) > 10-27 cmVs fljed04al : [jedoibh . 



This is illustrated in Figure 35: interestingly for the annihilation rate required 
to produce the correct WIMP DM thermal relic abundance, an O(IOO) GeV 
DM particle produces an amount of ^Li in the same range required to explain 
observational claims of ^Li. 



9.2 Catalyzed BBN 



The constraints derived in the previous Section assume that, before decay- 
ing, the meta-stable particle X is "inert". While DM can not be charged 
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or strongly interacting, there is no a priori reason why its parent particle 
should bring no electric or strong charge. Whenever this happens, the possi- 
bility arises that X particles may form bound systems with baryons, altering 
the nuclear reaction pattern and thus the yield of light elements. This sce- 
nario is currently known as catalyzed BBN (CBBN). The cosmological role 
of ch arged inassive particles (CHAMPs) was already considered in the late 
80's flDeR89l : IPimSol : iRafSOh . but the influence of bound states in BBN has 



only been fully appreciated recently. In 2006, within a few weeks three pa - 
pers appeared pointing out the importance of CBBN (jPosOGi : iKohOTI : iKapOGi ) 
and identifying the main physical mechanisms responsible for the alteration 
in the nuclear network. In the last two years, several articles have followed, 
clarifying the physical ingredients regulating this complex scenario, including 
refin ement in the calculation of cata l yzed r e actions, late-tirn e nucleosynthesis, 
etc. (ICvbOel: iHamOTl: iBirOTi : iKawOTi : IjitOTi : IJedOSal : ICumOTi : IJedOSbl : iKamOS : 



Kus07al: lKus07bl : lJit08l ). A quite complete review of the physics can be found in 



f lJedOSal ) , which we shall mainly follow for the present summary. We limit our- 
selves to consider singly charged CHAMPs, although it has been argued that 
negative doubly charged and stable CHAMPs bou nd to ^H e"*"^ may be viable 
as DM candidates in walking technicolor theories (IKhlOSal ). It is worth noting 
that, at least qualitatively, one expects similar catalytic mechanisms if the 
particle X is strongly interacting, rather than being electrically charged. One 
physically motivated scen ario of this kind is the long-lived gluino in split-SUSY 
( lArkOSal : lGiu04l : lArkOSbl ) . Some calculations of the primordial nucleosynthesis 
in presence o f massive, stro ngly interacting particles have been performed in 
the past (see (jPicSOl : |Pla95l )). where bound states with "ordinary" nuclei were 
considered. The main difficulty with this scenario is that the nuclear physics 
of such bound states is very hard to treat reliably. For ex ample, analogies 
with toy-models used to describe hypernuclei are employed flMohQSl ) and the 
description of the bound systems is at best parametric. An accurate treat- 
ment is made highly non-trivial by the effects of non-perturbative physics, 
and we shall not consider them further. We want to remark, however, that the 
BBN cascade bounds that are sometimes considered in the literature for such 
particles (see flArv05h ) may be altered, since catalytic effects are completely 
neglected. 



9.2.1 Early time CBBN: formation of bound states and catalysis 

A negatively charged, long lived particle X (or CHAMP) with mass Mx ^ 
nip would form bound-states with nuclei and alter the network of reactions 
leading to the synthesis of light elements. Compared with a generic nucleus 
A of charge Z and mass m^, its corresponding bound state with a CHAMP, 
AX, has a mass higher by ~ mx, a charge lower by one unit {Z — 1), and 
it is characterized by an AX binding energy given by Eb ~ Z"^ amA/'^ in 
the limit <^ mx- Another interesting quantity is the photo-dissociation 
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temperature Tph of the AX system, defined as the temperature at which the 
photodissociation rate, TphiT), for the bound nucleus becomes smaller than 
the Hubble rate, H(T). Roughly speaking, below Tph the AX system is stable 
against photodestruction. In Table 16 we report binding energies in the Bohr- 
like atom approximation (corrections are of order ~ 10 % for ^ HeX, up to 



50% for ^BeX) and the Tph, as originally calculated in (jPosOGi ) . A precise 
calculation of the fractional abundance of bound nuclei requires solving the 
Boltzmann equation, including all relevant reactions for the AX state. The 
inadequacy of simple approaches based on the Saha equation t o deter mine 



the evolution of the AX abundance was realized immediately in (jPosOGi ) and 



further analyzed e.g. in (jKoh07 : Icvboel i However, even the latter analysis 



had to rely on strong approximations on the reaction rates for bound nuclei 
(that the bound state nucleus would be destroyed in the interaction, and that 
standard BBN is over at the time CBBN takes place, etc .). Figu re 36, which 



we take from the more recent and accurate analysis in (jJedOSal ). shows the 
evolution of bound state fractions f\ = uax / n^A of nuclei A bound to X as a 
function of temperature T, for a model with mx = 100 GeV and Vlxh? = 0.1. 
It confirms that the behavior of is significantly different than that expected 
from simple estimates by the Saha equation due to nuclear destruction of 
bound states and slow recombination rates (with respect to the Hubble time) . 
This is particularly relevant in fi^^ due to the ^LiX(p, X) 2'^He reaction. 

Once CHAMP-nuclei states XA are formed in the plasma, each nuclear reac- 
tion will have its CHAMP homologue: 



standard BBN 
CBBN 



A + Ai^ 
AX + Ai 



A2 + A3, 
-^A2 + As 



X . 



(153) 



At first sight, the main advantage of CBBN reactions is the smaller Coulomb 
barrier. Howev er, the re is a more subtle effect — whose importance was already 
recognized in fjPos06l ) — which acts on homologues of radiative captures , i.e. 
of reactions of the kind A{Ai,'~f)A2. If we denote by the wavelength of the 
emitted photon, in general of the order of 100 fm, electric dipole (El) reaction 
rates scale as A~^, whereas electric quadrupole (E2) ones scale as A~^. The 
introduction of the photonless state in the CHAMP-mediated reaction replaces 
A^ with the Bohr radius of the bound system, approximately 5 fm for ^He. 

Table 16 

Binding energies Ef, in the Bohr approximation \Ei,\ ~ a mA/2 and photo- 
dissociation decoupling temperatures Tph in keV (as calculated in ( PosOd )) for exotic 
bound states AX. 



bound state 


pX 


2HX 




^HeX 


^HeX 


■^LiX 


^BeX 


^BeX 


\Eb\ 


25 


50 


75 


299 


397 


1566 


2787 


3178 


Tph 


0.6 


1.2 


1.8 


6.3 


8.2 


21 


32 


34 
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Fig. 36. Bound state fractions f\ = riAx/n^A nuclei A bound to CHAMP X~ 
as a function of temperature T, for a model with Mx = 100 GeV and Qxh^ = 0.1. 
Shown are for ^Be solid (red), ''Li long-da shed (gr een), ^Li short -dashed (blue), 
and ^He dotted (purple), respectively. From (jjedOSal ). 



Table 17 

Dipole amplitude, Q-values, and catalyzed S -factor e nhancement in the cross-section 
for relevant (0,7) rea ctions. M odified from ( Cvbod ) . taking into account the values 
recently reported by ( KamOSi ) . 





EM 


iV(a,7)C 


NX{B,C)X- 


Enhancement 


Reaction 


Transition 


QsPN (MeV) 


QcBBN (MeV) 


5'cbbn/5'sbbn 


2H(a,7)6Li 


E2 


1.474 


1.124 


~ 10^ 


^R{a,jy Li 


El 


2.467 


2.117 


~30 


3He(a,7)^Be 


El 


1.587 


1.237 


~30 



In terms of the usual parameterization of low-energy nuclear reaction (see 
Eq. (122)), the former Coulomb effect enters the exponential barrier-factor, the 
latter — when present — influences the S'-factor, which can be indeed enhanced 
by orders of magnitude, as reported in Table 17. Another, usually sub- leading 
effect in the direction of decreasing 5* is a decrease of phase-space due to the 
reduction of the reaction Q- value. 

The most important alteration to the standard scenario in catalyzed BBN is 
mainly due to the enhancement of the single ^Li producing process. 



standard BBN : + ^ ^Li + 7; Q = 1.47 MeV 



(154) 
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replaced by the process 



CBBN : ^HeX 



'Li + X; Q ~ 1.13 MeV 



(155) 



which in ( IPosOGI ) was identified to be very effective even for a small fraction 
of X particles bound with nuclei. The usual BBN process of Eq. (154) is in- 
deed only allowed at the quadrupole level (due to the almost identical mass 
to charge ratio of and Yp), which is the reason for the very small value of 
®Li/''Li, as already discussed in Section 4. Anot her possible path to enhance 
the ®Li yield in CBBN was proposed in (IKapOGl ). who noticed that the decay 
of X when still in a bound state with Yp could result in a break-up of the 
Yp nucleus, producing '^He and that would eventually fuse into ®Li when 
reacting with Yp. However, the possibility that this would happen is estimated 
to be very low and the ^B.j'^B.e X, ^Li)X appe ars to be in all cases still dom- 
inating the production of ^Li (IKapOGl : Icvboel i The enhancement of the ^Li 
yield in CBBN due to the process in Eq. (155) has been confirmed by all the 
publis hed anal ysis as the most remarkable effect of CBBN and in particu- 
lar by flKamOSl ) . who performed an accurate quantum three-body calculation 
of the cross-section of the reactions involved in CBBN. The observational 
hints of a plateau in ^Li at a value well above standard BBN predictions, 
as well as the persisting discrepancy between ''Li observations in the Spite 
Plateau and the (apparently overproduced) ^Li yield has motivated several 
authors to explore the CBBN scenario further, trying to also explai n the ^L i 
"problem". A mechanism to address this issue was pointed out in (jKohOTI ). 
Since significant fractions of ^Li and (mostly) ''Be are in bound states with 
CHAMPs, ^Be can be depleted by the enhancement of the CBBN a nalogues 
of ^Li(p,a)^He, ^Be(n,p)^Li, and ''Be(n,Q;)^He. The authors of Ref. (ICumOTi ) 
performed a CBBN analysis also adding the effect of X decay cascades. They 
concluded that in presence of strong showers from decaying relic particles, 
bound-state effects on nucleosynthesis are negligible, and both Li problems 
are solved (i f at all ) in a way very similar to the cascade BBN case in absence 
of catalysis. fiBirOTh proposed a more elaborate solution of the ^Li problem: the 
'^BeX{p, ■j)^BX, and the subsequent beta-decay of ^B — > ^Be-|-e"'"-|-t'e would 
deplete the final ^Li abundance, with little consequence on Yp. The ^BeX(p, 7) 
reaction would happen through a shifting of the resonance as an effect of the 
X presence, which would lead to a huge rate enhancement. 



To a large extent, the reason why settling these issues is far from trivial is 
that significant uncertainties remain in the estimates of binding energies and 
CBBN reaction rates, due to the use of very simplified nuclear models and of 
the Bo rn appro ximation. An account of the situation has been given in Section 
III of (|Jed08a] l. which we address the reader to for further technical details. 
(jJedOSal ) also contains the most systematic and up-to-date analysis of CBBN, 
including calculations of rates of AX recombination-photodisintegration and 
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CBBN analogues of BBN nuclear reaction rates. As long as the range 10 keV > 
T > 0.8 keV is concerned, the author confirms that the most relevant role is 
played by re action (1 55), even when using the detailed calculation for its rate 
obtained in flHamOTjl, Moreover, in this regime only nine reactions (reported 
in Table 11 of ( iJedOSal )) are sufficient to describe the physics of CBBN. The 
evolution of light nuclide abundances in presence of bound-state reactions is 
shown in Figure 37, where the enhancement of ^Li abundance is clearly visible 
(note however that the role of X decays has been "switched off"). 
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Fig. 37. Limits on the primordial CHAMP to entropy ratio = nx/s for 
CHAMPs with intermediate lifetimes. Shown are constraint Unes for CHAMPs of 
mass mx=100 GeV and a variety of hadronic branching ratios, = 10~^ — 1, as 
labelled in the Figure. Solid (red) lines correspond to the conservative limi t ^Li/^L i 
<0.66, whereas dashed (blue) lines correspond to ^Li/^Li <0.1. From ( JedOGbl ) . 
updated (courtesy of K. Jedamzik). 



9.2.2 Late time CBBN: HX bound states, CHAMP-exchange, and decays 



The authors of (jKoh07l ) suggested the possibility that the formation of bound 
state nuclei of CHAMPs with '^H, and p at very late times might induce 
the suppression of synthesized ^Li. Bound states of Z=l nuclei with X form at 
1-2 keV temperatures, see Table 16. These bound states behave essentially as 
"long-liy ed" neutrons, which can dissociate Li and ^Be Coulomb-unsuppressed. 



f lJedOSal ) carried on an extensive analysis of these effects, making use of reac- 
tion rates derived in the Born approximation, identifying nineteen reactions 
as relevant when late time CBBN effects are taken into account. These early 
results indicated that bound states of CHAMPs with Z = 1 WOULD in- 
duce at late times the destruction of most of the synthesized ^Li and some 
^Li. Initially, it appeared that the stringent constraints initially put on the 
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0.001 



Fig. 38. Evolution of light-element number ratios ^Be/H (solid - red), ^Li/H (long- 
dashed - green), ^Li/H (short-dashed - blue), and ^H/H (dotted - purple), for a 
CHAMP model with Mx = 100 GeV, nxh2 = 0.01, and tx = 10^°s. Neither ef- 
fects due to electromagnetic and hadronic energy release duri ng CHA MP decay 
nor charge exchange effects have been taken into account. From ( JedOSal ). updated 
(courtesy of K. Jedamzik). 



abundance of CHAMPs, e.g. in (IPosOGI ) and (iKawOTI ) loosen significantly. To 
add ano ther lay er of complication, an additional late time effect was pointed 
out in fj JedOSal ). which somewhat compensates the previous one. It is due 
to exothermal transfer of CHAMPs from Z = 1 nuclei into heavier nuclei 
(CHAMP-exchange reactions), i.e. HX('^He, H)^HeX. By lowering the abun- 
dance of neutral states B.X, more of the ^Li and ''Be produced can sur vive. 
The detailed study of bound state reaction rates performed by (IKamOSl ) has 
finally shown that CHAMP-exchange reactions dominate the late time CBBN, 
thus "protecting" the abundances of Li and ^Be produced at earlier times, as 
it can be seen in Figure 38. The late-time drop in ''Be is due to electronic 
capture decays. 



9.2.3 ^Be from CBBN 

Finally, we want to point out here that modifications to the yields of heavier 
than ''Li elements can take place in CBBN. As we have summarized in Sec- 
tion 3.4, one robust prediction of the standard BBN is the absence of sizeable 
yields of A >7 elements. Roughly speaking, this is due to the lack of stable 
A=8 elements, and to the inefficiency of the 3 a — C, which would respec- 
tively allow a slow light element chain to produce heavier elements or "bridge" 
it, as it happens in stars. The very short lifetime of *Be is a problem that can 
be overcome in CBBN, where sizeable amounts of meta-stable bound states 
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^BeX can be c reated through the mechanism ^HeX(a,7)^BeX at T<30 keV, 
as showed by (jPosOTI ). Once this bound state has been produced, the neu- 
tron capture reaction ^BeX(n,^Be)X can take place and ^Be be efficiently 
produced in CBBN, whereas the analogue SBBN mechanism is suppressed as 
a consequence of the lifetime of ^Be, shorter than a femtosecond. Although 
the absolute abundance of ^Be produced b y this mechanism is sensitive to 
the CHAMP abundance in the plasma, Yx, fjPosOTl ) argued that an enhanced 
production of ^Be and ^Li are peculiar signatures of CBBN, and found that 
their ratio is independent of Yx- The author derives the conclusion that the 
primordial ratio ^Be/^Li ~10~^ should therefore be a "signat ure" of CBBN, 
as it dramatically differs from ^Be/^Li ~10"^ ( ^Be/H~ 10~^^, (llocOTl )) which 
one obtains in SBBN. In their detailed study, (IKamOSl ) have c halleng ed this 
result by arguing that the rate estimate of ^BeX(n,^Be)X by ( PosOTI ) is too 
high; although strong conclusions would require a careful examination of all 
the nuclear effects in volved, they hint toward a much lower efficiency of the 
process proposed by (jPosOTl ). While this issue is yet to be clarified, the kind 
of constraints to supersymmetric scenarios that might follo w from a pecu liar 
^Be production in CBBN has preliminarily been studied in (jPraOSi : iBaiOSi ). 



9.2.4 Constraints in CHAMP parameter space 

If the lifetime of the CHAMP is tx <3x 10^ s, i.e. before that bound states can 
form, clearly its only effect on BBN is equivalent to injecting electromagneti- 
cally and hadronically interacting particles into the plasma. For intermediate 
lifetimes, 3x10^ s< tx <5x10^s, the main novel constraint is due to the pos- 
sible overproduction of ^Li via ^HeX(^H, X)^Li. Fortunately, as this reaction 
i s now k nown within a factor three in light of the dedicated calculations of 



(iHamOTI ). this range of lifetimes is relatively well constrained. However the 
decay of a CHAMP at this time induces again a cascade nucleosynthesis; in 
hadronic cascades from a neutral particle decaying at the same time, it is 
well known that ^Li overproduction is the main effect. Thus, unexpectedly, 
for hadronic branching ratios Bh >0.01 (and e.g. Mx=l TeV) the effects 
of a cha rged particle decay do not differ much from those of a neutral one 
(jJedOSbl ). The results of the previous Section apply since the ^Li produced by 



the hadronic shower induced reactions dominates over the effect of catalyzed 
^Li production. Only for Bh < 10~^ and for sufficiently large decay times tx 
the CBBN mechanism provides new bounds. 

Finally, for longer lifetimes, tx >5x 10^ s, th e conservative limits on charged 
decaying particles initially derived in ( JedOSbl ) appeared to be no stronger than 



those on neutral particles, due to the large uncertainties in nucleosynthesis 
at T<3 keV. However, the nuclear rate calculations for CBBN presented in 



(IKamOSl ) strongly reduce those uncertainti es. A recent application of these 



updated calculations has been performed in (IBaiOSf l with an analysis restricted 
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to a minimal SUGRA model with heavy gravitino, which we address for further 
details. 



10 Conclusions 

In this review we have reported the current status of BBN, focusing in the 
first part on precision calculations possible in the standard scenario, which 
provide a tool for current cosmological framework, in the second part on the 
constraints to new physics, which become particularly important in the forth- 
coming LHC era. The "classical parameter" constrained by BBN is the baryon 
to photon ratio, 77, or equivalently the baryon abundance, VLsh?- At present, 
the constraint is dominated by the deuterium determination, and we find 
VLsh"^ = 0.021 ± 0.001(1 a). This determination is consistent with the up- 
per limit on primordial ^He/H (which provides a lower limit to 77), as well 
as with the range selected by ^He determinations, which however provides a 
constraint almost one order of magnitude weaker. The agreement within 2 
a with the WMAP determination, fi^/i^ = 0.02273 ± 0.00062, represents a 
remarkable success of the Standard Cosmological Model. On the other hand, 
using this value as an input, a factor > 3 discrepancy remains with ^Li deter- 
minations, which can hardly be reconciled even accounting for a conservative 
error budget in both observations and nuclear inputs. Even more puzzling are 
some detections of traces of ^Li at a level far above the one expected from the 
Standard BBN. If the observational determinations are solid, both nuclides 
indicate that either their present observations do not reflect their primordial 
values, and should thus be discarded for cosmological purposes, or that the 
early cosmology is more complicated and exciting than the Standard BBN 
lore. Neither a non-standard number of massless degrees of freedom in the 
plasma (parameterized via A^eff) or a lepton asymmetry (all asymmetries 
assumed equal) can reconcile the discrepancy. Current bounds on both quan- 
tities come basically from the ^He measurement, Ncr = 3.2 ± 0.4 (la) and 
= -0.008 ±0.013 (la). 

On the other hand, other exotic proposals have been invoked to reconcile this 
discrepancy. Typically they involve massive meta-stable particles with weak 
scale interactions, which should be soon produced at the LHC. In Super- 
symmetric scenarios, long-lived particles are possible whenever the Next to 
Lightest Supersymmetric Particle (NLSP) decays into the Lightest Supersym- 
metric Particle (LSP) are gravity-mediated, or "disfavored" by phase space 
arguments, with a modest mass splitting between NLSP and LSP. Cases fre- 
quently considered in the recent literature are neutralino — > gravitino decays, 
for example, or stau — > gravitino. The phenomenology associated with the 
catalysis of reactions due to bound states of charged particles (as the stau) 
with ordinary nuclei is a particularly new topic in recent investigations. Also, 
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the importance of a possible primordial origin of the ^Li measured in a few 
systems of the ^Li plateau has been recognized: first, the bounds in param- 
eter sp a ce tighten sig n ificantly if lithium constraints are used, especially ^Li 



(jHol96l : iJedOOl : lHol99l : lJed04a] ): second, because these exotic BBN scenarios 
may accommodate for a cosmological origin for ^Li while solving the ^Li ex- 
cess problem as well, their phenomenology is very appealing. Although these 
links among primordial nucleosynthesis, dark matter, and perhaps SUSY phe- 
nomenology are quite fascinating, it is worth stressing that BBN bounds on 
cascade decays or annihilations of massive particles apply well beyond the re- 
stricted class of SUSY-inspired models. For example, the electromagnetic cas- 
cades following heavy sterile neutrino decays are constrained by these kinds 
of arguments, as well as decays of massive pseudo Nambu-Goldstone bosons. 



as considered in (jMasQTl : lMas04l ) . 



There are two directions along which we can expect the BBN field to develop 
in the future. On one hand, BBN is an important tool for precision cosmology, 
especially if its priors are used in combination with other cosmological observ- 
ables. Already BBN provides the best bounds on parameters as A'cff and .^e 
(and bounds on rj comparable to the CMB); yet, since theoretical uncertainties 
are at the moment well below observational ones, there is surely room to refine 
its power, provided that significantly greater efforts are devoted to determine 
light element abundances, and in particular Yp. It is instr uctive in this sense 



to look back to what S. Sarkar wrote in his review (jSar96l ) thirteen years ago: 



Thousands of person years of effort have been invested in obtaining the pre- 
cise parameters of the resonance in — e~ collisions, which measures 
the number of light neutrino species (and other particles) which couple to the 
Z". In comparison, a modest amount of work has been done, by a few small 
teams, on measuring the primordial light element abundances, which provide 
a complementary check of this number as well as a probe of new superweakly 
interacting particles which do not couple to the Z^. 

Despite the improvements reported in this article, we feel that unfortunately 
insufficient attention has been devoted to this problem, if compared to other 
areas of observational cosmology. In particular, the ^He determination is still 
plagued by systematic uncertainties. Although their importance has been re- 
cently recognized and assessed more carefully, the fact that this reanalysis 
was triggered after the independent determination of rj from CMB (and its 
agreement with the "low deuterium determinations" in QSO spectra) shows 
that there is still a long way to go towards a precision era for primordial ele- 
ments. On the other hand, a significant improvement has taken place in assess- 
ing and reducing theoretical uncertainties, mostly related to nuclear reaction 
data. BBN has benefit from a wealth of new nuclear astrophysics measure- 
ments at low energies and covering large dynamical ranges. Given the much 
larger observational uncertainties, in this sector an effort in reassessing the 
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systematic errors in older datasets might be more useful in reducing remain- 
ing discrepancies in the nuclear rates error budget. This is in particular the 
case for reactions involving ^Be. 

The other direction of development follows from the interplay with Lab exper- 
iments. Neutrinos have reserved many surprises, and it is not excluded that 
exotic properties may show up in future experiments with important impli- 
cations for BBN, as we illustrated in Section 6. However, it is in particular 
from LHC that one expects a better understanding of high energy scales, and 
thus of the cosmology at earlier times and higher temperatures. Most theories 
that go beyond the Standard Model of Particle Physics require new states to 
appear at or above the electroweak scale and, as already reported, they might 
have implications for the phenomenology at the BBN epoch. If the LHC should 
provide indication for the existence of the SMPP Higgs and nothing else, there 
will be no natural scale to explore. In this case, albeit sad, BBN and other 
cosmological tools might be the only practical means to explore very high 
energy phenomena leaving their imprint on the cosmos. One example treated 
here is the effect of variations of fundamental "constants" on cosmological 
time-scales that emerge in extra dimensional scenarios, possibly embedded in 
grand unified theories or string theories. If, as hopefully more likely, the LHC 
will reveal new dynamics above the electroweak scale, we might be able to 
infer from the empirical evidence the presence of cosmological effects before 
the BBN epoch. A new Standard Cosmological Model would emerge as well, 
perhaps making the BBN one more step in the ladder back to the Big Bang, 
rather than the first one. 
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